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→GGGgpmA2,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.

GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  NC_000913/786611‑787105
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
gCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTgg                                                                                                                                                                                                                                                        <  1:527675/251‑1 (MQ=255)
                                gttacgaaaaCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTAAGAGTGGCGTATCCGATAAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGAAACAAAAAAGCCGACTCTCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                                                                                                                       <  1:328186/242‑1 (MQ=255)
                                gTTACTTCAACTGGAGAATGCAGTAAAAGCGCGTCGCGGGTAGGGTTTAATGGTCGGATGCGCCGCCAGAGGGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTGGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                                                                                      <  1:225234/251‑1 (MQ=255)
                                   aCGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                                                                                                                                   >  1:301420/1‑251 (MQ=255)
                                         gcTGGCGAATGCAGAAACAGCGCGGCGCGGGTAAGGTGTAATAGGCGGATGCGCCGAGAGGGGGGCGTATCCGATGAATCACCACTGGCCTGATAAGTCGCGCAGCGCCGCTTCAGGCAATGTGCCCCATTGTTAGAAAAAAAAAAGCCGACCCACTTGCAGTCCGCTTTCTCTTTTAAAACGAATTTCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAtt                                                                                                                                                                                                             <  2:143712/250‑1 (MQ=255)
                                           gggggaaagCAGAAAAAGCGCGCCGGGGGTAAGGTTTAATTGCGGGATGCGCCTACAGGGTGGCGTATCCGATGATACACCACAGGCCTGAAAAGTCGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTGGCAGTCGTCTTTCTCATTTAAAACGAATGACGTTTTCTTCGTTTTACCCTGGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                                                                                                                           <  1:354249/243‑1 (MQ=255)
                                            ggggAATGCAGAAAATGGGCGTCGCGGGAAAGTTTTAATGGTGGGATGCGCCGGCAGAGTGGCGTACCCGAAGAATCACCACAGGCCGGAAAGGCCGCGCAGCGTCGCATCGGGAAGGGTGCTCCATGGTTAGAAAAAAAAAAGCCGCCTCCCTTGCAGTCGTCTTTCTCATTTAAAACGAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                                                                                                                          <  2:287654/248‑1 (MQ=255)
                                              agAATGCAGTAAAAGCGCGTCGCGGGTAAGGGTTAATTGTCGGATGCGCTGTCAGAGTGGCGTATCCGATGAATAACCACAGGCCTGATAAGTCGCGCAACGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGACGACACACTTGAAGTCGGCTTTCTCATTTTAAACGGATGAAGTTTACTTCGCTTTACCCTGGCTTGCAACCGCGGATGCTTTCACTGCTCTCTGGTCCGCATTAccc                                                                                                                                                                                                         >  2:302134/2‑250 (MQ=255)
                                               gAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGACGGATGCGCCGTGAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCATCATTAccc                                                                                                                                                                                                         >  2:112834/1‑249 (MQ=255)
                                                aaTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                                                                                                                         >  2:542457/1‑248 (MQ=255)
                                                    aagaaaaagggggTCGCGGTTAGGGTTAAATTGTCGGATGCCCCGTCAGAGTGCCTTACCCGTTAAACCCCCACAGGCCTGATAAGTGGCGCAGCTCCGCATCAGGCAATGTGCTCCATTGTTATCAACAAAAAAGCAGCCTCCCTTTCATTCGGCTTTCTCATTTTAAACGAATGACTTTTTTTTCGCTTTTCCCGGGTTTTCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                                                                                                                                  <  1:547050/239‑1 (MQ=255)
                                                      gTAAAAGCGCGTCGCGGGTAAGGTTTACTTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGGTAGCAACAAACAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCGGCTTTCGCTTCGATCTCGTCAGCATTAACCAGataat                                                                                                                                                                                                  >  2:327877/1‑249 (MQ=255)
                                                         aaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                                                                                                                             >  2:4126/1‑251 (MQ=255)
                                                          caGCGCGTCGCGGGTAAGGATTTATTGTCGGCTGCGCCGTCCGATTGGCGTATCCGCTGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACACAAAAACCTACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATTACGATTACTTCTCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATATCGTCAACATTACCCCGATAATCGCGtt                                                                                                                                                                                            >  2:382097/2‑251 (MQ=255)
                                                             cgcgTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTcc                                                                                                                                                                                         >  1:182839/1‑250 (MQ=255)
                                                                 gcgcgGGAAAGTTTTAATTGTGGGAGGCGCCGGCAGAGTGGCGTATCCGAAGAACCCCCACGGGCCTGAAAAGTCGCGCAGCGTCGCATCGGGAAATGTGCCCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTCCGTTTTCCCCTGTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                                                                                                                                     <  2:484049/250‑1 (MQ=255)
                                                                  ggctggTAAGGTTTATTTGCCGGATGGGCCGTCAGAGTGGCGTAGCCGATGAATCACCACAGGCTTGATAAGCCGCGAAGGGTGGCACCGGGAAAGTTGCTCATTTGTAAGAAAAAAAAAAGCCGACTCATTTGCATGCGGTTTTCTCTTTTTAAACGAATGGCTTTTCCTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGGTTCAGCGGt                                                                                                                                                                                    <  1:409036/247‑1 (MQ=255)
                                                                  cgcgGTTAAGGTTTAATGGGCGGAGGCGCCGTCAGGGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGCCGCGCGGCGTCGCACCAGGAAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                                                                                                                                    <  1:100530/251‑1 (MQ=255)
                                                                    ggggTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCCCCACAGGCCTGATAAGTCGCGCAGCGTCGCACCAGGCAATGGCCTCCATGGTTAGCAACAAAAAAGCCGACTCACTTGCAGGCGGCTTTCTCATTTAAAACAAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                                                                                                                                  <  1:442783/250‑1 (MQ=255)
                                                                      ggTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                                                                                                                                >  1:83547/1‑219 (MQ=255)
                                                                                aaTTGGCGGAGGCGCCGTCAGGTGGGCTTATCCGAAGAATCCCCACAGGCCTTATAAGTCGCGCAGCGTCGCACCAGGCAATGTCCCCCATTTTTAACAACAAAAAAGCCTACTCTCTTGCAGTCGCCTTTCTCATTTTAAACGAATCACGTTTACTTCGCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                                                                                                                      <  1:560974/251‑1 (MQ=255)
                                                                                 agtgtcgcggtttgcttcgggggggCGTCTCCGATGAATCCCCACCGGCCTATTAGGTCGCGCAAGTTCTCTTCAGGCAGTGTCCCCCTTTGTTACCAAAAAAAAAGCGAACTCCCTTGCATCCGTCTTTTTCATTTAAAACGTATGAGGTTTCCTTCTTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                                                                                                                                      <  2:358449/227‑1 (MQ=255)
                                                                                    gTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                                                                                                                                  >  1:213362/1‑251 (MQ=255)
                                                                                     tCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                                                                                                                                 >  2:315954/1‑251 (MQ=255)
                                                                                        gATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGCCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGCCTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                                                                                                                              <  1:55392/251‑1 (MQ=255)
                                                                                          tGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                                                                                                                            >  1:587089/1‑251 (MQ=255)
                                                                                            cgcggtcaggtgggCGTAACGGATGAATCCCCACTGGCCTGATAAGCCGCGCGGGCTCGAATGAGGCATTGTGCTCATTTGTTAGAAAAAAAAAACCGGCCTCACTTGCAGTCGTCTTTCTAATTTAAAACGATTGCCGTTTCCTTCTCTTTCCCTTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                                                                                                                          <  2:121487/239‑1 (MQ=255)
                                                                                               cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                         >  1:391691/1‑249 (MQ=255)
                                                                                               cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                                                                                                                         >  1:59008/1‑249 (MQ=255)
                                                                                                    gagTGGCGTATACGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTATCGTCGAACTCATACACCAGCGGc                                                                                                                                                   >  2:101992/1‑250 (MQ=255)
                                                                                                     gggggCGTATCCGAAGAATCACCCCAGGCTTGATAAGTCGCACGGCGCCGCATCAGGCAATGTGCTCCTTTGTTAGCAAAAAAAAAGCCGCCTCACTTGCGGTCGGTTTTCTCATTTAAAAGGAATGACTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCGGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 <  2:407455/248‑1 (MQ=255)
                                                                                                     agTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 >  1:301394/1‑251 (MQ=255)
                                                                                                      gggggcgtacgcgaaaaacccccacAGCCCGGATAAGTCGAGCACCGTCGCTTCAGGCATGGTGCTCCATTGTTCAAAACAAAAATGCCAACCCACTCGCGGTCGGCTTTCTCATTTTAACCGAATGACTTTTCCTCCCCTTTCCCCGGTTTGCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                                                                                                 <  2:159953/230‑1 (MQ=255)
                                                                                                      gTGGCTTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGTCTTTCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                                                                                                                                <  1:531242/251‑1 (MQ=255)
                                                                                                      gTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACg                                                                                                                                                >  1:331305/1‑251 (MQ=255)
                                                                                                           cggagcgtatcacaccccacAGGCCTGAGAGTGCGCGGAGCGTCTCAGCAGGCATTGTCCTCCATGGTGAGCAAAAAAAAAGCGGCGTCACTGGCAGTCGGTTTTCAAATTTAAAAGGAATGACGTGAACTGCTCTTCCCCTGGTGTTGCAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                                                                                                                            <  2:22532/235‑1 (MQ=255)
                                                                                                               cccgttaaacccccaaaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGTAATGTCCTCCTTTGTTAGCAAAAAAAAGCCCGCCTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTCCTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                                                                                                                        <  1:368846/235‑1 (MQ=255)
                                                                                                                ccGATGAATCCCCCCAGCCCTGATAATGCGCGCAGCTCCGCAGGAGGTCTTGTCCTCCATTGTTAGCAACAAAAAAGCCGCCTCCCTTGCAGCCGTCTTTCTCATTTTAAACGAATGACGTTTACTCCCCTTGCCCTTGGTTTTCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                                                                                                                         <  2:64363/249‑1 (MQ=255)
                                                                                                                     aaaTCACCACAGGCCTGATAAGCCGCGCAGCGTCGCATCAGGCAATGTTCTCCATTGTTAGCAACAAAAACGCCGACTCCCTTGCAGTGGGCTTTCTCATTTTAAACAAAGGACGTTTACGTCGCTTTCCCCTGTTTTACAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                                                                                                                                   <  1:150739/249‑1 (MQ=255)
                                                                                                                         ccgccaacaactttataaggcgccaaGCGACGCAACAGCCAATGGCCTCCATTGTTACAAACAAAAAAGCCGACTCACTTGCTGTCTGCTTTTCCATTTTAAACGAATCCTTTTTACTTCGCTTTCCCCTTGTTTGAATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATATCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                                                                                                                               <  1:571210/226‑1 (MQ=255)
                                                                                                                         ccccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTTTTACCAACAAAAAAGCGGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                                                                                                                              <  1:605540/249‑1 (MQ=255)
                                                                                                                           ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                                                                                                                            >  1:496338/1‑251 (MQ=255)
                                                                                                                                       tAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc                                                                                                                >  1:117446/1‑251 (MQ=255)
                                                                                                                                         aGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt                                                                                                              >  2:120840/1‑251 (MQ=255)
                                                                                                                                                        gCATCAGCGAATGTGATCCCTTGTTAGCAACAACAAAGCCGACTCACTTGAAGTCGGCTTTCTATTTTTAAACGAATGAAGTTTACTTCGCTTTACCCGGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCGCGTCAGCATTACACAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGACGAACTCATACACCAGCAGCACGCCCGTCGGGCTATTAAGCCCAAGAATATCTTCTTCGCTCATGTTATc                                                                                                 >  2:537427/1‑250 (MQ=255)
                                                                                                                                                             aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                                                                                            >  1:285608/1‑250 (MQ=255)
                                                                                                                                                                 aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                                                                                       >  1:382731/1‑251 (MQ=255)
                                                                                                                                                                          cATTGTTCGCAACAAAAAAGCCGAATCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTCTACCCTGGTTTGCAACCGCCGCTGCTTTCGCAGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGAAATTTCAACAgcgc                                                                               >  1:148167/1‑250 (MQ=255)
                                                                                                                                                                           aTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCAc                                                                             >  1:429154/1‑251 (MQ=255)
                                                                                                                                                                                 aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaa                                                                         >  1:445448/1‑249 (MQ=255)
                                                                                                                                                                                        aaaagGCCGACTCACTTCCAGTCGGTTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTAcc                                                                <  1:315954/251‑1 (MQ=255)
                                                                                                                                                                                         aaaaGCCGACTCATTTCCAGGCCTCTTCCTCATTTAAAAAGAATGCTGTTTCATTCTCTTCCCCCTGTTGTGCCACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCg                                                               <  2:301420/251‑1 (MQ=255)
                                                                                                                                                                                                      aCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatg                                                   >  1:442116/1‑251 (MQ=255)
                                                                                                                                                                                                        ttGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCATTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCgatgat                                                 >  2:471374/1‑251 (MQ=255)
                                                                                                                                                                                                             gtaggtttcctaaTTTGAAACGAATTAGGTTTCCTTGGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcg                                             <  1:133557/238‑1 (MQ=255)
                                                                                                                                                                                                             gTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcgc                                            <  1:394732/251‑1 (MQ=255)
                                                                                                                                                                                                              tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCAcgcg                                           >  1:56429/1‑251 (MQ=255)
                                                                                                                                                                                                                gccTTTCCCATTTAAAACATATGACGTTTCCTTCCCTTTCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCt                                         <  1:361374/249‑1 (MQ=255)
                                                                                                                                                                                                                              aaaCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCAt                            >  1:279351/1‑250 (MQ=255)
                                                                                                                                                                                                                                aCGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcg                         >  2:436287/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATAcgcg                       >  1:191133/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aCGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  1:52411/1‑248 (MQ=255)
                                                                                                                                                                                                                                       aCGTTTACTTCGCTTTACCCTGGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGc                     >  1:551782/1‑248 (MQ=255)
                                                                                                                                                                                                                                              cTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGttt           >  1:560547/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCAt        <  1:322458/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGt  >  2:561377/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
GCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAAGAGTTACCGTGTGCAGCGATGATCACGCGCTCACCGCTCTTCATACGCGGCAGAATAGTTTCATTCCAGT  >  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: