Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG29.3% 81.4 / 16.6 60D492A (GAC→GCC) ybhFABC exporter ATP binding subunit YbhF
Reads supporting (aligned to +/- strand):  ref base T (7/34);  new base G (16/1);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 4.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGC  >  NC_000913/827276‑827746
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
gggATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTGAGCCGCCTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATGTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGGGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                                                                                                                                                                                      <  2:441819/250‑1 (MQ=255)
    tccgttctttAAATGCTTGCTCCATGGTGGAAGCGGGTTGCTCATCGTTACCCGACTGTTCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCGGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTGCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGCGGTCAACGCCGGAAGTCGGTtcgtc                                                                                                                                                                                                                                   <  2:469963/241‑1 (MQ=255)
    tCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGGAACGCCCGAACACGGTtcgtcg                                                                                                                                                                                                                                  >  2:402526/1‑251 (MQ=255)
     cAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGTTCATCGTTAGCCGACTGTGTTTGCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGGGGCGGCTGTCTCGACCATGCTGTTGATGTGCAGTCAAAATTCACGGCGGGTGAGGGGGTCAGCGCCGGAAGTTGGGtcgtcg                                                                                                                                                                                                                                  <  2:567909/250‑1 (MQ=255)
      aaCTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGGGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCgag                                                                                                                                                                                                                               <  1:35354/251‑1 (MQ=255)
          gAATAAAGGCTTGCTCCATGGTGGGATCGGGATGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATCTTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAcggcgggaggcggctcggcgggaaa                                                                                                                                                                                                                             >  1:544435/1‑231 (MQ=255)
            aTAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGTCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAACTCACGGCGGGGGAGGGGGGCCACGCCGGagggtggtgcgtgcaaagggg                                                                                                                                                                                                                           >  1:469963/1‑230 (MQ=255)
                aGGCTTGATCCATGGTGGGATCGGGTTGCTACTCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGGGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGACGGAAGGCGGGTTGTCGAGAAACAGAATg                                                                                                                                                                                                                     >  1:467061/1‑251 (MQ=255)
                    ttGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTAACACCCGGAGTAGGTTTGTTGAGGAAAAGAATGTCCg                                                                                                                                                                                                                 >  1:530555/1‑251 (MQ=255)
                          cATGGTGGGATCGGGTTGCTCATCGGTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                                                                                                                                                                                          <  2:368122/251‑1 (MQ=255)
                           aTGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTACGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGCGGGGGTCAACGCCGGAAGTCGGGTCGTTGAGATACAAAAAGTCCgggtacgg                                                                                                                                                                                                          >  1:433401/1‑245 (MQ=255)
                             ggtggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                                                                                                                                                                                        <  1:432562/250‑1 (MQ=255)
                               tggGTTCGGGTTGCTCATCGTTAGCCGACTGGGCTTTCAAAGCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTGCACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCTAGAAACAGAATGTCCGGTTCATGCAt                                                                                                                                                                                                      <  2:471329/250‑1 (MQ=255)
                                 ggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAAAGCCCGACGTCGGTTCGGCGAGAAACAAAATATCCGGtggatcgggagg                                                                                                                                                                                                   >  1:552066/1‑240 (MQ=255)
                                 ggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCAGGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGt                                                                                                                                                                                                              <  2:552066/240‑1 (MQ=255)
                                             caaTCGTTAGCCGCCTGGGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGTTGCGTTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTTACCATCACCGTGACGCCTTTCTCTACCAGGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccac                                                                                                                                                                                       <  2:305951/249‑2 (MQ=255)
                                               aTCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTAACTCCCGAAGGCGGTTCGTCGACGAAACAAAAGTCCCGGTCATGCATCAGCCGACAGgccca                                                                                                                                                                                       >  1:119499/1‑247 (MQ=255)
                                                tCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGgccagcgc                                                                                                                                                                                    >  1:562532/1‑251 (MQ=255)
                                                    tAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACATgccagcgccag                                                                                                                                                                                 <  2:344233/250‑1 (MQ=255)
                                                      gCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGAGGGGGGCAACGCCCGGAGGCGGTTCGTCGAGGAACAGAATGTCCGGGTCCTGCATCCGCGAAAAGgccagcccccgccc                                                                                                                                                                              >  1:506886/1‑243 (MQ=255)
                                                          aCTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGt                                                                                                                                                                          <  1:402526/251‑1 (MQ=255)
                                                            tgtgCTCTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGTCGGGTCGTCGGGGAAAAGAATGTCCGGGTCATGCATCCGCGAACAGGCCAGCGCCCGCCGCTGttt                                                                                                                                                                        >  1:155411/1‑251 (MQ=255)
                                                            ctggctgtCAAATCGTCCGGCGTGCTGCTGGCTATTAATTTCCCGCGGGACCCCAGGCCGATGCGGGCGCAATATTCCGCTTCATCCATAAAGGGGGTGGTGACCATCCCCGTGATGCCTTTCGCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAATAGGTCAGCGCCAGCCGCTGt                                                                                                                                                                          <  2:504092/242‑1 (MQ=255)
                                                             gtgCTTTCAAATCGACCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGTTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTGGATGCATCAGCTAACAGGCCAGCGCCAGCCGCTGTTTa                                                                                                                                                                       <  1:89764/251‑1 (MQ=255)
                                                              tgCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATCTTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGGAGTCGGGTTGTAGAGAAACCGAAAGTCCGGTTTATGCATCCGCGAACAGGCCAGCGCCCGCCCCCGTTTa                                                                                                                                                                       >  1:580875/1‑250 (MQ=255)
                                                              ggCTTGCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCTCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTGTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGTGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATTTCCGGTTCATGCAGCAGGGAACAGTCCAGCGCCAGCCGCTGTTTa                                                                                                                                                                       <  2:580875/249‑1 (MQ=255)
                                                                             cGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGtt                                                                                                                                                      <  1:547433/251‑1 (MQ=255)
                                                                                      gCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTgg                                                                                                                                           <  1:117614/251‑1 (MQ=255)
                                                                                            gATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAAAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCCGGTCATCGGTGGCGgggg                                                                                                                                     >  1:99652/1‑251 (MQ=255)
                                                                                                           gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAATCGGTTCGTCGAGAAACAGAATGTCCGGGTCATACATCAGCGCACAGGCCAGCGCCAGCCCCTGGTTAAAAACTAAAGGCAGTTCACTCGGGGCGTGGGGGGCGAGACATTTc                                                                                                                      >  1:388669/1‑250 (MQ=255)
                                                                                                                   ggCCGATGCGGTCGCAATATGCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAAGTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAAc                                                                                                             <  2:464956/251‑1 (MQ=255)
                                                                                                                       gATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAGCCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATTCTTTTCAGGCCGATCGCCt                                                                                                         <  2:593407/251‑1 (MQ=255)
                                                                                                                                gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCCAAGCCGGAAGTCGGGTCCTCGAGAAACAGAAAGTCCGGGTCATGCCTCAGCGAAAAAGCCAGCGCCAGCCGGTGTTTAAAAACCAATAGCCAGTCATACGGGGCCGGGGAGGCGGTACTTTTTATGCCGAACGCCCCGCTCaacc                                                                                                >  2:203347/1‑248 (MQ=255)
                                                                                                                                       tCCGCTTCATCCATAAAGTGGGTGGTGACCATCTCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGgaga                                                                                          <  1:337792/250‑1 (MQ=255)
                                                                                                                                        ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGAt                                                                                         <  1:180990/250‑1 (MQ=255)
                                                                                                                                          gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGATGGGGGGGTCAACGCCGGGAGTGGCTTTTTCCGGAAACAGAATGTCCGGTTCCTGCATCAGCGAAAAGGCCAACGCCAGCCGCCGTTTAAAACCTAATGGCCGGTTATTGGGGGCGGGGGAAGGGGGACTTTTCAAGCCGCACGCCTCCCCCATGCGGGGGAtttt                                                                                      >  1:573593/1‑251 (MQ=255)
                                                                                                                                              cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCg                                                                                  <  1:10429/251‑1 (MQ=255)
                                                                                                                                               attcATAAATTGGGGGGTGACCAGCACCGTGACGCCTGTCTCTACCATTCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGTGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTGCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAGACTGGTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGt                                                                                 <  1:46323/248‑1 (MQ=255)
                                                                                                                                               aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGCGGGGGGCAACGCCGGAAAGCGGGGCGGCGAGAAACAGAAAGGCCGGGTCATGCATCAGCGGACAAGCCAACGCCCCCCCCTGGTTTAAACCCAATGGCAGGTCATCGGTGGGGGGGGGGGCGCTAATTTTTATGCCGGAAGCCTCGCTCATGCGGGAGGtttttttgt                                                                                 >  1:305466/1‑248 (MQ=255)
                                                                                                                                                           ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGAc                                                                    <  1:304020/251‑1 (MQ=255)
                                                                                                                                                               gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACcgcg                                                                <  2:254228/251‑1 (MQ=255)
                                                                                                                                                                     aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGcc                                                          <  1:538875/251‑1 (MQ=255)
                                                                                                                                                                      tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAATCGGGTCGTCGGGAAACAGAATGTCCGGTTCATGCATCAGCGGACAGGCCAGCGCCCGCCGCCGGTTAAAACCCAATTGCAGTTTCTCGGTGGCGTGGGGGGCGGTACTTTTCAGGACGAAACCCCCGCTCACGCGGGGGATTTTTTTGTTTTGCGCCCGCCCGCGTaaaccc                                                         >  1:156080/1‑247 (MQ=255)
                                                                                                                                                                          cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATaca                                                     <  1:270696/251‑1 (MQ=255)
                                                                                                                                                                             gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATacaca                                                   <  2:562532/250‑1 (MQ=255)
                                                                                                                                                                                  cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGGAGACGGGTCGTCGAGAAACAGAAAGGCCGGGTTCTGCATCAGCGAAAAAGACAACGCCCGCCGCTGTTTAAAAACTCAAGGCGGATTATCGGGGGGGGGGGGGGGGAGACTTTTTCAGGCGGACGCCTCCCTCCAGCGGGGGGGTTTTTTGTTCTGCGCCCCACCCCGTAAGCCCTacaccacag                                               >  1:307037/1‑244 (MQ=255)
                                                                                                                                                                                    ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAgaga                                           <  2:155411/250‑1 (MQ=255)
                                                                                                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGTCAGTTCATCGGTGGCg                                                                                                                                         <  2:348212/153‑1 (MQ=255)
                                                                                                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGACGGCGGGGATAACGACGGCCGTCGGTTCGTCAAGAAAACGAATTTCCGGTTACTACATCCACGAACACGCACCCGACAGCCGCTGGTTAAAAACCAATCGAAGTTTATAGGAGGGGGGGGGGGCGATACTTTTTAAGCCCAACCCCTCGCTCCTGAGGGGGGGTTTTTTGTTCTGCGCCCGACCCCGTAAGCCATACACACCCCAGaaaaaacg                                    >  1:588673/1‑248 (MQ=255)
                                                                                                                                                                                             aTGCTGTTGATGTGCAGCCATAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTaa                                 <  2:530555/251‑1 (MQ=255)
                                                                                                                                                                                                 tGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTAGAGAAACAGAATGTACGGGTTATTCCTCAGCGAACAGGGCCGCGCCAGCCGCTGTTTTAAACCTAATGGGCGTTCACCGGGGGCGGGGGGGGGGGGACTTTTTAGGCCGGACCCCCCGCCCCTGCGGGGGGGTTTTTTGTTCTGCGCCCCGCCCGGTAAACCCTACACACCCGAGaaaaaaagtaaaatt                             >  2:31171/1‑241 (MQ=255)
                                                                                                                                                                                                    ggATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCGCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTg                           <  2:359987/249‑1 (MQ=255)
                                                                                                                                                                                                              ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCa                <  1:21419/251‑1 (MQ=255)
                                                                                                                                                                                                                aaaaTTCACGGGGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGCCGAGAACCAGAATGTCCGGTTCATGAAGCAGCGAACAGGCCAGCGCCAGCCGCTGGTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGTATATTTTTTCGGTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAATAAGCGTAAATTCTGTTCGACCTTCAgg              <  2:374571/251‑1 (MQ=255)
                                                                                                                                                                                                                      cACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGCCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGTCCGAACGCCTCGCTCATGCGGGAGATTTTGTCGTTCTGCGCCCGACCGCTTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTAcc         <  2:527442/250‑1 (MQ=255)
                                                                                                                                                                                                                            gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGc  >  1:2874/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAGGGGGTCAACGACGGCAGTCGGTTAGTCGAGAAACCGAATGTCCGGTTCATGCATCCGCGCACAAGCCAGCGCCAGACGCTGTTTAAAACCTAATGGCAGTTACTCGGTGGTGTGGGCGGCGAAACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACCCACGAGAAAAAGCGtcacatatgtt                         >  2:552487/1‑214 (MQ=255)
                                                                                                                                                                                                                                gAGGGGGCCAACGCCGTTAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGGGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCTTAAGCCGAACACACCAGAGAAAAAGCGTAAATTCTGtt                         <  1:552487/224‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGAttt                                                                                       <  2:405717/163‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGAttt                                                                                       >  1:405717/1‑163 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGC  >  NC_000913/827276‑827746

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