Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,2350TA32.3% 81.5 / 19.0 64M292L (ATG→TTG) ybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/8);  new base A (0/20);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 3.38e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCC  >  NC_000913/825010‑825448
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGACAAaaac                                                                                                                                                                                                              >  1:335071/1‑238 (MQ=255)
          aaaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGGAAAGTGGCGAATAGGGGTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGACCACAACGACAAACACGCCCATaaaa                                                                                                                                                                                                    >  2:994431/1‑237 (MQ=255)
            aCATCGCGTGCGCCGATGCCCCTGTCGTGGCCGTTATCCCCAGTAGCGGCCACAAACTTTCCCACACAATATCCAGACTCGCTTCCTTCAAATAAATCTGCTGGGTAAGGCCCGAAAAGTGGCAAATAGGGTTAACCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTCCCACCGAAGAAACGTAACCGGAAAGGAGAATGGGGGGCATCAAAAAGACAAACACGCCGATAAACGcc                                                                                                                                                                                                 <  2:1513112/239‑1 (MQ=255)
                cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCccgcc                                                                                                                                                                                             >  1:298078/1‑238 (MQ=255)
                 gcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAAAACGCCGAAAAACACct                                                                                                                                                                                                >  1:1424218/1‑235 (MQ=255)
                   gTACGCCGCTGACCATGCCGGGGCCTTTATCACTGGTGGCGCCAACAAACTATCCCAAACAATACCCAGGCACGCATCCTTCAAAAAATTCTGCTTGTTAAGGTCAGTAAAGTGGCAAATAGGGTTAACCCACGCCAGTTTTTGCAGCCAAACCGGCATTTTTTCCACCGGAGAAACGTAACGGAAAGGGAAAAGGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                          <  2:619837/239‑1 (MQ=255)
                      cgccgcTGCCCCTGTGTTGGCCGTTCTCACCAGTAGCGGCCACAAATTATCCCAAACAATATCCAGACTCGCATCCTTCAAAAAAATCTGTTGGGTATGGCCCGAAAAGTGGCGAAGGGGGTAAATCCACGTCGGGTTTTGCAGCAATACCGGCATGTTTCCCACCGGAGAAACGTAACCGGAAGGGGGAGGGGCGGGCATCATAAAGAAAAACACGCCGATAAACGCCTGCtgttgtg                                                                                                                                                                                       <  2:1178858/239‑1 (MQ=255)
                                   gTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAAAAAACACGCCGATAAAAGCCTGCTTTTGGGTTTGACAAAGtga                                                                                                                                                                          >  1:34693/1‑239 (MQ=255)
                                        ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAAAACCGATAAAAGCCTGCTGTTGTGTTGAAAAGAGTGAGGaaa                                                                                                                                                                     >  1:1244594/1‑239 (MQ=255)
                                                cccccAGAAGCGACCAAAAACTATCCCACAAAATACCCAGACTCGCATCCTAAAAATAAATCTGCTTGGTAATGTCAGAAAAGGGGCAAATAGGGTTAAACCACGTCAGGTTTTGCAGCCATACCGGCATTTTTCCCACCGGAGAAACGAAACCGAAAAGGAAAAGGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCctgctg                                                                                                                                                                                            <  2:764513/205‑1 (MQ=255)
                                                       tAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACCCCGATAAACGCCCGCTGTTGTGGTGAACAGAGTGGTGAAAACAACAGaaccaaa                                                                                                                                                       >  2:312357/1‑232 (MQ=255)
                                                          cGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCCATAAACGCCTGCTGTTTTGTTGaa                                                                                                                                                                                  >  1:764513/1‑208 (MQ=255)
                                                            gCCACAAACTATTCCACACAATATCCAGACTCGCCTCCTTAAAATAAATCTGCCTGGTAATGTCCGTAAAGTGGCGAATAGTGTTAATCCACGTCAGGTTTTGCAGCCATACCCGCCTGTTGTCCACCGGAGAAACGCAACCGGTAAGGCGAATGGCGGGGTTCAGACAAGCCAAAAAGCCCAGAAAACCCCGGCTGTTGGGTGGAACAGCGGATGAAATAAAAAGaaccaaatcaca                                                                                                                                                  >  2:991239/1‑227 (MQ=255)
                                                                        tCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCCGAAAAACGCCTGCTGTTGGGTTGAACAGAGGGATGAAATCAAAAGACCGAATCCCACCAACGATAAACCa                                                                                                                                    >  1:723068/1‑239 (MQ=255)
                                                                                 tatCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAACGATAAACCATAAATCAcc                                                                                                                           >  1:980471/1‑239 (MQ=255)
                                                                                 tatCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACGCCGATAAAAGCCTGCTGGTGGGTTGAAAAAAGGGATGAAAACAACAAACCGAATCCCACCAACGGTAAAACATAAAACAcc                                                                                                                           >  1:295339/1‑239 (MQ=255)
                                                                                         cTCCCATCCTTAAAATAAATCTGTTGGGTAATTTCCGAAAAGTGGCAAATAGGTTTAATCCACGTCAGGTTTTGCAGCCATCCCGGCATTTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCg                                                                                                                       <  2:295339/235‑1 (MQ=255)
                                                                                         atgcaaTCCTTAAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCCCGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTACCCGGAAAGGGGAATGGCGGGCATAATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCg                                                                                                                       <  2:1384347/230‑1 (MQ=255)
                                                                                           cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAAAAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATAACCATCGTAAAGt                                                                                                                 >  1:1603637/1‑239 (MQ=255)
                                                                                           cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCCCGCTGTTGTGTTGAACAGAGTGATGAAATCAAAAGACCGAATCCCACCAGCGATAAACCATAAAACACCATCGGaaaat                                                                                                                 >  1:82084/1‑237 (MQ=255)
                                                                                             cATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGCTAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCCTCGTAAAGTAg                                                                                                               >  1:266504/1‑239 (MQ=255)
                                                                                                 cTTAAAAAAAATCGGCTTGTAAATGCCCGAAAGGGGGCAAAAGGGTTAAATCCCGGTCAGGTTTTGCAGCCATACCGGCATTTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAAt                                                                                                                                                       <  2:1433664/196‑1 (MQ=255)
                                                                                                     aaataaatCTGTTGGTTAATGCCCGAAAAGTGGCAAAGAGGTTAAACCCACGTAGGTTTTGGCAGCCATACCGCATTGTTTTCCACCGGAAAAACGTACCGGGAAAGGAAAGGGGGGGGCATCAAAAAAACAAACACGCCGATAAACGCctg                                                                                                                                                                                               <  2:683735/152‑1 (MQ=255)
                                                                                                       ataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAAACAAACACGCCCATAAAAGCCCGCtgttgtg                                                                                                                                                                                       >  1:748194/1‑158 (MQ=255)
                                                                                                           atCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACAAGCCGATAAAAGCCCGCTGGTGTGTTGAACAGAGGGATTAAATCAACAGACCGAATCCCACCAGCg                                                                                                                                             >  1:1433664/1‑196 (MQ=255)
                                                                                                            tCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGGATCATAAAGACAAAAACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATACCACCAGCGATAAAACATAAAACACCAACGTAAAATAGAAAAACACCCGCGa                                                                                                >  2:379544/1‑238 (MQ=255)
                                                                                                               gCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCCTGCtgttgtgtt                                                                                                                                                                                     >  1:683735/1‑152 (MQ=255)
                                                                                                                cTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATaaa                                                                                                                                       >  1:474189/1‑196 (MQ=255)
                                                                                                                         gTCCGTAAAGTGGCAAATAGGGTTAATCCACGTCGGGTTTTGCAGCCATACCGGCATGTTTGCCACCGGAGAAACGTAACCGGAAAGGAAAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCCCCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGa                                                                                  <  2:337994/239‑1 (MQ=255)
                                                                                                                          tCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAAACATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGAt                                                                                 >  1:1500213/1‑239 (MQ=255)
                                                                                                                            cgaaaaGTGGCAAATAGGGTAAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAACACGCCGATAAACGCctg                                                                                                                                                                                               <  1:31675/126‑1 (MQ=255)
                                                                                                                                agtgggCGAAGAGGTTAAATCCAGTTCAGTTGTTGCACCCATCCCGGCATGTTTTCCACGGGAGAACCGTAACCGGAAGGGAGAATGGCGGGCATCAAAAAAACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATa                                                                           <  2:1244594/235‑1 (MQ=255)
                                                                                                                                   tGGCAAATAGGGTTAATCCCCGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCccac                                                                                                                                                  <  1:668436/167‑1 (MQ=255)
                                                                                                                                     gagAATAGGGTAATACCTCGTAGGGTTTGGCAGCCAAACCGGCATTTTTTCCACCGGAGAAACGAAACGGGAAAGGAGAATGGGGGGCATCATAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGcc                                                                       <  2:1424218/236‑1 (MQ=255)
                                                                                                                                      cGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAAAACCCCGATAAACGCCTGCtgttgtgtt                                                                                                                                                                                     >  2:31675/1‑129 (MQ=255)
                                                                                                                                             ggTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAAACAACAGACCGAATCCCCCCAGCGATaaa                                                                                                                                       >  2:668436/1‑167 (MQ=255)
                                                                                                                                                   acccccGTCGGGTTTTGCAGCCATCCCGGCATGTTTCCCACCGGAAAAACGTAACCGGAAAGGAAAAGGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGc                                                        <  2:249888/234‑1 (MQ=255)
                                                                                                                                                   aTCCACTGCAGTTTTAGCAGCGATGCCGACATGTTTTCCACCGGAAAAACGTAACCGAAAGGGAAAATGGCGGCCATCAAAAACACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGc                                                        <  2:324807/239‑1 (MQ=255)
                                                                                                                                                   aTCCACGTCAGGTTTTGCTGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGAAAAACACGCCGAAAAACGCCCGGTGTTTTGTTGAACAGAGGGGTGAAATCAACACACCGAAACCCACCAGCGCTAAAACATAAAACAACATCGCAAAGGAGAATAGCGCCAGCGATCCGGCGAAGGGGAATTGATACCACCCGATACCAAACg                                                         >  2:1555472/1‑238 (MQ=255)
                                                                                                                                                      cACGTCAGGTTTTGCAGCCATCCCGGCATTTTTTCCACCGGAAAAACGTAACCGGAAAGGAGAATGGCGGGCATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAg                                                     <  1:312357/239‑1 (MQ=255)
                                                                                                                                                        cGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATAAGCCCAGATACCCATCGCCCGc                                                    >  1:249888/1‑238 (MQ=255)
                                                                                                                                                               ttgtGCAGCAATCCCGGCATTTTTCCCACCGGAGAAACGTAACCGAAAAGGAGAATGGCGGGCATCAAAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAAtgg                                             <  2:298078/235‑1 (MQ=255)
                                                                                                                                                                  tGCAGCCATACCGGCATGTTTCCCACCGGAAAAACGTACCCGGAAAGGAGAATGGCGGGAATCAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACag                                                                                                                                                                               <  1:1616016/108‑1 (MQ=255)
                                                                                                                                                                        ccATCCCGGCATGTTTCCCACGGGAAAAACGTACCCGAAAAGGAGAATGGCGGGAATAAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGa                                    <  1:1437289/239‑1 (MQ=255)
                                                                                                                                                                           tcccGGCATGTTTCCCACCGGAAAAACGTAACCGGAAAGGAGAAGGGCGGGCATCAAAAAAAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGtg                                                                                                                                                                           <  2:1024705/102‑1 (MQ=255)
                                                                                                                                                                            cccGGCATGTTTCCCACGGGAAAAACGTAACCGGAAAGGAGAAGGGGGGGCATAAAAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCccacca                                                                                                                                                <  2:1466855/129‑1 (MQ=255)
                                                                                                                                                                             ccGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGGTGTGTTGAACAGAGTGATGaaa                                                                                                                                                                     >  2:1616016/1‑108 (MQ=255)
                                                                                                                                                                                   agTTTTCCACGCGAAAAAACGGAACCGAAAAGGAGAAGGGGGGGAATCAAAAAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGa                           <  2:266504/237‑1 (MQ=255)
                                                                                                                                                                                     ttttCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCaa                                                                                                                                                                 >  1:1024705/1‑104 (MQ=255)
                                                                                                                                                                                      tttCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATTAAGACAAACACCCCGATTAACGCCTTTCGTTGTGTTGTAAAGGGAGATGAAAAAAAAAGAACGAATACCACCAGCGATAAAACCATAATCAACCCCGTAAAATAAACCACCGCACAGGGACCGGCGAAGGGGGTTTGTATCTCCCAGGATACCATAGCCCGCACAAAAGAGGCCTTGAAGGTggggctcatt                      >  2:880797/1‑231 (MQ=255)
                                                                                                                                                                                      tttCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTTAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAAcg                                                                                                                                     >  1:1466855/1‑129 (MQ=255)
                                                                                                                                                                                           ccaccgaaaaaCGTTAACGGGAAAGGGAAATGGAGGGCATCAAAAAGAAAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAAc                   <  2:723068/231‑1 (MQ=255)
                                                                                                                                                                                                gagaAACGTAACCGGAAAGGAGAATGGCGGGCCTCATAAAGACAAAAACCCCGAGAAACACCTGCCGTTGTGTTGAACAGAGTGATGAAATCAACAGACCCAATCCCACCAGGCATAAACCATAAAACACCAAACGAAAGGAGAACAAGGCCAAGGATTCCGCGAAGGGGATTTGATACCCCCCGAAACCAAACGACAGCACACTAGGGGCCTGGAAGGGGGCGACAAATAAAACCGGt             >  2:206566/1‑239 (MQ=255)
                                                                                                                                                                                                     aCGTAACCGGAAAGGAGAATGGCGGGCATCCTAAAGGCAAACACCCCGATAAACGCCCGCTGTTGTGTTTAACAAAGTGATTAAATCAACCGACCGGATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAAACGCCCGCCCCATGGTGGCCTGGAGAGTGGCGACAATTAAACCCGGGAcc          >  2:793193/1‑236 (MQ=255)
                                                                                                                                                                                                      cGTACCCGAAAAGGAGAAGGGCGGGCTTAAAAAAGACAAACACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                          <  2:1149615/63‑1 (MQ=255)
                                                                                                                                                                                                          aaCCGGAAAGGAAAATGGCGGGCATCAAAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCttt     <  1:691290/238‑1 (MQ=255)
                                                                                                                                                                                                              ggAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAAAAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACCGCGCCCGCGAACCGGCGAAGGGGGATTGAAAAGCCCCGATACCAAACGCCCGCACAATGGTGGCCTGGAAGGTGGGGACAATTAAAACCGGTCCGGCCTTGcc  >  2:928279/1‑237 (MQ=255)
                                                                                                                                                                                                               gAAAGGAGAATGGCGGGCATCATAAAGCCAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCa                                                             <  1:1264219/177‑1 (MQ=255)
                                                                                                                                                                                                                 aaGGAGAATGGCGGGCATCATAAAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACa                                                                                                                                                                                >  1:1149615/1‑63 (MQ=255)
                                                                                                                                                                                                                         tGGCGGGCATCATAAAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAAAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAAAAGCGCCAGCCAACCGGCGAAAGGGAATTGATACCCCCCGAAAACAAACGCCCGca                                                   >  2:1264219/1‑177 (MQ=255)
                                                                                                                                                                                                                          ggCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                  >  2:1290816/1‑52 (MQ=255)
                                                                                                                                                                                                                          ggCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                  <  1:1290816/52‑1 (MQ=255)
                                                                                                                                                                                                                              ggCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCg                                                                                                                       <  2:1469851/106‑1 (MQ=255)
                                                                                                                                                                                                                              ggCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCg                                                                                                                       >  1:1469851/1‑106 (MQ=255)
                                                                                                                                                                                                                                     tAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGc                                                        >  1:32332/1‑160 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
CCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAAAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGATAAACCATAAATCACCATCGTAAAGTAGAACAGCGCCAGCGATCCGGCGAAGGGGATTTGATACGCCCAGATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCC  >  NC_000913/825010‑825448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: