Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,8550AC39.1% 109.4 / 16.5 70G247G (GGT→GGG) gpmA2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (38/4);  new base C (1/26);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 9.59e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTC  >  NC_000913/786611‑787101
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
gCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCGTTACCCtt                                                                                                                                                                                                                                                   <  2:512675/250‑2 (MQ=255)
   tGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGttt                                                                                                                                                                                                                                               >  2:208662/1‑251 (MQ=255)
     cATCGGCTGGGCAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTACGGGTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGCTGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCCTCAAGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CTTTTTAAACGCATGACGTTTACTTAGCTTTCCCCTGGTTTGc                                                                                                                                                                                                                                             >  2:71503/1‑251 (MQ=255)
      aTCGGCTGGGGAGATACCGATGCTCTGTTACGTCAACTTGCGAATGCAGTAAAAGC‑GCGTCGCGGGGAAGGTTTAAGTGTCG‑GAAGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑TGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc                                                                                                                                                                                                                                             >  2:291273/1‑250 (MQ=255)
                              ttGTTACGTCAACGGGGGAATGCAGAAAAAGC‑GCGTGGGGGGTAAGGTTTAATTGTCG‑GATGGGC‑CGTCGGGGGGGCGTTCCCGAGGAATCCCCACAGGCCTGATAGGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAAAAAAAAAGCCGACT‑CACTTGCAGCCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                                                                                                                                     <  2:552033/249‑1 (MQ=255)
                                gTTACGTCAACTGGCAAATGCAGTAAAAGC‑GGGTGGCGGGTAAGGTTTAATTGTCG‑GAGGCCC‑CGTCAGGGTGGCGTACCCGATGAATCACCCCAGGCCGGATAAGTCGCGCGGCGTCGCATCAGGCAAG‑GTGC‑TCCATTGTT‑AGAAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CC‑CATTTAAAACGAATGACGTTTACTTCGCTTTCCCCGGGTTTGCAGCCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                                                                                  <  2:175887/251‑1 (MQ=255)
                                    cgtaagcgggCGAATGCAATAAAGGC‑GCGTGGCGGTTAGGGTTTAAT‑GGCGTGCGGCGC‑CTTCAGAGTGGCGTCTCCGATGAATCACCACAGGCCTGAAAGGTCGCGCAGCGTCGCATCAGGCAGT‑GTGC‑TCCATTTTT‑AGCAACAAAAAAGCCGCTT‑CACTTGCAGTGGTCTTT‑CT‑CATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                                                                                                                              <  2:249420/243‑1 (MQ=255)
                                      tCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCAACAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                                                                                                                            >  1:80609/1‑251 (MQ=255)
                                         acgggcaaaTGCAGTAAAGGC‑GCGTGGGGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGGGGCGTATCCGATGAATCCCCACAGGCCTGATAGGGCGCGCGGCGTCGCATCGGGCAAT‑GTGC‑CCCATTTTT‑AGAAAAAAAAAAGCCGCCT‑CCCTTGCAGTCGGTTTT‑CT‑CATTTAAAACGAATGCCTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                                                                                                                                         <  2:38570/244‑1 (MQ=255)
                                               gAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                     >  1:266570/1‑249 (MQ=255)
                                                aaTGCAGAAAAAGC‑GCGGCGGGGGTAAGTTTTAATTGGCG‑GATGCGC‑GGTCAGAGGGGCGTACCCGATAAATCCCCACAGGCATGAAAAGCCGCGCGGCGTCGAATCGGGCAAT‑GTGC‑CCCATTGTA‑AGCAAAAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACTTTTCCTTCGCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGa                                                                                                                                                                                                  <  2:393755/251‑1 (MQ=255)
                                                       aaaaaGC‑GCGTGGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAAAAAAAAAGCCGACT‑CACTTGCAGTGGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                                                                                                                           <  1:551493/250‑1 (MQ=255)
                                                        aaaaGCCGCGTCGCGGTAAAGGTTTAATGGTGG‑GATGCCC‑CGTCAAGGTGGCTAATCCGATAAATCCCCACAGGCCGGATAAGTCGCGCAGCGACTCATCAGCTATT‑GTGC‑CCCATTGTG‑AGCAAAAAAAAACCCGCCT‑CCCTGGGAGCCGGCTTT‑CT‑CTTTTAAAACAAATGACGTTTCCTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                                                                                                                           <  2:522834/251‑1 (MQ=255)
                                                            ggcgggTCGCGGTAAAGGTTTAATTGTCG‑GA‑GGGCGCGTGAGGGGGGCTTATCCGATGAATCCCCCCAGGCCTGATAGGCCGCGCGGCGTCGCAGCGGGAAAT‑GTTC‑TCCTTTGTT‑AGCAAAAAAAAGGCCGCCT‑CACTTGCAGTCGGCTTT‑CT‑CATTTAAAACGAATGACGTTCCCTTCGCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                                                                                                                       <  2:564339/246‑1 (MQ=255)
                                                            gc‑gcGTCGGGGGTAAGGTTTAATTGTCG‑GATGCCC‑CGTCAGAGTGGCGTATCCAATGAATCCCCCCAGGCGTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGTTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                                                                                                                      <  1:544746/251‑1 (MQ=255)
                                                                cgTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                                                                                                                   >  1:177965/1‑251 (MQ=255)
                                                                 gTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGGTTATCCGAGGAATCCCCACAGGCCTGATAAGTCCCGCAGCGTCGCATCAGGCATG‑GTGC‑TCCATTGTT‑AGCAAAAAAAAAGCCGCCT‑CACTTGCAGTCGGCTTT‑CT‑TATTTAAAACGAATGACGTTTACTTCTCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                                                                                                                                  <  1:298660/251‑1 (MQ=255)
                                                                  ccgcgGTTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGATGGGCGTATCCGATGAATCCCCACAGGCCTAATAAGTCGCGCGGCGCCGCATCAGGCAAT‑GGGC‑TC‑AATGTTGAGCAACAAAAAAGCCGACT‑CCCTTGCAGT‑GGGTTT‑CTTCTTTTGAACCAAATTACTTTTACTTCGCTTTGCCCGGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                                                                                                                                 <  1:280901/250‑1 (MQ=255)
                                                                   ggcgGGTAAGTTTTAATTGCCG‑GATGCGC‑CGTCAGAGGGGCTTTCCCGATAAACCCCCACAGGCCTGATAAGTCGCGCGGCGCCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGAAAAAAAAAAGCCTCCC‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGCCTTTTACTCCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                                                                                                                                <  2:392019/250‑1 (MQ=255)
                                                                             gTTTAATTGCCG‑GATGCGC‑CGTCAGAGTGGCGTATACGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCCTCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                                                                                                                      >  2:295403/1‑251 (MQ=255)
                                                                              tttAATTGTCG‑GATGCGC‑CGTCCGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCCTCAGGCACT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGCATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtctc                                                                                                                                                                     >  2:518056/1‑251 (MQ=255)
                                                                                tAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                                                                                                                   >  1:84755/1‑251 (MQ=255)
                                                                                aatttGTCG‑GAGGCGC‑CGCCAGGGTGGCGTATCCGATAAATCCCCACAGGCCTGATAGGCCGCGCAGCGTGGCACCAGGCATT‑GTCC‑TCCATTGTT‑AGCAACAAAAAAGCCGTCC‑CGTTTGGAGTGGGCTTT‑CT‑CTTTTTAACCGATTGAGGTTTCGTTCGCTTTACCCGTTTTTGCAACCGCCGCTGCTTTCGCTGCGATCGCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                                                                                                                    <  1:577236/247‑1 (MQ=255)
                                                                                              gcgccgtcagggggggtaacccgttGAATCCCCACAGCCCTGATAGGTCGCGCAGCTTCGCATCAGGCA‑T‑TTGCCTCCATGGTT‑AGCAAAAAAAAACCCCCCC‑CACTGGCGGTCGGTTTT‑CT‑CATTTAAAACGAATAGCTTTTCCTCCGCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                                                                                                                                       <  1:411376/227‑1 (MQ=255)
                                                                                               gc‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                     >  2:207312/1‑251 (MQ=255)
                                                                                               gc‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                     >  2:434850/1‑251 (MQ=255)
                                                                                                    tCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                                                                                                                 >  1:61881/1‑251 (MQ=255)
                                                                                                     cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                                                                                                                >  2:545687/1‑251 (MQ=255)
                                                                                                     cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                                                                                                                >  1:215214/1‑251 (MQ=255)
                                                                                                       gagTGGCGTACCCGATGAATCCCCACAGGCCTGATAAGCCGCGCGGCGTCGCATCAGGCATT‑GTGC‑TCCATTGTA‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCa                                                                                                                                              <  1:278430/251‑1 (MQ=255)
                                                                                                        agTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                             >  1:432896/1‑251 (MQ=255)
                                                                                                             cGTATCCGATGAATCACCACAGGCTTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCAATTGTT‑AGCAACAAAAATGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                                                                                                                        <  1:158382/251‑1 (MQ=255)
                                                                                                              gTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCCTCAGGCAAT‑GTGC‑TCCCTTGTT‑AGCAACAAAAAAGCCGCCT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                                                                                                                       >  1:265239/1‑251 (MQ=255)
                                                                                                                    gATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                                                                                                                 >  1:210008/1‑251 (MQ=255)
                                                                                                                    gATGAACCCCCACGGGCTTGATAGGTCGCGCAGCGTCGCATCAGGAAAT‑GTGC‑CCCATTTTA‑AGAAAAAAAAAAGCCGACT‑CCCTTGCAGTGGGCTTT‑CT‑CATTTAAAACGAAGGACGTTTCCTTCGCTTTCCCCGTGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACGCCAGCGGCACGCCAGTCGGGat                                                                                                                                 <  1:235324/251‑1 (MQ=255)
                                                                                                                           ccccccgggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGAAAAAAAAAAGCCGCCT‑CACTTGCAGTCGGCTTT‑CT‑CATTTAAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                                                                                                                          <  1:212731/244‑1 (MQ=255)
                                                                                                                            ccaccagccttGATAAGTCGCGCAGCGTCGCAT‑AAGCAATTGTCC‑TCATTTGTT‑AGAAAAAAAAACGCCGACT‑CACTGGCA‑TCCTCTTTTCT‑CATTTTAACCGATTGCTGTTTCCTCCCTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                                                                                                                         <  2:132179/241‑1 (MQ=255)
                                                                                                                                caGGCCTGAAAGGTCGCGCGGCGCCACATCAGGCAAT‑GTGC‑TCCATTGTA‑AGAAAAAAAAAAGCCGCCT‑CACTTGCGGTCGGCTTT‑CT‑CATTTTAACCGAATGACGTTTCCTTCGCTTTCCCCTGTTTTCCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAg                                                                                                                     <  2:177965/251‑1 (MQ=255)
                                                                                                                                   gCCTGATAAGTCGCGCAGCGCCGCATCAGGCAAT‑GTCC‑TCCATTGTT‑AACAACAAAAAACCCGCCT‑CCCTTGCAGTCGGCTTT‑CT‑CATTTAAAACGAATGCCTTTTCCTTCGCTTTCCCCTGGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                                                                                                                  <  2:37078/251‑1 (MQ=255)
                                                                                                                                    ccTGATAAGTGGCGCAGCGTCGCATCAGGCAAT‑GTCC‑TCCATTGTT‑AGCAAAAAAAAACCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTAAAACGAATGACGTTTACTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                                                                                                                 <  2:186611/251‑1 (MQ=255)
                                                                                                                                       gATAAGTCGCGCAGCGTCGCATCAAGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGACGGGATATTAAGCTCAAGAAtccct                                                                                                              >  2:554082/1‑248 (MQ=255)
                                                                                                                                            gTCGCGCAGCGTCGCATCAGGCATT‑GTGC‑TCCTTTGTT‑AGCAACAAAAAAGCCGACT‑CCCTTGCAGTCGTCTTT‑CT‑CATTTTAAACGAATGTCTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                                                                                         <  1:290960/251‑1 (MQ=255)
                                                                                                                                                      gTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGtt                                                                                                >  2:84099/1‑250 (MQ=255)
                                                                                                                                                            tCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCAAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCa                                                                                            >  1:194169/1‑248 (MQ=255)
                                                                                                                                                                   aT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                                                                                   >  2:436452/1‑250 (MQ=255)
                                                                                                                                                                     gttgC‑TCCTTTGTA‑AGCAAAAAAAAACCCGACT‑CGTTTGCAGTCGGCTTT‑CT‑CATTTAAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGATACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCAc                                                                                 <  1:484970/249‑1 (MQ=255)
                                                                                                                                                                         c‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAgc                                                                             >  2:554612/1‑251 (MQ=255)
                                                                                                                                                                                      gCAACAAAAAAGCCGACTGCAAGTG‑AGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                                                                                                                >  2:140648/1‑173 (MQ=255)
                                                                                                                                                                                       caacaaAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAag                                                                   >  1:501507/1‑249 (MQ=255)
                                                                                                                                                                                        aacaaAAAAGCCGACT‑CACTTGCAGTCGGATTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGt                                                                >  2:482445/1‑251 (MQ=255)
                                                                                                                                                                                                   cGACT‑CACTTGCAGTCGGTTTT‑CT‑CATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCa                                                     <  1:107003/251‑1 (MQ=255)
                                                                                                                                                                                                   cGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGc                                                      >  2:183053/1‑250 (MQ=255)
                                                                                                                                                                                                      cT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGCCGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCg                                                  >  2:270143/1‑251 (MQ=255)
                                                                                                                                                                                                         cccttccgtccggtttt‑CT‑CATTTAAAAGGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgat                                                <  1:295403/237‑1 (MQ=255)
                                                                                                                                                                                                         cACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgat                                                >  2:258897/1‑251 (MQ=255)
                                                                                                                                                                                                           cTTGCAGAAGGCTTT‑CT‑AATTTTAAACGCCTGACGCTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCAAGATAATCGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatg                                               >  2:339684/1‑250 (MQ=255)
                                                                                                                                                                                                              gCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATc                                            >  2:9861/1‑250 (MQ=255)
                                                                                                                                                                                                                aGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcg                                         >  1:378853/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggcggtttgct‑cATTTAAACCGAACAACGTTTCCTTCGCTTTCCCTTGGTTTCAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCGCATGTTATCAAGATATGTCACCAGCGCGCGTAAAGAGTTACCGTGTGCAGCGATGATCAcgc                                        <  1:518056/242‑1 (MQ=255)
                                                                                                                                                                                                                        tt‑CT‑CATTTTCAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACAAGCGGCACGCCAGTCGGGATATTACGCTCAAGCATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCAcc                                 >  2:353991/1‑251 (MQ=255)
                                                                                                                                                                                                                           ct‑cATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGc                               >  1:396541/1‑251 (MQ=255)
                                                                                                                                                                                                                              cATTTTAAACGAAGGACTTTTACTCCTCTTTCCCCGGGTTGGCAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGct                              <  1:246303/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ttAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTc                          >  1:251105/1‑250 (MQ=255)
                                                                                                                                                                                                                                  taaaaggattgcgttttccTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTc                          <  1:297690/232‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg                     >  1:553769/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg                     >  2:309224/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aCGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                 >  1:511588/1‑248 (MQ=255)
                                                                                                                                                                                                                                                tttACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAAt            >  1:124869/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      tCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGttt       <  1:103274/249‑1 (MQ=255)
                                                                                                                                                                                                                                                         cTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTc  >  1:259164/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGC‑GCGTCGCGGGTAAGGTTTAATTGTCG‑GATGCGC‑CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAT‑GTGC‑TCCATTGTT‑AGCAACAAAAAAGCCGACT‑CACTTGCAGTCGGCTTT‑CT‑CATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTC  >  NC_000913/786611‑787101

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