Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,8550AC30.4% 347.5 / 24.9 149G247G (GGT→GGG) gpmA2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (79/24);  new base C (1/44);  total (80/69)
Fisher's exact test for biased strand distribution p-value = 5.61e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 153 total aligned reads displayed.

ACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786605‑787105
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
aCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTa                                                                                                                                                                                                                                                               <  1:1013266/251‑1 (MQ=255)
  cGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTAcc                                                                                                                                                                                                                                                             >  2:584788/1‑251 (MQ=255)
  cGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTAcc                                                                                                                                                                                                                                                             <  1:837445/251‑1 (MQ=255)
    aTGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCt                                                                                                                                                                                                                                                           >  1:809736/1‑251 (MQ=255)
      gCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTgg                                                                                                                                                                                                                                                         >  1:942331/1‑251 (MQ=255)
      gCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTg                                                                                                                                                                                                                                                          <  1:854604/250‑1 (MQ=255)
          gCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑ttt                                                                                                                                                                                                                                                     <  1:237076/250‑1 (MQ=255)
           cATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACCCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGc                                                                                                                                                                                                                                                   <  1:532024/251‑1 (MQ=255)
              cGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑c                                                                                                                                                                                                                                              <  1:957639/251‑1 (MQ=255)
              cGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑c                                                                                                                                                                                                                                              <  2:574499/251‑1 (MQ=255)
               ggCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑cc                                                                                                                                                                                                                                             <  2:89190/251‑1 (MQ=255)
               ggCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAACGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑cc                                                                                                                                                                                                                                             >  2:1045179/1‑251 (MQ=255)
                      aaGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCgctg                                                                                                                                                                                                                                      >  1:378364/1‑251 (MQ=255)
                              gttgctttgttacgcaagcgggCGATTGCAGTAAAGCCGCGTGGGGGTTAAGGTTTAATTGGCGGATGCGC‑GTGCAGGGTGGCGTACCCGATGAAT‑CACCACAGGCCGGAAAAGTCGCGCGGCGTCCAACCAGGCAAGTTGCTCTA‑TGGTTAGCAACAAAAAAGCC‑GCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGATTGCCTTTTCCTTCG‑CTTT‑CCCTGGGTTTGC‑AT‑CCGCCGCTGCTTTCGCt                                                                                                                                                                                                                              <  2:325123/231‑1 (MQ=255)
                              gATGCTCTTTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGGGGCGTATCCGATAAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCt                                                                                                                                                                                                                              <  1:156233/251‑1 (MQ=255)
                                tGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGc                                                                                                                                                                                                                            >  2:304476/1‑251 (MQ=255)
                                      tttACGACAACGGGCGAATGCAGTAAAAGCGCGTCGCGGGAAAGTTTTATTTGTCGGATGCGCCGG‑CAGAGTGGCTTACCCGATGAAC‑CCCCGCAGGCATGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑GTGTTAGCAACAAAAAAGCC‑GCCCCCCTTGCAGTCGGTTTTCTCATTTAAAACGAAGGTCGTTTTCTGCG‑CTTTACCCTGG‑TTTCC‑AA‑CCGCCGCTGCTTTCGCTGCGAt                                                                                                                                                                                                                         <  1:892574/247‑1 (MQ=255)
                                        gacgccaacgggagaaggCAGTAAAAGCGCGTCGCGGGTAAGTTTTAATTGTCGGAGGCGCCGT‑CAGGGGGCCGTATCCGATGAAC‑CACCACAGGCCTGATAAGCCGCGCAGCGCCGCATCAGGCAATGTGCCCCA‑TTGTTAGCAACAAAAAAGCC‑GCCTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCG‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                                                                                                                                     <  2:378364/233‑1 (MQ=255)
                                         ccGTCAACGGGCGAATGCAGTAAAGGCGCGCCGGGGGTAGGGTTTAATGGTCGGTTGCCCGGT‑CAGAGTGGCGTACCCGATGATC‑CACCACGGCCCTGATAAGGCGCGCGGCGTCCCATCAGGCAATGTGCTCCA‑TTTTTAGCAACAAAAAAGCC‑GACTCTCTTGCAGTCGGTTTTCTCTTTTAAAACGATTGACTTTTCCTGCC‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                                                                                                                                   <  1:513530/250‑1 (MQ=255)
                                                   gCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑CA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:957224/1‑251 (MQ=255)
                                                    cGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  1:140238/1‑250 (MQ=255)
                                                     gAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:1111226/1‑249 (MQ=255)
                                                     gAATGCACTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGTCGTATCCGCTGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:844827/1‑249 (MQ=255)
                                                       aTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCCGat                                                                                                                                                                                                     >  2:754783/1‑251 (MQ=255)
                                                        ggCAGAAAAGGCGCGTCGCGGTTAAGGTTTAATTGTCGGATGCGCCGT‑CAGTGTGGCGTATCCGATAAAT‑CACCACGGGCCTGATAAGTCGCGCAGCGTCGCATCAGGAAATGTCCTCCA‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAACCGAAGGAGGTTTACTTCG‑CTTTACCCTGT‑TTTGC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                                                                                                                    <  1:1111226/250‑1 (MQ=255)
                                                         gcattAAAAGCGCGTCGCGGGTAAGTTTAAATTGGCGGAGGCGCCGT‑CAGAGTGGCGTATCAGATGAAT‑CCCCACGGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTTCTCCA‑TTTTTAGCAACAAAAAATCC‑GCCTCCCTTGCATTCGGTTTTCTCATTTTAACCGAATGCCTTTTACTCCG‑CTTTCCCCTGT‑TTTGC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataa                                                                                                                                                                                                   <  1:647916/247‑1 (MQ=255)
                                                            gTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CCCCCCAGGCTTGATAAGTCGCGCAGCGTCGCACCAGGCAATGTGCTCCA‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                                                                                                                                <  1:396878/251‑1 (MQ=255)
                                                             aaaaaGCGCTTCGCGGTTAAGGTTAAATTGTCGGGTGCGCCGT‑CAGATTGGCGTATCCGAGGAAT‑CCCCACAGGCCTGAAAAGCCGCCCAGCGTCGCATCAGGAAATGTGCTCCA‑TTGTTAGCAAAAAAAAACCC‑GACTCCTTTGCATTCGGCTTTCTCTTTTAAAACGAATAACTTTTACTTCT‑CTTTCCCCGGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                                                                                                                               <  1:1025889/250‑1 (MQ=255)
                                                              aaaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                                                                                                                              >  1:573207/1‑251 (MQ=255)
                                                                   gggggggggggTAAGTTTTATTTGTGGGATGCGCCGT‑CAGAGTGGCGTACCCGATGAAT‑CCCCACAGGCTTGATAGGGCGCGCAGCGTCGAATCAGGCAATGTCCTCCT‑TTGTTAGAAAAAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCCCATTTTAAACGAATGCTGTTTCCTTCG‑CTTTCCCCTG‑‑TTTGG‑AATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                                                                                                                         <  2:724795/243‑1 (MQ=255)
                                                                   ggcggcggggttAAGGTTTAATGGGGGGATGCGCCGT‑CAGGGTGGCGTATCCGATGAAT‑CCCCACAGGCCTGATAAGTCGCGCAGCGCCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAAAAAAAAAGCC‑GCCTCACTTGCAGTCGGTTTTCTCTTTTTAAACAAATGCCTTTTACTCCG‑CTTTACCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                                                                                                                         <  1:1090115/240‑1 (MQ=255)
                                                                          ggggTAAGGTTTAATTGTCGGAGGCGCCGT‑CAGAGTGGCGTACCCGATGAAT‑CCCCCCAGGCCGGATAAGCCGCGCAGCGTCGCATCAGGCAATGTGCCCAA‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCGGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                                                                                                                                  <  1:304476/250‑1 (MQ=255)
                                                                           gggtaagggtaaATTGTGGGAGGGGCCGT‑CGGAGGGGCGAAGCCGATGAAC‑CCCCACAGGCCTGATAAGCCGCGCGGCGTCGAAGCAGGAAATGGGCTCCA‑TTTTTAGAAACAAAAAAGCC‑GTCTCACTGGCAGTCGGTTTTCTCTTTTTAAACAAATGACTTTTACTTCG‑TTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                                                                                                                                 <  2:53870/251‑1 (MQ=255)
                                                                            ggTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                                                                                                                                 >  2:283110/1‑250 (MQ=255)
                                                                             gtagggTTTAATTGTCGGAGGCGCCGG‑CAGTGTGGCTTATCCGAT‑AATACGCCCCAGGCCTGAAAAGTCGCGCAGGGTGGCATCAGGCAATGTGC‑CCAGTGGTGAGAAACAAAAAAGCC‑GACTCACTTGCTGTCGGTTTTCTCTTTTAAAACAAATGACGTTTACT‑CGTCTTTACCCTGG‑TTTGA‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCGTTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                                                                                                                               <  1:286009/247‑1 (MQ=255)
                                                                             gTAAGGTTTAATTGTCGGCTGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                                                                                                                                >  1:1008387/1‑250 (MQ=255)
                                                                                agttttAATTGTCGGTTGCGCCGT‑AAGAGTGGCGTACCCAATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGCAAATGTGCTCCA‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAtt                                                                                                                                                                            <  1:77722/248‑1 (MQ=255)
                                                                                   tttAATTGTCGGAGGCGCCGT‑CAGAGTGGCGTATCCGATGAAC‑CACCACAGGCCTGATAAGTCGCGCAGGGTCGCATCAGGAAATGTGCTCCT‑TTGTTAGCAAAAAAAAAGCC‑GCCTCACTTGCAGCCGGTTTTCTCATTTTAAACTAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                                                                                                                          <  1:317726/250‑1 (MQ=255)
                                                                                     tAATTGTCGGATGCCCCGT‑AAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTAATAAGTCGCGCGGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                                                                                                                        <  1:624160/250‑1 (MQ=255)
                                                                                          gTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                                                                                                                                  >  1:595730/1‑251 (MQ=255)
                                                                                              gATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                                                                                                                              >  2:770551/1‑251 (MQ=255)
                                                                                              gATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCAAGAGAATAGCGTTTCAGCGGTTTGAAATTATCGTCGAACTCAt                                                                                                                                                              >  2:598818/1‑251 (MQ=255)
                                                                                                                  gTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAAGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCCGCTCTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                                                                                                                           >  1:937618/1‑251 (MQ=255)
                                                                                                                  gTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTATCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTGCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGAGGCACGCCAGt                                                                                                                                           >  2:473965/1‑251 (MQ=255)
                                                                                                                        gATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGa                                                                                                                                      >  2:167561/1‑250 (MQ=255)
                                                                                                                         gagaaaccaccacAGGCTTGATAAGCCGCGCAGCTGCGCATCAGGCATTGTGCTCCA‑TGGTAAGCAAAAAAAAAGCC‑GCCTCACTGGCGGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑TTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                                                                                                                     <  1:214718/244‑1 (MQ=255)
                                                                                                                           gAAT‑CACCAAGGGCCTGATAAGTCGCGCAGCGTCGCACCAGGCAATGTGCTCAA‑TTTTTAGAAAAAAAAAAGCC‑GACTCACTTGCATTCGTCTTTCTCATTTTAAACGAATGCCGTTTACTCCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATAtt                                                                                                                                  <  1:98460/251‑1 (MQ=255)
                                                                                                                             aT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                                                                                                                                 >  2:310286/1‑250 (MQ=255)
                                                                                                                                 ccccaaggccggaaaagtccggcAGCGTCCCACCAGGCATTGGCCTCCA‑TTTTTAAAAAAAAAAAGGCC‑GACTCCCTTGCATTCGGCTTTCTCATTTAAAACGAATGTCGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                                                                                                                             <  2:192154/230‑1 (MQ=255)
                                                                                                                                   cacgggcctgtaaaggcgcgCCGCGTCGCATCAGGAAGTGCGCTCCA‑TTGTTA‑CAAAAAAAAAGCCTCACTCACTGGCAGTCGGCTTATTCATTTAAAACGAATGTGTTTTCCTTCG‑CTTTCCCTTGG‑TTTAC‑AT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                                                                                                                           <  2:777760/235‑1 (MQ=255)
                                                                                                                                        gCCTAATAAGGCGCGAGGCGTCCCATCAGGCAATG‑GCCCCAATTGTTGACAAAAAAAAACCC‑GCCCCGCTTGCGGGCGGTTTCCTCATTTAAACCAAATGCCGTTTACTCCG‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                                                                                                                      <  2:216871/251‑1 (MQ=255)
                                                                                                                                         ccTGATACGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCCGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTAGAACTCATACACCTGCGGCACGCCAGTCGTGCTATTACGCTCAAGAAtc                                                                                                                     >  2:789052/1‑251 (MQ=255)
                                                                                                                                         ccTGATAAGCCGCACAGCGACGCATCAGGCAATGTTCTCCT‑TTGTTAGCAACAAAAAACCC‑CTCCCCCTTGCAGTCGGCTTCCTCATTTTAAACGAATTACGTTTACGTCG‑TTTTCCCCGGG‑TGTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                                                                                                                     <  2:407529/251‑1 (MQ=255)
                                                                                                                                                 gccgcgCAGCGTGGCAGCGGGAAATTTGCTCAA‑TTGTAAGAAACAAAAACGCC‑GACTCCTTTGCGGCCGGCTTTCTAATTTTAAACAAAGGACTTTTCCTTCG‑TTTTCCCCTG‑‑TTTGCCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                                                                                             <  2:220096/249‑1 (MQ=255)
                                                                                                                                                 gTCGCGCGGCGTCGCATCAGGCAATGTGCCCAA‑TTGTTAGCAACAAAAAAGCC‑GACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                                                                                             <  1:1039399/251‑1 (MQ=255)
                                                                                                                                                     cgcAGCGTCCCATCAGGAAATGTCCTCCA‑TTGTAAGCAACAAAAAAGCC‑GCCTCTCTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCT‑CTTTCCCCTGG‑TTTGC‑AG‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                                                                                         <  2:1029401/251‑1 (MQ=255)
                                                                                                                                                      gcAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                                                                                         >  1:731144/1‑250 (MQ=255)
                                                                                                                                                      gcAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                                                                                         >  2:1153132/1‑250 (MQ=255)
                                                                                                                                                      gcAGCGCCGCAACAGGCATTGTGCTCCA‑TTTTTAGCAAAAAAAACGCC‑TCCTACCTTGCAGTCGGTTTTCCCTTTTTAAACGAATGCCTTTTTCTTCG‑CTTTCCCCTGG‑TTTAC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCa                                                                                                        <  2:1052613/251‑1 (MQ=255)
                                                                                                                                                        aGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg                                                                                                      >  1:232816/1‑251 (MQ=255)
                                                                                                                                                          cGTCCCATCAGGAATTGTGCTCCA‑TTGTTAGAAACAAAAAAGCC‑GACTCCCTTGCAGTCGGCTTTCTCATTTTAAACTAATGCCGTTTCCTCCG‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGtt                                                                                                    <  2:518790/251‑1 (MQ=255)
                                                                                                                                                                aTCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCa                                                                                                >  2:135749/1‑249 (MQ=255)
                                                                                                                                                                  cAGGCAATGTGCTCCT‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCG‑CTTTACCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                                                                                            <  1:413647/251‑1 (MQ=255)
                                                                                                                                                                  cAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCa                                                                                                >  1:376038/1‑247 (MQ=255)
                                                                                                                                                                        aTGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                                                                                       >  1:107217/1‑250 (MQ=255)
                                                                                                                                                                          gtgCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCAcc                                                                                    >  1:818338/1‑251 (MQ=255)
                                                                                                                                                                           taccccA‑TTGTTAGCAAAAAAAACCCC‑GACTCCCTTGCGGCGGGTTTTCTCATTTAAAAGGATGGGCTTTTCCTCCG‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCAcc                                                                                    <  2:1027958/246‑1 (MQ=255)
                                                                                                                                                                               ccA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcgc                                                                               >  1:802576/1‑251 (MQ=255)
                                                                                                                                                                                 a‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                                                                                                                               >  2:1086757/1‑153 (MQ=255)
                                                                                                                                                                                 a‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgc                                                                                 >  1:779951/1‑247 (MQ=255)
                                                                                                                                                                                 a‑TTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                                                                                                                               <  1:1086757/153‑1 (MQ=255)
                                                                                                                                                                                    tGTTCGCAACAAAAAAGCC‑GCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGt                                                                           >  1:558255/1‑251 (MQ=255)
                                                                                                                                                                                    tGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGt                                                                           >  2:36254/1‑251 (MQ=255)
                                                                                                                                                                                       tAGCAAAAAAAAACCC‑GACTCCCTTGCAGTGGGTTTTCGCATTTAAACCGAATGACTTTTCCTTGG‑TTTTACCCTGT‑TTTGC‑AT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaa                                                                         <  1:1020438/250‑1 (MQ=255)
                                                                                                                                                                                          aaaaaaaaaaGCG‑GACTCCCTTGCAGTCGTTTTTCTCATTTTAAACGAATGCCTTTTCCTCCG‑CTTTCCCCTGT‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAagag                                                                     <  1:754783/247‑1 (MQ=255)
                                                                                                                                                                                           aaaaaaaaaGCC‑GACTCGTTTGCAGTCGTTTTTCTCATTTTAAACGAATGATTTTTACTTCG‑TTTTCCCTTGT‑TTTGC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGt                                                                    <  2:908298/248‑1 (MQ=255)
                                                                                                                                                                                                  aaGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCgtg                                                             >  1:1062122/1‑251 (MQ=255)
                                                                                                                                                                                                    gCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGtgtg                                                           >  1:259634/1‑251 (MQ=255)
                                                                                                                                                                                                            catttgcggtgggCTTTCTCATTTTAAACGAAGGACGTTTACTTCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCg                                                      <  1:1045179/238‑1 (MQ=255)
                                                                                                                                                                                                               ttGCAGGCGGTTTTCCCATTTTAAAGGAATGTCTTTTACTCCG‑CTTTCCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatga                                                  <  2:53884/250‑1 (MQ=255)
                                                                                                                                                                                                               gtgcggtcgtttttCTCATTTAAAACGAATGCCTTTTCCTTCG‑TTTTCCCCTGT‑TTTCC‑AC‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatgat                                                 <  2:7361/240‑1 (MQ=255)
                                                                                                                                                                                                                 gCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATc                                                >  1:677139/1‑250 (MQ=255)
                                                                                                                                                                                                                   aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcg                                             >  1:884520/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCt                                         >  1:1132893/1‑251 (MQ=255)
                                                                                                                                                                                                                        gcttcctatatgaaaccGAATGAGTTTTACTCCT‑CTTTCCCCTGT‑GTTCC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTGCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTc                                        <  2:836753/235‑1 (MQ=255)
                                                                                                                                                                                                                           ttCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgcg                                                                                >  1:754805/1‑208 (MQ=255)
                                                                                                                                                                                                                            tctcATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCg                                    >  2:44198/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            >  1:428138/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaaCGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            >  2:617278/1‑250 (MQ=255)
                                                                                                                                                                                                                                          aaTGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCgg                      >  1:348250/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTGACGATTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATGTTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGCTGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  1:497060/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cGTTTTCTTCC‑CTTTACCTTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGaa                 <  2:408006/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGa                  >  1:12959/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAAt                >  1:1120859/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ttACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGCGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCTGCAGAAt                >  2:589141/1‑249 (MQ=255)
                                                                                                                                                                                                                                                       tcccgtttACCCTTG‑TTTGA‑CG‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGttt           <  2:159449/245‑1 (MQ=255)
                                                                                                                                                                                                                                                        cG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGTTCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCAt        >  1:389842/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             ttACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                                                                                                                           <  1:582927/132‑1 (MQ=255)
                                                                                                                                                                                                                                                               aCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGt  >  1:347877/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               aCCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACACGCTCACAGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGt  >  2:509136/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             
ACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGT‑CAGAGTGGCGTATCCGATGAAT‑CACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCA‑TTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG‑CTTTACCCTGG‑TTTGC‑AA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786605‑787105

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