Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 948,211 G→A P18L (CCG→CTG)  unknown ← 3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1948,2110GA96.3% 137.7 / NA 56P18L (CCG→CTG) unknown3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (17/35);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCGGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCCATTCATTACGCTGCTCTCG  >  1/947977‑948457
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
aCGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAgg                                                                                                                                                                                                                                          <  1:276103/250‑1 (MQ=255)
aCGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGc                                                                                                                                                                                                                                         <  1:82033/251‑1 (MQ=255)
            tcgtcgGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAAATCGTTATCGGGGGGGATAATGACGCGGGCGGCACCGGGTTGTTTGGCTTGTTTTAAAAAGTGTTGAATCAACGGCTTGGCGGCAATATC‑GGGTAA‑GGGGTGGCCCCGTAAAGGGCG‑GGAGGCGGAAcggggggg                                                                                                                                                                                                                             >  1:245738/1‑251 (MQ=255)
                       caCCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑c                                                                                                                                                                                                                   <  1:541812/249‑1 (MQ=255)
                              ttttCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATtg                                                                                                                                                                                                          <  1:241221/251‑1 (MQ=255)
                                 tCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTgtgt                                                                                                                                                                                                       <  1:270342/251‑1 (MQ=255)
                                      cacTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCatat                                                                                                                                                                                                  <  1:183609/251‑1 (MQ=255)
                                           ccGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATattaatt                                                                                                                                                                                             <  1:197226/251‑1 (MQ=255)
                                                         gTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTTTCGGGGGCGGTAATGACGCGGGCGGCACCGGGTTGTTTTGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCCATATC‑GGCTAAGGGGTTTGCCCGG‑AAACGGCTGGGGGGGGAAAGGGGGGGGAAGGA‑CA‑CGGGGGATTGTGGCATATTTATTATTCTGTTCTTGAt                                                                                                                                                                               >  1:293855/1‑251 (MQ=255)
                                                          tGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGAt                                                                                                                                                                               <  1:566682/250‑1 (MQ=255)
                                                                aaTTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACgg                                                                                                                                                                        <  1:78674/251‑1 (MQ=255)
                                                                       gTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTgtgt                                                                                                                                                                                                       <  1:851137/213‑1 (MQ=255)
                                                                       gTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTgcgc                                                                                                                                                                  <  1:314531/250‑1 (MQ=255)
                                                                                 gTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTTTATGTGGCGATAATGACGCTGGCGGCACCCGGTTGTTTTGCTTGTTTATAAAAGTGTTGAATCATCCGCTTCCCCGCCATATCGGGTTAG‑GGGTTTTCCCGTTAA‑CGGCT‑GGCGGCGTCACGGGCGGGAATGATTACCGGGGAATT‑TGTCATATTTATTTTTCTGTTTTTTATTTTACGGCTGCGCCTGTTc                                                                                                                                                            >  1:806394/1‑246 (MQ=255)
                                                                                      gctcacTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCGCCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATattaatt                                                                                                                                                                                             <  1:882109/205‑1 (MQ=255)
                                                                                      gCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATcc                                                                                                                                                                                         <  1:630611/212‑1 (MQ=255)
                                                                                           cTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGGGGGGGTAATGACGCGGTCGGCACCGGATTGTTTTGCTTGTTTAAAAAAGTGTTGAATTATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGGTTGCCCCGTAAACGGCT‑GGGGGGGGAAAGGGCGGGGATGATCA‑CGGGGAATTGTGTTATATTTATTATTCCGTTCTTTGTTTAAAGGCTGCGccccgtctctttttaaaatct                                                                                                                                              >  1:317945/1‑231 (MQ=255)
                                                                                                          ttCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTTGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTc                                                                                                                                                            <  1:444454/221‑1 (MQ=255)
                                                                                                          ttCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGc                                                                                                                                        <  1:174577/241‑1 (MQ=255)
                                                                                                                         ccaccaATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAAt                                                                                                                                                                                                            <  1:818281/158‑1 (MQ=255)
                                                                                                                          caccaATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGcc                                                                                                                                                                 <  1:129590/200‑1 (MQ=255)
                                                                                                                             caATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACggg                                                                                                                                             >  1:504468/1‑217 (MQ=255)
                                                                                                                                              cGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACg                                                                                          <  1:60621/251‑1 (MQ=255)
                                                                                                                                              cGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACg                                                                                          <  1:761323/251‑1 (MQ=255)
                                                                                                                                                aTAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATtg                                                                                                                                                                                                          <  1:102378/137‑1 (MQ=255)
                                                                                                                                                aTAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAAt                                                                                                                              >  1:585816/1‑213 (MQ=255)
                                                                                                                                                aTAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGtt                                                                                        >  1:638867/1‑251 (MQ=255)
                                                                                                                                                 tAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTgtgt                                                                                                                                                                                                       <  1:121593/139‑1 (MQ=255)
                                                                                                                                                      aCGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGAtgt                                                                                   <  1:451525/250‑1 (MQ=255)
                                                                                                                                                       cgcgGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGAtgttg                                                                                 <  1:404586/251‑1 (MQ=255)
                                                                                                                                                               gCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAAc                                                                                                                                                                          >  1:643974/1‑154 (MQ=255)
                                                                                                                                                                cACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCG‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCAGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGt                                                                         <  1:590503/250‑1 (MQ=255)
                                                                                                                                                                 aCCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTAt                                                                       >  1:262111/1‑251 (MQ=255)
                                                                                                                                                                     gATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTcc                                                                    <  1:511271/250‑1 (MQ=255)
                                                                                                                                                                      aTTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCaa                                                                  <  1:47215/251‑1 (MQ=255)
                                                                                                                                                                        tGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAAt                                                                                                                              <  1:747688/189‑1 (MQ=255)
                                                                                                                                                                        tGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCaaa                                                                 >  1:714062/1‑250 (MQ=255)
                                                                                                                                                                         gTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCg                                                                                                                                                                                                                                                   <  1:388688/73‑1 (MQ=255)
                                                                                                                                                                            tGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTcc                                                            <  1:701453/251‑1 (MQ=255)
                                                                                                                                                                                 tGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGAc                                                       <  1:714175/251‑1 (MQ=255)
                                                                                                                                                                                         aaCGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTggg                                               <  1:391171/251‑1 (MQ=255)
                                                                                                                                                                                           cGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCaaa                                                                 <  1:662397/231‑1 (MQ=255)
                                                                                                                                                                                           cGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCa                                             <  1:522480/251‑1 (MQ=255)
                                                                                                                                                                                            gtgtTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCAt                                            >  1:581575/1‑251 (MQ=255)
                                                                                                                                                                                                 gAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTc                                                                                                                                                            >  1:256867/1‑134 (MQ=255)
                                                                                                                                                                                                              gCCGGCAATATC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATcc                                                                                                                                                                                         >  1:254251/1‑92 (MQ=255)
                                                                                                                                                                                                                     atatC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCcat                   >  1:126938/1‑251 (MQ=255)
                                                                                                                                                                                                                      tatC‑GGCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCcatt                  <  1:693195/251‑1 (MQ=255)
                                                                                                                                                                                                                           ggCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGc                                                                                                                     >  1:820831/1‑148 (MQ=255)
                                                                                                                                                                                                                           ggCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGAtgt                                                                                   >  1:460262/1‑182 (MQ=255)
                                                                                                                                                                                                                           ggCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCcattcatt              >  1:683967/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCcattcatt              >  1:477702/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCcattcatt              >  1:869895/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCTAA‑GGGTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCCATTCATTa             >  1:188706/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gggTTTGCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTc                                                             <  1:651399/198‑1 (MQ=255)
                                                                                                                                                                                                                                        gCCCAGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCCATTCATTACGCTGCTCTCg  >  1:487792/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
ACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAATCATCGGCTTGCCGGCAATATC‑GGCTAA‑GGGTTTGCCCGGTAAGCGGCT‑GGAGGCGTAACGGGCGGGAATGATCA‑CGGTGAATTGTGTCATATTAATTATCCTGTTCTTGATTTAACGGCTGCGCCTGTTCCGCCGGTAATACGGGAATGCCGTTCTTAATCGGGTAAGCGATTTTTTCAAAACGACAGATTAAACGTTGATGTTGTTCATCGTATTCCAAACGTCCGTGACAGCTTGGGCATGCCATGATTTCTAATAATTTTCCATTCATTACGCTGCTCTCG  >  1/947977‑948457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑