Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,486 T→C 21.7% S135P (TCG→CCG)  torC → trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,4860TC21.7% 193.1 / 12.8 113S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (54/27);  new base C (5/22);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 1.50e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATT  >  NC_000913/1058252‑1058710
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGACGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATCCACC‑TGAAAAATTCCAAGCCAAAAC‑CCCGGAACTTACCGACGAGGAA‑TGGGGCGG‑AAGGAAAGAAAACAA‑CAC‑gca                                                                                                                                                                                                                               >  1:1135358/1‑237 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTTCCGAGCGTGAA‑TAGGCGCG‑AATGAAAGAAAACAA‑CTA‑gga                                                                                                                                                                                                                               >  2:656251/1‑238 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGc                                                                                                                                                                                                                               >  2:83012/1‑239 (MQ=255)
        aaaaTATAAAACGGCGGTGCATTTCCAGAACGCCGCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGc                                                                                                                                                                                                                        <  2:1092534/236‑1 (MQ=255)
          aaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGCCCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAAAG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑ACTGAAAGAAAAAAA‑CCC‑GGCAACCtccag                                                                                                                                                                                                                      >  1:999451/1‑234 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGc                                                                                                                                                                                                                     >  2:804429/1‑238 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CCC‑GGCAACCTGCCGc                                                                                                                                                                                                                     >  2:1281129/1‑238 (MQ=255)
                   ctgctgTGCATTTCCAGAGCGCCTCGGGGGTGAGAGCTAAATGCCAGAACGGTCATATCCGGCCGGATATGCCAGGCGTGGGGAAGGGAAAACGGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑GTGATCCACC‑TGAAAAATGCGAAGCCAAACG‑CGCGGAACTTGCCGGGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCAt                                                                                                                                                                                                           <  1:620182/234‑1 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCAATCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTTCCGATAGTGAA‑TGGGCGCG‑AATGAAAGAAAACCA‑CCC‑GGCAAACGCCGGC‑CCCTGCCaaa                                                                                                                                                                                                          >  1:718160/1‑237 (MQ=255)
                     tCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCCACCGTGAA‑TGGGCCCG‑AATGAAAAAAAACAA‑CCC‑gccccctcccttaacaaacacgg                                                                                                                                                                                                           >  1:524224/1‑216 (MQ=255)
                      cGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAAc                                                                                                                                                                                                        >  2:1189540/1‑239 (MQ=255)
                       ggTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGCTACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGCCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAAc                                                                                                                                                                                                        >  2:271471/1‑238 (MQ=255)
                          gCATTTCCAAAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCAAATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    <  2:654458/239‑1 (MQ=255)
                           cATTTCCAGACCCCCCCCGGCGTGGGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    <  2:513296/238‑1 (MQ=255)
                           attttCCAAACCGCCTCCGGCGGGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACCCC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    <  2:941854/236‑1 (MQ=255)
                               tccgaaccGCCTGCGGGGGGCGAGCTGAATGCCCGGACTTTACTATCCCGCCGGAATTTCCAGGCTGGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTGCCAGACCTTTAT‑TGCTCACTCCA‑TTGATACCCC‑TGAAAAATTCGAAGCCAAACG‑CGCGGACCTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcgatgcg                                                                                                                                                                                               <  2:987024/232‑1 (MQ=255)
                               tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcgatgcg                                                                                                                                                                                               >  2:945513/1‑239 (MQ=255)
                               tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcgatgcg                                                                                                                                                                                               >  2:42209/1‑239 (MQ=255)
                                         cTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGCTACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGGGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑ATC‑GGCAACTAGCCGC‑ACCTGCCAGAACTACGATGCAATGGAtaaag                                                                                                                                                                                     >  2:1171072/1‑235 (MQ=255)
                                              gCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑ACCGGAACTTGCCTACGGGGAA‑TGGGGGCC‑AATGAAAGAAAACAA‑CTC‑GGAAACCTCCCGC‑CCCTGCCATAACCACGACGCCATGGaaccagaacaa                                                                                                                                                                                >  1:1038012/1‑229 (MQ=255)
                                                  gggacatGAACGCCCTGACTGTCATATCCAGCCGGATATGCCAGGCATGGTGAAGCGCAAGCTGGAAGCGAGAAA‑TGATATCTACCAGCCCTTTAT‑TGCTCATTCCA‑TTGATCCACA‑TGAAAAATTCGAAGCCAAACG‑CGCGGACTTTGCCGAGCGTAAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                             <  2:960140/232‑1 (MQ=255)
                                                  gCGGGCTGAATGCAATGCTTTTCATTCCCCGCCGGATATTCCGGGATGGGTGAACCGCACCATTGAAGC‑TGCAATTGATACCTACAAGACCTTTCT‑TGCTC‑CTCCAGTTGAT‑CACCTTGAAAATTCCAAAGCAAAACG‑CGCGGAGCTTGCCGGGCGTGAA‑TGGGGCCA‑AATGAAAAAAAACAC‑CCC‑GGCACCCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                             <  1:656251/238‑1 (MQ=255)
                                                   cGAGCTGAATGCCATGACTGTCATATCCCTCCGGATATTCCAGGCATGGTCAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGCTACACC‑TGCAACATTCGAAGCCAAACG‑CGCGGAACTTGCCGTGCGCGTA‑TGGGCGCG‑AATGAAAGCAACCAA‑CTC‑GGCCACCTGCCCC‑CCATGCCCTCACCACCATGCACTAGTCCATGCgcagcacc                                                                                                                                                                            >  1:352573/1‑236 (MQ=255)
                                                          aaTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATccctg                                                                                                                                                                     >  2:1076633/1‑235 (MQ=255)
                                                                   aCTGTTATTCCCCGGCGGCTATTCCAGGACAGGGAAAGCCCAACCTGGAAGCGACAAA‑TGAAATTGCCCAGACCTTTAT‑TGCTCACCCCA‑TTAAACCCCC‑TGAAAATGTCGAAGCAAAACG‑CGCGGAACTTGCCGGCGGTGAA‑TGGGCGCA‑AATGAAAGAAAACAA‑CCC‑GGCACCCTGCCGC‑CCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGt                                                                                                                                                           <  1:833144/239‑1 (MQ=255)
                                                                    cTGTCATATCCCGCGGGATATTCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAA‑TGATATCTCACAGCCCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTACCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGt                                                                                                                                                           <  2:1135358/238‑1 (MQ=255)
                                                                               cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CCCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑ACTGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGCAGCAGCACGTCAGATGAAGGTg                                                                                                                                               >  2:1233533/1‑239 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    >  1:434548/1‑184 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    <  2:434548/184‑1 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAaacag                                                                                                                                                                             >  1:999582/1‑207 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCACCATCCTGAAGCAGCACGTCAGATGAAGGTGGc                                                                                                                                             >  2:41832/1‑239 (MQ=255)
                                                                                         cccAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTCCCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑GGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAg                                                                                                                                     <  2:154244/238‑1 (MQ=255)
                                                                                                   gTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCGGCATCCTGAAGCAGCACGCCAGAGTAAGGTGGCAGCGAAAGATaatcaat                                                                                                                            >  2:837641/1‑238 (MQ=255)
                                                                                                   gTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGCCAGATGAAGGGGGAAGCGAAAGATaaacaaa                                                                                                                            >  1:194929/1‑237 (MQ=255)
                                                                                                        gcgcAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATACCTACGATTCGATGGATCAGGCGAAGCACCATCCTGAAGCCACACGACAGCTCAAGGGGGCAGCGAAAACAAAACAATCCGGc                                                                                                                       >  1:206495/1‑238 (MQ=255)
                                                                                                             aatgggaagggcacAA‑TGATCTTCACCAAACCTTTAT‑TGCTCCCCCCA‑TTGATCCACC‑TGAAAAATGAGAAGCAAAACG‑CGCGGAACTTGCC‑AGTGGGGA‑AGGGCCAGAAAAGAAAGAAAACCA‑CCC‑GGCACCCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCGGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                                 <  1:83012/226‑1 (MQ=255)
                                                                                                             aaggggAAACGAGAAA‑AGATATATACGAGGCCTTGAT‑TGCCACCCCCA‑TTGATACGCC‑GAAAAAATTCGAAGCCAAACG‑CGCGGAACTGGCCGGGCGTGAA‑GGGGCGCG‑AAAGAAAGAAAACAA‑CGC‑GGCAACCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGAGGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                                 <  1:1189540/235‑1 (MQ=255)
                                                                                                             aaCTGGAAGCGAGCAA‑GTATATCCACCAGCCCTTTAT‑TGCTCACTCCA‑TTGATTCCAC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTTAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                                 <  2:999451/239‑1 (MQ=255)
                                                                                                                  tggaaccagagaatgtatCTTCCCGACCTTTATCGGCTCCTTCCA‑TGGATCCCCC‑GAAAAAATCCGAAGCAAAACGACGC‑GAAGTTGCTGAGGTGGAG‑TGAGCGCG‑AAGGAAAGAAACCAA‑CCC‑GGCAACGTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGt                                                                                                              <  1:971996/224‑1 (MQ=255)
                                                                                                                    aGCGAGCAA‑TGATATCACCCAGACCTTTAT‑GGCTCACTCCA‑TTGACACCCC‑TGAAAAATTCAAAGCCAAACG‑CGCGGAACGTGCCGGGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                                                                                          <  2:352573/239‑1 (MQ=255)
                                                                                                                            a‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGAGGCGCTGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGAAAATCAATCCGACATCGACTGTAATAAAGGTAt                                                                                                   >  1:43762/1‑238 (MQ=255)
                                                                                                                                      aCCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGCAACTTGCCGACCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCCTAACTACGATGCGATGGATCATGCGAAGCCGCACCCTGAGGCAGCACGCCAGAGGAAGGGGCCAGAAAAAAATAATCAAGCCGGCATCCACAGTCAAAAAAGTaggccccaacat                                                                                         >  1:24448/1‑228 (MQ=255)
                                                                                                                                        cAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGtt                                                                                       >  2:968269/1‑239 (MQ=255)
                                                                                                                                        cAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATCGATCATGCGAAGCAGCATCCTGAAGAAGCACGTCAGATGAAGGTGGCAGCGAACGAAAATCAATCCGCAATCAACTGTCATAAAAGAATTGCCACCCAgga                                                                                       >  1:1113179/1‑237 (MQ=255)
                                                                                                                                            ccTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGGATTGCCCACCAGTTACCg                                                                                   >  2:502477/1‑239 (MQ=255)
                                                                                                                                             cTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCgg                                                                                  >  2:592834/1‑239 (MQ=255)
                                                                                                                                               ttAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGGGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGat                                                                                >  2:481001/1‑239 (MQ=255)
                                                                                                                                                        cACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGCAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAg                                                                        >  2:1322125/1‑239 (MQ=255)
                                                                                                                                                        cACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGGGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACAAATTACCGGATATGAGAAg                                                                        >  1:1110767/1‑239 (MQ=255)
                                                                                                                                                           tCCA‑TTGATACACC‑GAAAAAAGTCGAACCCAAACG‑GGCGAAACGTGCCTAGCGAGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CCC‑GGAAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  2:1172626/239‑1 (MQ=255)
                                                                                                                                                           cccA‑TTGATACCCC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGGGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CCC‑GGCACCCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  2:1154541/238‑1 (MQ=255)
                                                                                                                                                           cccA‑TTGATACACC‑TGAAAAATTCAAAGCCAAACG‑CGCGGAACTTGCCGGGCGTGAA‑TGGGCGCA‑AATGAAAGAAACCAA‑CTC‑GGCACCCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  2:1223943/238‑1 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CTCGGAACTTGCCGAGCGTGAC‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATCCGCAGCAGCATCCTGAACCCGCACGGCAGCTGCAGGTGGCGGCGAACGATAATCAATCCTGCATCCGCTAGCATAAAGGCATTGCCCCCCAGTTAACGGctaactacagcagtttcccaaaa                                                           >  1:194651/1‑219 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGGGGAACCTGCCGGGCGGGAA‑TGGGCACG‑AATGAAAGAAAACAA‑CCC‑GGCAACCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                           <  2:482640/239‑1 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGc                                                                                                                                                                                                                               >  1:66519/1‑76 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGc                                                                                                                                                                                                                               <  2:66519/76‑1 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCCGTTACCGGATATGGGTAGCGGCTTCCGTAAg                                                           >  2:306132/1‑239 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                           >  2:835746/1‑239 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTCCATCGACTGTCATAAAGGCATTGCCCACCAGTTACCGGATATTAGTAGCGGCTTCCGTAAg                                                           >  2:338768/1‑239 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCGGCATCGACTGTCCTAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTaa                                                            >  2:209040/1‑238 (MQ=255)
                                                                                                                                                                      cacC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAAATGAAGGGGGCAGCGAAAGATAATCAATCCTGCATCGCCTGCCATAAAGGAATTGCCCACCAGTTACCGGAAATGAGAAGCGGCTTCCGAAAg                                                           >  1:1116831/1‑239 (MQ=255)
                                                                                                                                                                        cc‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATAAAGGATTTGCCCACCAGTTACCGGATATGAGGAGCGGCTccccaaag                                                           >  1:610951/1‑232 (MQ=255)
                                                                                                                                                                        cc‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAg                                                                                                                                                  >  2:199448/1‑150 (MQ=255)
                                                                                                                                                                        cc‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGAAGCGGCTTCCGTAAGCa                                                         >  2:70235/1‑239 (MQ=255)
                                                                                                                                                                        cc‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCa                                                         >  2:1107425/1‑239 (MQ=255)
                                                                                                                                                                        cc‑TGAAAAATTCGAAGCCAAACG‑CGAGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATCACTACGATGCGATGGATCATGCGAAGAAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGGACGGTACTCAATCCTACACCCACTGTTATAAAGGAATTGCCCCTCAGATTCCGGAAAGGGGTAGCGGCATCCGaaacga                                                         >  1:1079406/1‑236 (MQ=255)
                                                                                                                                                                         c‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAg                                                        >  2:1269335/1‑239 (MQ=255)
                                                                                                                                                                            gAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAcc                                                                                    >  2:1061043/1‑209 (MQ=255)
                                                                                                                                                                             aaaaaTTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGCTAATCAATCCTGCATCGACCTTCCTAAAGGAATTCCCCACAAGTTACCGGATAGGAGTAACGGATTcagaaaaaaagac                                                     >  1:29759/1‑227 (MQ=255)
                                                                                                                                                                             aaaaaTTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCccccca                                                                                          >  2:647779/1‑202 (MQ=255)
                                                                                                                                                                              aaaaTTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCg                                                                      >  2:911019/1‑221 (MQ=255)
                                                                                                                                                                                aaTCCGAAGCCAAACG‑CGGGGAACTTGCCGAGCGTGAA‑TGGGGGCG‑AATGAAAGAAAACAA‑CTC‑GGCACCCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGAt                                                  <  2:63599/239‑1 (MQ=255)
                                                                                                                                                                                       aacccaagcccgggactcgtcccGGGCGCGAA‑GGGGCGCA‑AATGAAAAAAAACAA‑CTC‑GGCACCTGGCCGT‑CCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAAGCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTg                                            <  1:70235/218‑1 (MQ=255)
                                                                                                                                                                                             aCG‑CGCGTACCTTGCCG‑GCGCGAATTGGGCGCG‑ATAAAAAGAAAACAC‑CGC‑GGCAACCTCCCTC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCa                                     <  2:1027597/239‑1 (MQ=255)
                                                                                                                                                                                              ccgcccctcacctgcgggggttGAG‑TGGGCGGA‑AAAAAAAGAAAACAC‑CTC‑TGCCACCTGCCGC‑GCCTGGCATAACTCCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGGAGATCAGATGAAGGTGGCAGCTAAAGATAATCAATCCTGGATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCCTCCGTAAGCAGTTCGATGAGCTGCTCGCCa                                     <  1:1248309/218‑1 (MQ=255)
                                                                                                                                                                                                    ggAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAACGGGATTGCCCACACGTAACCGGAAATGAGTAGCGGCTTCCCTAACCCGTCCGAtagatcgggggccagggaca                               >  1:813004/1‑220 (MQ=255)
                                                                                                                                                                                                        cTTGCGGGGGGTGAA‑TGGGCGCG‑AATGAAAAAAACCAC‑CCG‑GGAAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                           <  2:194929/239‑1 (MQ=255)
                                                                                                                                                                                                        cTTGCCGGCGGTGAA‑TGGGCGCA‑AATGAAAGAAAACAA‑CTC‑GGAAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                           <  2:1110767/239‑1 (MQ=255)
                                                                                                                                                                                                        cTTCCCGGGCGTGAA‑‑GGGGGCGAAATGAAGAAAAACAA‑GCC‑GGCAACCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGGACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                           <  2:194651/239‑1 (MQ=255)
                                                                                                                                                                                                          tttgccggggtggagggggccaaaaTGAAA‑AAAACAACCCC‑GGCACCCCGCCCC‑CCCGGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                           <  1:209134/216‑1 (MQ=255)
                                                                                                                                                                                                           gccgtggggaaa‑TGGGGGCG‑AATAAAAAAAAACAA‑CCC‑GGAAACCCGCCGC‑TCCTGCCATAACTTCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAg                        <  2:1260555/229‑1 (MQ=255)
                                                                                                                                                                                                           gccggggttGAA‑GGGGCGCA‑AAGGAAAGAAAACAA‑TCG‑GGCACCCTGCCCC‑CACTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAg                        <  2:854564/231‑1 (MQ=255)
                                                                                                                                                                                                            ccGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgg                       >  2:209134/1‑238 (MQ=255)
                                                                                                                                                                                                                gCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAGGGCATTGCCAACCAGTTACCGGATATGAGGAGCGGTTACCGAAACAAATTCAAAGAGCTacgccaccgagcaaaagacgaggggt                   >  1:41968/1‑217 (MQ=255)
                                                                                                                                                                                                                  gTGAA‑TGGTCGCG‑AATTAAAGAAAACAA‑CTC‑GCAATCCTGCCGC‑TCCTGCCATAACTTCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGAGAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac                 <  2:410061/239‑1 (MQ=255)
                                                                                                                                                                                                                      aaggggCGCG‑AATGAAAGAAACCAA‑CCC‑GGAAACCTGCCGC‑CCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg              <  2:157996/236‑1 (MQ=255)
                                                                                                                                                                                                                       aggggCGAA‑AATGAAAGAAAACAA‑CCC‑GGCACCCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg              <  2:1038012/236‑1 (MQ=255)
                                                                                                                                                                                                                       aatggggcgaaaggaaagaaaCAAC‑CTC‑GGGCACCTGCCCC‑CCCTGCCATAACTACGATGGGATGGGTCATGCGAAGCAGCATCCTGAAGCAGCAAGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg              <  1:306132/226‑1 (MQ=255)
                                                                                                                                                                                                                         tgggggaaaaataaagaaaAAAC‑CCC‑GGCACCATGCCTCGCCCGGCCACAACTACGAGCCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCGGCATCGGCTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg              <  1:945513/226‑1 (MQ=255)
                                                                                                                                                                                                                         tggCACG‑AAAGAAAAAAAAAAA‑CCC‑GGCACCGAGCCGC‑CCCTGCCAAAACTACGATTCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAAGCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt            <  1:338768/237‑1 (MQ=255)
                                                                                                                                                                                                                          gcgcag‑AATAAAAAAAAACAA‑CTC‑GGCAACCTGCCTG‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt            <  2:734854/232‑1 (MQ=255)
                                                                                                                                                                                                                            cgcg‑AATGAAAGAAACCAA‑CCC‑GGCAACCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:331763/238‑1 (MQ=255)
                                                                                                                                                                                                                             gcg‑AAAGGAAAAAAAAAA‑CTCGGGAAACCGGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:1079406/238‑1 (MQ=255)
                                                                                                                                                                                                                             gcaaaagtaagaaaaaCAC‑CCC‑GGCAACCGGCATC‑CCCTGCCAAAACGACGATGCGATGGCTCAGTCGAAGCAGTATCCTGAATCAGCACGTCAGATGAAGGGGGCAGCGAAAGATAATCAATCCGGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  1:271471/226‑1 (MQ=255)
                                                                                                                                                                                                                             cgcaaagaaaaaaaaaCAC‑CCC‑GCAAACCGGCCCC‑CCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:610951/230‑1 (MQ=255)
                                                                                                                                                                                                                             actaaaaaaaaaaaaaCAA‑CCC‑GGCATACAGGCGT‑CCCTGCCATAACTGCGATGGGATGGATCAGGCGAAGAAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAAGCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGGGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAc        <  1:929508/231‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTGAAAGAAAACAA‑CTC‑GGCACCCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    <  2:919664/47‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTAcg                                                                                                                                                                                                    >  1:919664/1‑47 (MQ=255)
                                                                                                                                                                                                                                  aTGAAAGAAAACAA‑CCC‑GGACACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTAt   <  2:29759/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tGAAAGAAACCAA‑CCC‑GGCAACCGCCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTAtt  <  2:1113179/239‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAA‑TGATATCTACCAGACCTTTAT‑TGCTCACTCCA‑TTGATACACC‑TGAAAAATTCGAAGCCAAACG‑CGCGGAACTTGCCGAGCGTGAA‑TGGGCGCG‑AATGAAAGAAAACAA‑CTC‑GGCAACCTGCCGC‑TCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATT  >  NC_000913/1058252‑1058710

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