Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,611,4550GA21.1% 108.3 / 10.6 54C93C (TGC→TGTyhhSputative transporter YhhS
Reads supporting (aligned to +/- strand):  ref base G (8/33);  new base A (11/0);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 1.25e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGTGATAAAACACGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGTTT  >  NC_000913/3611207‑3611704
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
cAGTGATAAAACACGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAagca                                                                                                                                                                                                                                                             <  1:89334/251‑1 (MQ=255)
      tAAAACACGAAGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTTATCCTGCCCGCCGATACCCCCGAACGCTCCAAAagaagacgcc                                                                                                                                                                                                                                                       >  1:104517/1‑243 (MQ=255)
                      aCGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGACATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCAAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTCATCCTGCCGTCACATAACCCCCACCGCGTAAAACACAGCCGCATAAACAGAAgacgac                                                                                                                                                                                                                                       >  1:480100/1‑250 (MQ=255)
                               ccagctcaatcgccCCGTAAGTTACAATGCCGGTCCACGAAATCACCCGCCCGATATGCAGCGAGCCATCCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTGTGCCCAATCCCAAGGACGACGCGCCCCAGGCTGAGTAATAACATGCTGATGGCAGTCAGTCTGGCTGTTAATCCTGCCGGCAGAGGCCCCAGGCCTCTCAACAATCAGCCGCATAAACCAAGGACGACGATCTTTTTg                                                                                                                                                                                                                            <  2:462556/243‑1 (MQ=255)
                                catcgccatcgcCCCGTAAGTTACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTTGCGGTTAATCCTGCCGTCAGATACCCCAGATCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTgg                                                                                                                                                                                                                           <  1:486993/251‑1 (MQ=255)
                                  ccgccatcgccCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGt                                                                                                                                                                                                                         <  2:138359/250‑1 (MQ=255)
                                                     gACAATGCCGTTCCACGACCTCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTTCCGTCAGAAACCCCCGACCCCTCCAAAAACAGCCGCCTAAAACGAAGACGACAACCTTTTTTGGGCCCAGCGAATCGGCggaa                                                                                                                                                                                                       >  1:373293/1‑247 (MQ=255)
                                                             cGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGGCCCAGCGAATCGGCGTAACGTCCGGCa                                                                                                                                                                                              <  2:312395/251‑1 (MQ=255)
                                                              gTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACACCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCATACTGGCTGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAATCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCAt                                                                                                                                                                                             <  1:335489/251‑1 (MQ=255)
                                                                 ccACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCAATCGGGTCGATCCCGTTCCGGCCAAACTTTGCCCAATCCCAAGGATGAAGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGATAATCCTGCCGCCAGGTTACCCCGACCGCTCAAAAAGGAGCCGAAAAAAACGAAAACGGCAAGTTTTTTGGGGCCCCGCGacacgggggaaaggccggcatg                                                                                                                                                                                            >  1:461642/1‑228 (MQ=255)
                                                                  cccGAAATCCCCCGCCCGATATGCATCGAGCCACCCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAACCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGTCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTGGGGTGCCAGCGAATCGTCGTAACGTCCGGCATGAgg                                                                                                                                                                                         <  2:42223/249‑1 (MQ=255)
                                                                             ccgcccgATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCAGAAACCGAAGACGACAATCTGTTTGGGTCCCAGCTAATCGGCTTAACGGCCTGCATGAGGGCGGCTCAGCa                                                                                                                                                                              <  1:417855/251‑1 (MQ=255)
                                                                                   gATATGCAGCGAGCCAACCCCGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTgg                                                                                                                                                                        <  2:261795/251‑1 (MQ=255)
                                                                                    atatGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAg                                                                                                                                                                                <  1:113911/242‑1 (MQ=255)
                                                                                             gAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCCGCTACCCCAGACCGCTCCAAAAACAGCCGCATAAACCGAAGAAGACAATCTTTTTTGGTCCCCGCGAAACGGCGGAACGTCCCGCATGGGGGGGGGTCAGCCAGGGGGCGAAATAttt                                                                                                                                                              >  1:153407/1‑250 (MQ=255)
                                                                                              aGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCATCAAGGTGGCGAAATATTGc                                                                                                                                                             <  1:177407/251‑1 (MQ=255)
                                                                                                aacgccACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGAGGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCTAAGACGACAAGCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGTTGGCGAAATATTGCAg                                                                                                                                                           <  2:281650/247‑1 (MQ=255)
                                                                                                   aCCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGGGGCGAAATATTGCAGGCt                                                                                                                                                       <  1:379447/251‑1 (MQ=255)
                                                                                                    ccACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAATGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTg                                                                                                                                                      <  1:249712/251‑1 (MQ=255)
                                                                                                    ccACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAAATTTGCCCAATCCCAAGGATGACGAGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTTCCCCCCGATACCCCAGACCGCCCCAAAAACAGCCGCATAAACCGCAGAACACCAACTTTTTTGGTCCCAGCGCAACGGCGGAACGTCCGGCACGAGGGGGGCTCAGCAAGGTGGCGGAATAtttgcggg                                                                                                                                                         >  1:281650/1‑243 (MQ=255)
                                                                                                          cAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCCGAAAACCCCGACCCCCCAAAAAACAGCCGCCTACAACGCAGACGACAATCTTTTTTGGGTCCCGCGAATCGGCGTAAAGTCCGGCATGAGGGCGGCTCAGCACGGGGGCGAAAAAATGCAGGCTGATAAcc                                                                                                                                                >  1:352571/1‑251 (MQ=255)
                                                                                                           aaCGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGGGGCGAAATATTGCAGGCTGATAACCa                                                                                                                                               <  1:122297/251‑1 (MQ=255)
                                                                                                           aaCGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGCATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCCGCCGTCAGATAATCCCGACCCCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTTGGACCCAACGAATCAGCGGAACATCCGCGATTAGGGGGGCGCCGCAAGGTGGCGGAatcttttcggctggtaacc                                                                                                                                                 >  1:329860/1‑233 (MQ=255)
                                                                                                           aaCGCCCCATAGGGTCGATCCCGTTCAGGCAAAACTTTGCCCCATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACCGGCAGACTGGCGGTTTATCCCGCCGTCAGATTACCCAGACCGCTCCAAAAACAGCCGGCTAAAACGAAGACGTCAATCTTTTTGGGTGCCAGCGAATCGGCGGAAAGTCCGGGCTGTGGGCGGGCCAGCAAGGGGGCGGCAAATTGCCGGCTGGTAAccc                                                                                                                                               >  1:410732/1‑250 (MQ=255)
                                                                                                                cccATAGGGTCGATCCCGTTCCGGCAAAGCTCTGCCCAATCCCAAGGATGACGCGCCCCAGGCCAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCTGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCATTTTGGGGCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGGAATATTGCAGTCTGATAACCCATCCt                                                                                                                                          <  2:332681/251‑1 (MQ=255)
                                                                                                                cccATAGGGTCGATCCCGTTCCGGCAAAACGTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGATGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCt                                                                                                                                          <  1:32793/251‑1 (MQ=255)
                                                                                                                cccAGAGGGTCGATCCCATTCCGGCAAAACTTTGCCCAATCCCAAGGATGACCCGTCCCAGGCAAAGTAATAACAGGCGGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCACAAAAAGCAGCCGCATAAGCCGAAGACGACAATCTTTTTGGGTCCCTGCGAATCGGCGTAACGTCCGGCAGGAGGGCGGCTCAGCGTGGTGGCGAGATATTGCAGGCTGATAACCAATCCt                                                                                                                                          <  2:85813/251‑1 (MQ=255)
                                                                                                                 ccATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGGCCGATAACCCAGACCGCCCAAAAAACAGCCCCAAAAACCGAAGAAGACCATCTTTTTGGGTACCACCGAATCGGAGTCAAGGCCGGCATGAGGGGGGGTCAGCAAGGTGGCGGAATATTTCAGGCTGATAACCAAACCtt                                                                                                                                         >  1:29345/1‑250 (MQ=255)
                                                                                                                    caggggccaTCCCGATCCGGCAAAACTTTGCCCAATCCCAAAGATGACGCGCCCCAGGCACAGTAAGAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGGCAGAGACCCCAGACCGCGCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCGAGCGAATCGGCGGAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGccc                                                                                                                                     <  2:318126/243‑1 (MQ=255)
                                                                                                                     aGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGccc                                                                                                                                     <  1:387561/250‑1 (MQ=255)
                                                                                                                                 gTTCCGGCAAAACGTTGACCAAACCCAAGGAGGACGCGCACCAGGCAAAGTAATATCAGGCTGATGACAGGCGGGCTGCGGGTTAAACCTACCGACAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAGGCCGAAGACGACAAGCTTTTTGGGACGCTGCGTATCAGCGTAACGTCGGGCATGAGGGCGGCTCTGCAAGTTCGCTAAAAATAGCAAGCTGATAACCAAGCCTGGCCAGAAGGCTCTaaa                                                                                                                        <  2:252737/251‑1 (MQ=255)
                                                                                                                                    ccGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAata                                                                                                                                                                 <  1:236360/209‑1 (MQ=255)
                                                                                                                                    ccGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGc                                                                                                                      <  2:373293/250‑1 (MQ=255)
                                                                                                                                           aaCTTTGCCCCATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCcatcaca                                                                                                              <  2:267359/251‑1 (MQ=255)
                                                                                                                                           aaCTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTGCCGTCAGATACCCCAAACCGCTCAAAAAACAGCCGCCTAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAAACGGCGTAAAGTCCGGCATGAGGGGGGCCCCGCAAGGGGGGGGAAAATTGCAGGCTGGTAACCAATCCTGCCCAGAAGGCGCTAAAGccccacac                                                                                                               >  1:1532/1‑245 (MQ=255)
                                                                                                                                                         ccAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGATAATCCCGCCCTCAGGTACCCCAAACCGCCCAAAAAGCAGCCGCCAAAAACGAAGACAACCATCTTTTTTGGTCCCAGCGAAACGGCGGACCGTACGGCATGAGGGCGGCCCAGCAAAGTGGCGAAAAATTTCAGGGTGATTAACAATCCTGACCACAAGACGCTCAAACCCATCACCTCATCGACATAGcc                                                                                                >  1:485909/1‑251 (MQ=255)
                                                                                                                                                           aaGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCg                                                                                              <  1:68243/251‑1 (MQ=255)
                                                                                                                                                              gATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTGCCGGCAGATACCCCCGGCCGCTCAAAAAACAGCCGCCTAAACCGAAGGCGACAATCTTTTTGGGTCCCCGCGAAACCGCGGAACGTCCGGCCTGAGGGGGGCCCCGCAAGGGGGGGAAATATTTCAGGCTGAATACCAATCCTTCACCCGAGGGGCTAAAGCCcatcaccatatgggcctaacccggg                                                                                            >  1:336608/1‑231 (MQ=255)
                                                                                                                                                                      gcCCCAGGCAAAGTAATAACAGGCGGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACagcg                                                                                   <  2:336608/251‑1 (MQ=255)
                                                                                                                                                                            ggCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGc                                                                             <  1:93481/251‑1 (MQ=255)
                                                                                                                                                                            ggCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCACCCCCCCAAACACCCCAGAACGCTCAAAAAACAGCCCCATAAACCGAAGGACACCAACTTTTTTGGTTCCCGCGAATCGGCGTAACGCCCGGCCTGGGGGGGGCCCCGGAAGGTGGCGCAATATTGTAGGCTGACAAACAATCACGCCCACAAGGGGCCAAAACCCACCCCATCATGGGCAAAGCTCGGGAATAAAGCGAGCGGc                                                                             >  2:159383/1‑251 (MQ=255)
                                                                                                                                                                             gCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAAGCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCa                                                                            <  2:50827/251‑1 (MQ=255)
                                                                                                                                                                              cAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTGCCGTCAGGTAACCCCGAACGCCCAAAAAACAGCCGGATAAAACGGAGAAGACAATATTTTTTGGTGCCAGCGAAACGGCGTAAAGTCCCGCCAGGGGGGGGGTCAGCGAGGGGGCGGAATATTTCAGGCTGGTAAACAAACCTGCCCCGAAGGGGCCAAAATACATCCCATTATCGACAGAGCCCGGTAAttccgagaggggcg                                                                            >  2:246640/1‑237 (MQ=255)
                                                                                                                                                                               aaaGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAAc                                                                          <  1:188425/251‑1 (MQ=255)
                                                                                                                                                                                    aataaCAGGCTGATGACAAGCAGACGTGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGCCCCAGCGAATCGGCGTAACGTCCGGCAGGAGGGCGGCTCAGCAATGTGGCGAAATATTGCAGGCTGATAACCAATCCAGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATATCCCGGTAAGACATCGAGCGTCAACCCGAt                                                                     <  2:461642/251‑1 (MQ=255)
                                                                                                                                                                                      taaCAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGCTACCCCCGACCGCTCCAAAAGCAAACGCATAAACCCAAGACGAAAATCTTTTTTGGTCCCCGCGAATAAGCGGAAAGTCCCGCCTGAGGGGGGGCCCGCAAGGGGGCGGAATATTGCAGAACGGAAAACAATCCCGCCCAGAAGGCGCCACAGCCCATCACATCATGGACATAACCAGGCACTACagcgcgccgaaaaccgctgg                                                                   >  1:448152/1‑235 (MQ=255)
                                                                                                                                                                                          aGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAGCCCGATGGTGAg                                                               <  1:121766/251‑1 (MQ=255)
                                                                                                                                                                                                   aaaGGCAGACTGGCGGTTAAGCCTGCCGTCAGATGCCCCAGACAGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTGGGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCAGGGAAATAGCCCGGGACTACTGAGAGCGGCAACCCGATGGTGAGGTAGCTGGc                                                      <  1:172620/249‑1 (MQ=255)
                                                                                                                                                                                                          gACTGGCGTGTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAAGCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGGGAGGTAGCGGGCGAAGTTa                                               <  2:196115/251‑1 (MQ=255)
                                                                                                                                                                                                            cTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACCTCATGGACATAGCCCGGGAATACAGCGAGCGGCAACCCGGTGGCGCGGTAGCTGGCGAAGTTaa                                              >  2:383887/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaTCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACtat                                    <  1:162579/251‑1 (MQ=255)
                                                                                                                                                                                                                         cTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGACATGAGGGCGGCCCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATag                                  >  1:465795/1‑249 (MQ=255)
                                                                                                                                                                                                                                               cGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGGCTATAGAGACAATTCGCAAATTCAGgcgc           >  2:6747/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     aaaaGCCGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAAc                                            >  1:125443/1‑212 (MQ=255)
                                                                                                                                                                                                                                                     aaaaGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAAc                                            <  2:125443/212‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaaGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGt    <  1:230757/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        agcagcCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGttt  <  2:104517/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
CAGTGATAAAACACGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGTTT  >  NC_000913/3611207‑3611704

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: