Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC37.0% 51.5 / 11.2 54intergenic (‑1005/‑707)yohF/yohPputative oxidoreductase YohF/putative membrane protein YohP
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base C (7/13);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGC  >  NC_000913/2228029‑2228517
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGc                                                                                                                                                                                                                                                   >  1:182322/1‑251 (MQ=255)
   gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAGGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                                                                                                                                <  2:418155/251‑1 (MQ=255)
   gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                                                                                                                                <  2:32654/251‑1 (MQ=255)
       tGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGGGGGCCGTAACCTCCGCCCGATTc                                                                                                                                                                                                                                            >  1:290770/1‑251 (MQ=255)
            gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCCCCCGTGGTTTAACCTTAACCCCGCGGGCAGTAAACCCCGCCCCGTTCttttc                                                                                                                                                                                                                                       >  2:10030/1‑249 (MQ=255)
            gCATTAATATCGGTTTCCACACCGTCAACTTAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGAAGTAACCTCCGCCCGATTCTTTcc                                                                                                                                                                                                                                       <  2:566583/251‑1 (MQ=255)
                      cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGGAACCCCCGCCCGATTCTTTCCCCCGataacc                                                                                                                                                                                                                             >  1:148902/1‑248 (MQ=255)
                      cGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTATAACCGTTCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATc                                                                                                                                                                                                                              <  1:552672/250‑1 (MQ=255)
                            ccACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCTATAATTGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGc                                                                                                                                                                                                                      <  1:6031/251‑1 (MQ=255)
                              acacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCCCGCGTTAACCGTACCCCCGGGGGCCGGAACCCCCGCCCCATTCTTTtcccccatatccagagggc                                                                                                                                                                                                                      >  2:4519/1‑235 (MQ=255)
                                    tCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                                                                                                                                <  1:577616/249‑1 (MQ=255)
                                    tCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccg                                                                                                                                                                                                               <  2:377746/250‑1 (MQ=255)
                                       aCTGAAAAAGTATCCCCGTTCTGATATTGCTGGTGGCCGTGCGCGTGAATGGTGTTTGGTTCTTTTTTAATTGGCGGCAAGCCCGTGTTTAACCCCCCCTGGGTTTGCCCCGCCAATTACGGGTGGGCTACGCTGCCACCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTAACCCCCGCCCTATTCTGTCCCCAGAAATCGGGATGCCGGCCCGCCGt                                                                                                                                                                                                           <  2:542498/251‑1 (MQ=255)
                                                             gATATAGCTGGCGATCGGTCGCGATAGTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGACCGCCAGCAATTGGCGGGAGTACTTGCGGTCCGCCGCCCGTGGTTTAAAAGTCACCAGGCGGGGAGGAACCCACGCCCCATTCTTTACCACGATATAGGGAAGCCAGCccgccgattagaccacatgacccatgggg                                                                                                                                                                                     >  2:232735/1‑228 (MQ=255)
                                                                      ggcgatcgggctgaaTAATGGTGTTCGGTTTTTTTTCTATTTTCTGCAATCCCGTTTTTAACGCCGCTTGGGTTTGCCACTCCCATTACGGGCGGGCTCCGCTGCCAGCCAGCAATGGGCGGTTTTGCTCGCGTTCCGCCGCCCGTTCTTTAACCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTTTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACg                                                                                                                                                                              <  2:459204/235‑1 (MQ=255)
                                                                        cgaccGGTCGCGATAATGTTGTCCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTACCGCCGCCGGGGTTTGCCACTCCCAGTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGTGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                                                                                                                          <  1:231047/247‑1 (MQ=255)
                                                                            cGGTCGCGATAGGGGTGTTCTGTTCTTTTTCTTTTGCCGCCAAGACCTTGTTTAAGCCCGCCTGGGTTTTCCCCCCCCAGTACAGGCGGGCTCCGCGGCCAGCCAGCAATTGGCGGTTTTGCTCGCGTTCCGCCGCCCGTGTTTTAACCGTACCCAGGCGGGTGGTACCCTCCGCCTGATTCTTTCCCCAGATATCGAGTTGCGATCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACTGGCCGGTc                                                                                                                                                                      <  1:248909/251‑1 (MQ=255)
                                                                             ggtcgcgatatggtttttcttccttttttCTTTTGGCTGCAAGCCGGTTTTAAACGCGGCCGGGTGTTGCCCCCCCCATTACAGGGGGGCTCCGCTGCCGGCCAGAAATGGGCGGGTTTGCCCGCGTCCCCCCGCCCGTGCTTTACCCGTCCCCGGGGGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCg                                                                                                                                                                     <  2:129300/227‑1 (MQ=255)
                                                                                      aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCCTTCTTTTAACCTACCCCCGCGGGCAGTAAACTCCGCCCGATTCTTTTCCCACATCTCGAGATGCCCTCCCGCCCTTAAACCCAAAATACCGTTGGGGTACCACGGGCCCGTCGTACCCGCgg                                                                                                                                                            >  2:452994/1‑251 (MQ=255)
                                                                                      aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTTGCGTTTCGCCCCCCGCGCTTTAACCGGATCCCGGCGGGCAGCAACCCCCCCCCGATTCTTTCCCCAGAGATCGAGATTGCCACCCGCCCGTACACCAAAAGTAACGTTGGTGTAcccccagcccgtctcacatccgg                                                                                                                                                            >  1:38074/1‑230 (MQ=255)
                                                                                      aaTGGTGTTCGGCCTTTTTTCTATTTGCTGCAAGCCCGTGTTCAACGCCCCCGGGGTTGGCCCCCCCCAGCACGGGCGGGCTCCGCCGCCAGCCAGCAATGGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCATATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                                                                                                                             <  2:305278/250‑1 (MQ=255)
                                                                                          gtttcttGCTTTTTTTCTATTTGCTGCAAGCCCGTTTTACCCCCCGCCGGGGTTTCCCACCCCCGTTCCAGGGGGGGTCCGTGGCCGGCCAGCATTTGGCGTTTTTCCCCGCTTCCCGCCCCCTGGGCTTTAACCTTCCCCGGGGGGGCAGTCACCCCCGCCCGATTCTTCCCCCAGATAGCGAGATGCGAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCgg                                                                                                                                                        <  2:553701/245‑1 (MQ=255)
                                                                                             ttcttcgttttttCTATTCGGTGCAAGCCCGTGTTTACCCGCGCTTGGTTTTGCCCCCCCCAGGAGGGGGGTGCTCCGCTGCCACCCAACAATTGGGGGGTTTGCTCGCGTTCCGCTGCCCTTGTTTTCACCATCCCAAGGCGGGCAGAACCCTCCCCCCAATTCTCTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                                                                                                                                      <  2:12324/242‑1 (MQ=255)
                                                                                             ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGGGCTTTAACCCGACCCAGGCGGGGAGTAAACTCCGCCCGATTTTTTCCCCAGATATCTGGATGCCAGCCCGCCGTTAAGCCAAAAAGACCGTTTGGGTACCACGGGCCGGGCGTAACTGCGGCCGGATc                                                                                                                                                     >  1:426398/1‑251 (MQ=255)
                                                                                                  ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCCCCCCCCGCGCTTTTACCCTACCCCGGCGGGCCGTAACCTCCCCCCCATTCTTTTCCCCCATATCGAGATGCCCGCCCCCCCTTTAGCCAAAAATACCGTTTGTGGTCCCCCGGCCGGGCGTACCCTCGGCCGGacccttc                                                                                                                                                 >  2:185515/1‑244 (MQ=255)
                                                                                                       ttttctattggaggaaaGACCGTGTTAACCGCCGCGTGGGTTGGCCACTCCCAGTAAAGGGGGGCGACGCGGCAAGCCAGCAAGTGGCGGGTGGTCTCGGGGTCCGCCGCCCGTGCTTTAACCTTACCCAGGGGGGCAGAACCCTCCGCCCGATTCTTCCCCAAGATATGGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                                                                                                                              <  1:153739/233‑1 (MQ=255)
                                                                                                         tCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCCCCCGTGCTTTAACCCTACCCAGGCGGGCAGGAACCCCCGCCCGATTCTTTCCCCCGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTTGTGGACCACGGGCCGGGCGTCCCTGCGGCCGGATTGTTCAGAGCaa                                                                                                                                         >  1:556118/1‑250 (MQ=255)
                                                                                                            ttgttcttgaagctccGTGTTGAGCGCCGTGTGGGTTTGCCCCTCCGAGTACTGGCGGGCTACGCTGCCCTCCGCCAGTTGGCGGTTTGTCTCGCTTTCCGCCGCCCGTGTTTTACCCGTACCCAGGCGGAAATTAACCTCCGCCCGATTCTTTCCCGTAATATCGATATGCTGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGTGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                                                                                                                       <  2:307185/236‑1 (MQ=255)
                                                                                                            aTTTGCTGCAAGACCGTGTTTAACGCCGCCGGGGTTTGCCCCCCCAATGACGGGGGGCCTCCGCTGCCGGCCAGCAATTGGCGGGTTTGTTCGCGTTCCGCCCCCCGTGCTTAACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATAGCGAGATGCCATCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                                                                                                                       <  2:177446/249‑1 (MQ=255)
                                                                                                                gctgcAAGACCGAGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTGCTCCCGTTCCGCCGCCCGCGCCTTAACCGTACCCCAGCGGGGCGTAACATCCGCCCCATTCTTTTCCCCGATAACGAGATTCCCGCCCCCCGCTAAACCCAAAATACCGGTGGAGTAACACGGGCCGGGCGGACCCGCGGCCGGGTCGTGTCACAGAAACGGcccc                                                                                                                                 >  2:479970/1‑251 (MQ=255)
                                                                                                                gcggcAAGACCGTGTTGAACGCCGCCGGGGTTTGCCCCTCCCAGTACGGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCCCCCGTGCTTTACCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                                                                                                                                 <  2:202865/248‑1 (MQ=255)
                                                                                                                 ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCCGGCGTACCTGCGGCCGGATCGGTCCGAGCAAAAGGccccc                                                                                                                                >  1:466749/1‑250 (MQ=255)
                                                                                                                       agcccGTGTTTACCCCCGCCTGGGTTTGCCACCCCCAGTACGGGCGGGCTACGCTGCACGCCAGCAAGTGGGGGGTGTGCCCCCGTCCCCCCCCCCGTTCTTTACCCGTCCCCAGGCGGTCATTAACCCCCGCCCGTTTCTTCCCCAAGATCTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                                                                                                                          <  2:530848/248‑1 (MQ=255)
                                                                                                                        gggtgtgatgatattctcccTGGGTTTGCCCCTCCCAGTCCGTGCGGCCTACGCTGCCAGCCAGAAGTGGGCGGGTTTTCTCGCGTCCCGCCGCCCGGGCTTTACCCGTCCCCGGGCGGGAAGTACCCCCCGCCCGTTGCTTCCCCCAGATACAGATATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                                                                                                                          <  2:541409/233‑1 (MQ=255)
                                                                                                                         cccGTGTTTAACGCCGCTTGGGTTTGCCCCCCCCAGTACGGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCTTCCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                                                                                                                        <  1:10030/250‑1 (MQ=255)
                                                                                                                                   accgcccccggggTTTGCCACCCCCAGTACAGGCGGGATCCGCGGCCAGCCAGAATTTGGCGGGTTTGCCCTCGTTCCCCCGCCCGGTCTTTACCCGTCCCCGGGCGGGCTGTACCCCCCCCCAGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGa                                                                                                               <  2:581457/241‑1 (MQ=255)
                                                                                                                                           cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCCTTAACCGTAACCCGGCGGGCAGTAACCCCCGCCCGAATTTTTCCCCCGGTATCGAGATGCCCGCCCCCCGTTTAGCCAAAAATACCGGTGGGGTACCCCCGGCCGGTCGCACCTGCGGGCGGATCGTTCCAAGCAAACGGCCCCCATTAGCCTTTTGCCGACCtttttt                                                                                                       >  1:520238/1‑251 (MQ=255)
                                                                                                                                           cTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGCCCGTGCTTTTAACCTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCACAAGTACCGTTTGTGTACCCCGGGCCGGGCGTACCTGCGGCCCGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttttt                                                                                                       >  1:259735/1‑251 (MQ=255)
                                                                                                                                                         ccccAGTACAGGCGGGTTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATc                                                                                          <  1:144949/250‑1 (MQ=255)
                                                                                                                                                            cggTACAGGCGGGCTACGCTGCCAGCCAGCAATGGGGGGGTTTGCTCGCGTTCCGCCGCCCGTTCTTTAACCGCCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATccgc                                                                                       <  1:322608/249‑1 (MQ=255)
                                                                                                                                                                      cgggctaggtggccatccatcatgtGGCGGGTTGGGGCGCGTGCCGCCGCCCGTGCTTTAACCGTCCCCAGGTGGGTAGGAACCCCCGCCGGATTCTCCCCCCTAATAGCGAGATGCCAGTCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGTAGCGTGCAGAGCAAACGGCCCCATCAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTc                                                                              <  1:500886/227‑1 (MQ=255)
                                                                                                                                                                           tACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCCTTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGCAAAGTCGATCt                                                                         >  1:288396/1‑251 (MQ=255)
                                                                                                                                                                               cTGCCAGCCAGCAATTGGCGGGTTTGCTCGCCTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCCCCCGGTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCCTTGGGGTACCACGGGCCGGTCGTACCTGCGGCCGGGTCGTTCTGAGCAAACGGCCCCACTAAGCCTTCTGCCGACCTTTTTTTCCGCTCCCTATCCCCCGAAAAGACGCTCTGccg                                                                     >  1:13531/1‑249 (MQ=255)
                                                                                                                                                                                              tggggggTTTGCTCCCGTCCCCCCCCCCGTGCTTTACCCGTCCCCAGGGGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTgg                                                       <  1:4519/246‑1 (MQ=255)
                                                                                                                                                                                                                   cgccgcCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCtt                                  >  1:57294/1‑251 (MQ=255)
                                                                                                                                                                                                                        cccGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagct                             >  1:508569/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggggttttAACCGTACCCGGGCGGGCAGTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagctag                           <  1:452994/246‑1 (MQ=255)
                                                                                                                                                                                                                             gCTTTAACCGTACCCAGGCGGGCAGTAACCTCAGCCCGATTCTTTCCCCAGCTATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTg                        >  2:55379/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cccgatccccAGGCGGGCATTACCCCCCGCCCGGTTCTTTCCCCAGATTGAGAGGTGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc                   <  2:132867/244‑1 (MQ=255)
                                                                                                                                                                                                                                   aaCCGTCCCCCGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATTTCGAGATTCGAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgc                   <  1:278970/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCt                >  1:2368/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccAGGCGGGAAGTATCCTCCGCCCGATTCTTCCCCCAGAAAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCa           <  1:15839/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGACGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc     >  1:370560/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     aGTAACCTCAGACCGATTCTTTCACCAGGTATCGAGAAGCCAGCCCGCCGTTACGCAAAAAGAGCCGTTGGAGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGGCATTTTTTGCCGCTCCATATCCGACGACAAGTCGAGCTGAGGACCATCCTGGGTGGCAAATGCCTTCGCCTGGGCTTCAGCTAGCTGAAGGCGCTGTTCAGCCACCTGc  >  2:565432/1‑250 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGC  >  NC_000913/2228029‑2228517

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