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

GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786611‑787105
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
gCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTgg                                                                                                                                                                                                                                                        <  1:527675/251‑1 (MQ=255)
                                gttacgaaaaCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTAAGAGTGG‑CGTATCCGATAAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGAAACAAAAAAGCC‑GACTCTCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                                                                                                                       <  1:328186/242‑1 (MQ=255)
                                gTTACTTCAACTGGAGAATGCAGTAAAAGCGCGTCGCGGGTAGGGTTTAATGGTCGGATGCGCCGCCAGAGGGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTGGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTACCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                                                                                      <  1:225234/251‑1 (MQ=255)
                                   aCGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                                                                                                                                   >  1:301420/1‑251 (MQ=255)
                                         gcTGGCGAATGCAGAAACAGCGCGGCGCGGGTAAGGTGTAATAGGCGGATGCGCCGAGAGGGGGG‑CGTATCCGATGAATCACCACTGGCCTGATAAGT‑CGCGCAGCGCCGCTTCAGGCAATGTGCCCCATTGTTAGAAAAAAAAAAGCC‑GACCCACTTGCAGTCCGCTTTCTCTTTTAAAACGAATTTCGTTTACTTCGCTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                                                                                                                                             <  2:143712/250‑1 (MQ=255)
                                           gggggaaagCAGAAAAAGCGCGCCGGGGGTAAGGTTTAATTGCGGGATGCGCCTACAGGGTGG‑CGTATCCGATGATACACCACAGGCCTGAAAAGT‑CGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTGGCAGTCGTCTTTCTCATTTAAAACGAATGACGTTTTCTTCGTTTTACCCTGGT‑TTAC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                                                                                                                           <  1:354249/243‑1 (MQ=255)
                                            ggggAATGCAGAAAATGGGCGTCGCGGGAAAGTTTTAATGGTGGGATGCGCCGGCAGAGTGG‑CGTACCCGAAGAATCACCACAGGCCGGAAAGGC‑CGCGCAGCGTCGCATCGGGAAGGGTGCTCCATGGTTAGAAAAAAAAAAGCC‑GCCTCCCTTGCAGTCGTCTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                                                                                                                          <  2:287654/248‑1 (MQ=255)
                                              agAATGCAGTAAAAGCGCGTCGCGGGTAAGGGTTAATTGTCGGATGCGCTGTCAGAGTGG‑CGTATCCGATGAATAACCACAGGCCTGATAAGT‑CGCGCAACGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGAC‑GACACACTTGAAGTCGGCTTTCTCATTTTAAACGGATGAAGTTTACTTCGCTTTACCCTGGC‑TTGC‑AACCGCGGATGCTTTCACTGCTCTCTGGTCCGCATTAccc                                                                                                                                                                                                         >  2:302134/2‑250 (MQ=255)
                                               gAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGACGGATGCGCCGTGAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TGGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCATCATTAccc                                                                                                                                                                                                         >  2:112834/1‑249 (MQ=255)
                                                aaTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:542457/1‑248 (MQ=255)
                                                    aagaaaaagggggTCGCGGTTAGGGTTAAATTGTCGGATGCCCCGTCAGAGTGC‑CTTACCCGTTAAACCCCCACAGGCCTGATAAGT‑GGCGCAGCTCCGCATCAGGCAATGTGCTCCATTGTTATCAACAAAAAAGCA‑GCCTCCCTTTCATTCGGCTTTCTCATTTTAAACGAATGACTTTTTTTTCGCTTTTCCCGGGT‑TTTC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                                                                                                                                  <  1:547050/239‑1 (MQ=255)
                                                      gTAAAAGCGCGTCGCGGGTAAGGTTTACTTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGGTAGCAACAAACAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCGGCTTTCGCTTCGATCTCGTCAGCATTAACCAGataat                                                                                                                                                                                                  >  2:327877/1‑249 (MQ=255)
                                                         aaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                                                                                                                             >  2:4126/1‑251 (MQ=255)
                                                          caGCGCGTCGCGGGTAAGGATTTATTGTCGGCTGCGCCGTCCGATTGG‑CGTATCCGCTGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACACAAAAACC‑TACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATTACGATTACTTCTCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATATCGTCAACATTACCCCGATAATCGCGtt                                                                                                                                                                                            >  2:382097/2‑251 (MQ=255)
                                                             cgcgTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTcc                                                                                                                                                                                         >  1:182839/1‑250 (MQ=255)
                                                                 gcgcgGGAAAGTTTTAATTGTGGGAGGCGCCGGCAGAGTGG‑CGTATCCGAAGAACCCCCACGGGCCTGAAAAGT‑CGCGCAGCGTCGCATCGGGAAATGTGCCCCATTGTTAGCAAAAAAAAAGCC‑GACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTCCGTTTTCCCCTGTT‑TGGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                                                                                                                                     <  2:484049/250‑1 (MQ=255)
                                                                  ggctggTAAGGTTTATTTGCCGGATGGGCCGTCAGAGTGG‑CGTAGCCGATGAATCACCACAGGCTTGATAAGC‑CGCGAAGGGTGGCACCGGGAAAGTTGCTCATTTGTAAGAAAAAAAAAAGCC‑GACTCATTTGCATGCGGTTTTCTCTTTTTAAACGAATGGCTTTTCCTTCGCTTTACCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGGTTCAGCGGt                                                                                                                                                                                    <  1:409036/247‑1 (MQ=255)
                                                                  cgcgGTTAAGGTTTAATGGGCGGAGGCGCCGTCAGGGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGC‑CGCGCGGCGTCGCACCAGGAAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                                                                                                                                    <  1:100530/251‑1 (MQ=255)
                                                                    ggggTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCCCCACAGGCCTGATAAGT‑CGCGCAGCGTCGCACCAGGCAATGGCCTCCATGGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGGCGGCTTTCTCATTTAAAACAAATGCCGTTTACTTCGCTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                                                                                                                                  <  1:442783/250‑1 (MQ=255)
                                                                      ggTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                                                                                                                                >  1:83547/1‑219 (MQ=255)
                                                                                aaTTGGCGGAGGCGCCGTCAG‑GTGGGCTTATCCGAAGAATCCCCACAGGCCTTATAAGT‑CGCGCAGCGTCGCACCAGGCAATGTCCCCCATTTTTAACAACAAAAAAGCC‑TACTCTCTTGCAGTCGCCTTTCTCATTTTAAACGAATCACGTTTACTTCGCTTTCCCCTGGT‑TGGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                                                                                                                      <  1:560974/251‑1 (MQ=255)
                                                                                 agtgtcgcggtttgcttcggggggg‑CGTCTCCGATGAATCCCCACCGGCCTATTAGGT‑CGCGCAAGTTCTCTTCAGGCAGTGTCCCCCTTTGTTACCAAAAAAAAAGCG‑AACTCCCTTGCATCCGTCTTTTTCATTTAAAACGTATGAGGTTTCCTTCTTTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                                                                                                                      <  2:358449/227‑1 (MQ=255)
                                                                                    gTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                                                                                                                                  >  1:213362/1‑251 (MQ=255)
                                                                                     tCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                                                                                                                                 >  2:315954/1‑251 (MQ=255)
                                                                                        gATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGC‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCC‑GCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGTTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                                                                                                                              <  1:55392/251‑1 (MQ=255)
                                                                                          tGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                                                                                                                            >  1:587089/1‑251 (MQ=255)
                                                                                            cgcggtcaggtggg‑CGTAACGGATGAATCCCCACTGGCCTGATAAGC‑CGCGCGGGCTCGAATGAGGCATTGTGCTCATTTGTTAGAAAAAAAAAACCG‑GCCTCACTTGCAGTCGTCTTTCTAATTTAAAACGATTGCCGTTTCCTTCTCTTTCCCTTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                                                                                                                          <  2:121487/239‑1 (MQ=255)
                                                                                               cGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                         >  1:391691/1‑249 (MQ=255)
                                                                                               cGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                         >  1:59008/1‑249 (MQ=255)
                                                                                                    gagTGG‑CGTATACGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTATCGTCGAACTCATACACCAGCGGc                                                                                                                                                   >  2:101992/1‑250 (MQ=255)
                                                                                                     ggggg‑CGTATCCGAAGAATCACCCCAGGCTTGATAAGT‑CGCACGGCGCCGCATCAGGCAATGTGCTCCTTTGTTAGCAAAAAAAAAGCC‑GCCTCACTTGCGGTCGGTTTTCTCATTTAAAAGGAATGACTTTTCCTTCGCTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCGGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 <  2:407455/248‑1 (MQ=255)
                                                                                                     agTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 >  1:301394/1‑251 (MQ=255)
                                                                                                      gggggcgtacgcgaaaaacccccacAGCCCGGATAAGT‑CGAGCACCGTCGCTTCAGGCATGGTGCTCCATTGTTCAAAACAAAAATGCC‑AACCCACTCGCGGTCGGCTTTCTCATTTTAACCGAATGACTTTTCCTCCCCTTT‑CCCCGGT‑TTGCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 <  2:159953/230‑1 (MQ=255)
                                                                                                      gTGG‑CTTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCCCTTGCAGTCGTCTTTCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                                                                                                                                <  1:531242/251‑1 (MQ=255)
                                                                                                      gTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                                                                                                                                >  1:331305/1‑251 (MQ=255)
                                                                                                           cggagcgtatcacaccccacAGGCCTGA‑GAGTGCGCGGAGCGTCTCAGCAGGCATTGTCCTCCATGGTGAGCAAAAAAAAAGCG‑GCGTCACTGGCAGTCGGTTTTCAAATTTAAAAGGAATGACGTGAACTGCTCTT‑CCCCTGGTGTTGC‑AGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                                                                                                                            <  2:22532/235‑1 (MQ=255)
                                                                                                               cccgttaaacccccaaaGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGTAATGTCCTCCTTTGTTAGCAA‑AAAAAAGCCCGCCTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTCCTTCGTTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                                                                                                                        <  1:368846/235‑1 (MQ=255)
                                                                                                                ccGATGAATCCCCCCAGCCCTGATAATG‑CGCGCAGCTCCGCAGGAGGTCTTGTCCTCCATTGTTAGCAACAAAAAAGCC‑GCCTCCCTTGCAGCCGTCTTTCTCATTTTAAACGAATGACGTTTACTCCCCTTGCCCTTGGT‑TTTC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                                                                                                                         <  2:64363/249‑1 (MQ=255)
                                                                                                                     aaaTCACCACAGGCCTGATAAGC‑CGCGCAGCGTCGCATCAGGCAATGTTCTCCATTGTTAGCAACAAAAACGCC‑GACTCCCTTGCAGTGGGCTTTCTCATTTTAAACAAAGGACGTTTACGTCGCTTTCCCCTGTT‑TTAC‑AGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                                                                                                                                   <  1:150739/249‑1 (MQ=255)
                                                                                                                         ccgccaacaactttataaggcgccaaGCGACGCAACAGCCAATGGCCTCCATTGTTACAAACAAAAAAGCC‑GACTCACTTGCTGTCTGCTTTTCCATTTTAAACGAATCCTTTTTACTTCGCTTTCCCCTTGT‑TTGA‑ATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATATCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                                                                                                                               <  1:571210/226‑1 (MQ=255)
                                                                                                                         ccccacAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTACCAACAAAAAAGCG‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                                                                                                                              <  1:605540/249‑1 (MQ=255)
                                                                                                                           ccacAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                                                                                                                            >  1:496338/1‑251 (MQ=255)
                                                                                                                                       tAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc                                                                                                                >  1:117446/1‑251 (MQ=255)
                                                                                                                                         aGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt                                                                                                              >  2:120840/1‑251 (MQ=255)
                                                                                                                                                        gCATCAGCGAATGTGATCCCTTGTTAGCAACAACAAAGCC‑GACTCACTTGAAGTCGGCTTTCTATTTTTAAACGAATGAAGTTTACTTCGCTTTACCCGGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCGCGTCAGCATTACACAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGACGAACTCATACACCAGCAGCACGCCCGTCGGGCTATTAAGCCCAAGAATATCTTCTTCGCTCATGTTATc                                                                                                 >  2:537427/1‑250 (MQ=255)
                                                                                                                                                             aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                                                                                            >  1:285608/1‑250 (MQ=255)
                                                                                                                                                                 aaTGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                                                                                       >  1:382731/1‑251 (MQ=255)
                                                                                                                                                                          cATTGTTCGCAACAAAAAAGCC‑GAATCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTCTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCAGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGAAATTTCAACAgcgc                                                                               >  1:148167/1‑250 (MQ=255)
                                                                                                                                                                           aTTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCAc                                                                             >  1:429154/1‑251 (MQ=255)
                                                                                                                                                                                 aGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaa                                                                         >  1:445448/1‑249 (MQ=255)
                                                                                                                                                                                        aaaagGCC‑GACTCACTTCCAGTCGGTTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTAcc                                                                <  1:315954/251‑1 (MQ=255)
                                                                                                                                                                                         aaaaGCC‑GACTCATTTCCAGGCCTCTTCCTCATTTAAAAAGAATGCTGTTTCATTCTCTTCCCCCTGTT‑GTGC‑CACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCg                                                               <  2:301420/251‑1 (MQ=255)
                                                                                                                                                                                                      aCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatg                                                   >  1:442116/1‑251 (MQ=255)
                                                                                                                                                                                                        ttGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCATTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatgat                                                 >  2:471374/1‑251 (MQ=255)
                                                                                                                                                                                                             gtaggtttcctaaTTTGAAACGAATTAGGTTTCCTTGGTTTTCCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcg                                             <  1:133557/238‑1 (MQ=255)
                                                                                                                                                                                                             gTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGT‑TTGC‑ACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcgc                                            <  1:394732/251‑1 (MQ=255)
                                                                                                                                                                                                              tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcgcg                                           >  1:56429/1‑251 (MQ=255)
                                                                                                                                                                                                                gccTTTCCCATTTAAAACATATGACGTTTCCTTCCCTTTCCCCGGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCt                                         <  1:361374/249‑1 (MQ=255)
                                                                                                                                                                                                                              aaaCGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            >  1:279351/1‑250 (MQ=255)
                                                                                                                                                                                                                                aCGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg                         >  2:436287/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcgcg                       >  1:191133/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aCGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  1:52411/1‑248 (MQ=255)
                                                                                                                                                                                                                                       aCGTTTACTTCGCTTTACCCTGGT‑TTAC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  1:551782/1‑248 (MQ=255)
                                                                                                                                                                                                                                              cTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGttt           >  1:560547/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cGCTTTCCCCTGTT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCAt        <  1:322458/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGt  >  2:561377/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGG‑CGTATCCGATGAATCACCACAGGCCTGATAAGT‑CGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCC‑GACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGT‑TTGC‑AACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786611‑787105

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