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.

TGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGC  >  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)
      aaCTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGGGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCga‑g                                                                                                                                                                                                                               <  1:35354/251‑1 (MQ=255)
          gAATAAAGGCTTGCTCCATGGTGGGATCGGGATGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATCTTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAcggcgggaggcggctcggcgggaaa                                                                                                                                                                                                                             >  1:544435/1‑231 (MQ=255)
            aTAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGTCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAACTCACGGCGGGGGAGGGGGGCCACGCCGGagggtggtgcgtgcaaagggg                                                                                                                                                                                                                           >  1:469963/1‑230 (MQ=255)
                aGGCTTGATCCATGGTGGGATCGGGTTGCTACTCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGGGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGACGGAAGGCGGGTTGTCGA‑GAAACAGAATg                                                                                                                                                                                                                     >  1:467061/1‑251 (MQ=255)
                    ttGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTAACACCCGGAGTAGGTTTGTTGA‑GGAAAAGAATGTCCg                                                                                                                                                                                                                 >  1:530555/1‑251 (MQ=255)
                          cATGGTGGGATCGGGTTGCTCATCGGTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCAt                                                                                                                                                                                                          <  2:368122/251‑1 (MQ=255)
                           aTGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTACGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGCGGGGGTCAACGCCGGAAGTCGGGTCGTTGA‑GATACAAAAAGTCCgggtacgg                                                                                                                                                                                                          >  1:433401/1‑245 (MQ=255)
                             ggtggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGc                                                                                                                                                                                                        <  1:432562/250‑1 (MQ=255)
                               tggGTTCGGGTTGCTCATCGTTAGCCGACTGGGCTTTCAAAGCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTGCACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCTA‑GAAACAGAATGTCCGG‑TTCATGCAt                                                                                                                                                                                                      <  2:471329/250‑1 (MQ=255)
                                 ggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAAAGCCCGACGTCGGTTCGGCGA‑GAAACAAAATATCCGG‑tggatcgggagg                                                                                                                                                                                                   >  1:552066/1‑240 (MQ=255)
                                 ggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCAGGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑t                                                                                                                                                                                                              <  2:552066/240‑1 (MQ=255)
                                             caaTCGTTAGCCGCCTGGGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGTTGCGTTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTTACCATCACCGTGACGCCTTTCTCTACCAGGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGgccac                                                                                                                                                                                       <  2:305951/249‑2 (MQ=255)
                                               aTCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTAACTCCCGAAGGCGGTTCGTCGACGAAACA‑AAAGTCCCG‑GTCATGCATCAGCCGACAGgccca                                                                                                                                                                                       >  1:119499/1‑247 (MQ=255)
                                                tCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGA‑GAAACAGAATGTCCGG‑GTCATGCATCAGCGAACAGgccagcgc                                                                                                                                                                                    >  1:562532/1‑251 (MQ=255)
                                                    tAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACATgccagcgccag                                                                                                                                                                                 <  2:344233/250‑1 (MQ=255)
                                                      gCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGAGGGGGGCAACGCCCGGAGGCGGTTCGTCGA‑GGAACAGAATGTCCGG‑GTCCTGCATCCGCGAAAAGgccagcccccgccc                                                                                                                                                                              >  1:506886/1‑243 (MQ=255)
                                                          aCTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGt                                                                                                                                                                          <  1:402526/251‑1 (MQ=255)
                                                            tgtgCTCTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGTCGGGTCGTCGG‑GGAAAAGAATGTCCGG‑GTCATGCATCCGCGAACAGGCCAGCGCCCGCCGCTGttt                                                                                                                                                                        >  1:155411/1‑251 (MQ=255)
                                                            ctggctgtCAAATCGTCCGGCGTGCTGCTGGCTATTAATTTCCCGCGGGACCCCAGGCCGATGCGGGCGCAATATTCCGCTTCATCCATAAAGGGGGTGGTGACCATCCCCGTGATGCCTTTCGCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAATAGGTCAGCGCCAGCCGCTGt                                                                                                                                                                          <  2:504092/242‑1 (MQ=255)
                                                             gtgCTTTCAAATCGACCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGTTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TGGATGCATCAGCTAACAGGCCAGCGCCAGCCGCTGTTTa                                                                                                                                                                       <  1:89764/251‑1 (MQ=255)
                                                              tgCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATCTTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGGAGTCGGGTTGTAGA‑GAAACCGAAAGTCCGG‑TTTATGCATCCGCGAACAGGCCAGCGCCCGCCCCCGTTTa                                                                                                                                                                       >  1:580875/1‑250 (MQ=255)
                                                              ggCTTGCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCTCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTGTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGTGCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATTTCCGG‑TTCATGCAGCAGGGAACAGTCCAGCGCCAGCCGCTGTTTa                                                                                                                                                                       <  2:580875/249‑1 (MQ=255)
                                                                             cGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGtt                                                                                                                                                      <  1:547433/251‑1 (MQ=255)
                                                                                      gCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTgg                                                                                                                                           <  1:117614/251‑1 (MQ=255)
                                                                                            gATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAAAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCCGGTCA‑T‑CGGTGGCGgggg                                                                                                                                     >  1:99652/1‑251 (MQ=255)
                                                                                                           gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAATCGGTTCGTCGA‑GAAACAGAATGTCCGG‑GTCATACATCAGCGCACAGGCCAGCGCCAGCCCCTGGTTAAAAACTAAA‑GGCAGTTCACT‑CGG‑GGCGTGGGGGGCGAGACATTT‑c                                                                                                                      >  1:388669/1‑250 (MQ=255)
                                                                                                                   ggCCGATGCGGTCGCAATATGCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAAGTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAAc                                                                                                             <  2:464956/251‑1 (MQ=255)
                                                                                                                       gATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAGCCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATTCTTTT‑CAGGCCGATCGCCt                                                                                                         <  2:593407/251‑1 (MQ=255)
                                                                                                                                gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCCAAGCCGGAAGTCGGGTCCTCGA‑GAAACAGAAAGTCCGG‑GTCATGCCTCAGCGAAAAAGCCAGCGCCAGCCGGTGTTTAAAAACCAATAGCCAGT‑CA‑TACGG‑GGCCGGGGAGGCGGTACTTTT‑TATGCCGAACGCCCCGCTCaacc                                                                                                >  2:203347/1‑248 (MQ=255)
                                                                                                                                       tCCGCTTCATCCATAAAGTGGGTGGTGACCATCTCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGgaga                                                                                          <  1:337792/250‑1 (MQ=255)
                                                                                                                                        ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGAt                                                                                         <  1:180990/250‑1 (MQ=255)
                                                                                                                                          gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGATGGGGGGGTCAACGCCGGGAGTGGCTTTTTCCG‑GAAACAGAATGTCCGG‑TTCCTGCATCAGCGAAAAGGCCAACGCCAGCCGCCGTTTAAAACCTAAT‑GGCCGGTTA‑T‑TGGGGGCGGGGGAAGGGGGACTTTT‑CAAGCCGCACGCCTCCCCCATGCGGGGGAtttt                                                                                      >  1:573593/1‑251 (MQ=255)
                                                                                                                                              cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCg                                                                                  <  1:10429/251‑1 (MQ=255)
                                                                                                                                               attcATAAATTGGGGGGTGACCAGCACCGTGACGCCTGTCTCTACCATTCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGTGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TGCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGAGACTGGT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGt                                                                                 <  1:46323/248‑1 (MQ=255)
                                                                                                                                               aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGCGGGGGGCAACGCCGGAAAGCGGGGCGGCGA‑GAAACAGAAAGGCCGG‑GTCATGCATCAGCGGACAAGCCAACGCCCCCCCCTGGTTTAAACCCAAT‑GGCAGGTCA‑T‑CGGTGGGGGGGGGGGCGCTAATTTT‑TATGCCGGAAGCCTCGCTCATGCGGGAGGtttttttgt                                                                                 >  1:305466/1‑248 (MQ=255)
                                                                                                                                                           ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GAc                                                                    <  1:304020/251‑1 (MQ=255)
                                                                                                                                                               gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACcgcg                                                                <  2:254228/251‑1 (MQ=255)
                                                                                                                                                                     aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGcc                                                          <  1:538875/251‑1 (MQ=255)
                                                                                                                                                                      tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAATCGGGTCGTCGG‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGGACAGGCCAGCGCCCGCCGCCGGTTAAAACCCAAT‑TGCAGTTTC‑T‑CGGTGGCGTGGGGGGCGGTACTTTT‑CAGGACGAAACCCCCGCTCACGCGGGGGATTTTTTTGTTTTGCGCCC‑GCCCGCGTaaaccc                                                         >  1:156080/1‑247 (MQ=255)
                                                                                                                                                                          cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATaca                                                     <  1:270696/251‑1 (MQ=255)
                                                                                                                                                                             gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATacaca                                                   <  2:562532/250‑1 (MQ=255)
                                                                                                                                                                                  cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGGAGACGGGTCGTCGA‑GAAACAGAAAGGCCGGGTTC‑TGCATCAGCGAAAAAGACAACGCCCGCCGCTGTTTAAAAACTCAA‑GGCGGATTA‑T‑CGGGGGGGGGGGGGGGGAGACTTTTTCAGGC‑GGACGCCTCCCTCCAGCGGGGGGGTTTTTTGTTCTGCGCCC‑CACCCCGTAAGCCCTacaccacag                                               >  1:307037/1‑244 (MQ=255)
                                                                                                                                                                                    ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑gaga                                           <  2:155411/250‑1 (MQ=255)
                                                                                                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GTCAGTTCA‑T‑CGGTGGCg                                                                                                                                         <  2:348212/153‑1 (MQ=255)
                                                                                                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGACGGCGGGGATAACGACGGCCGTCGGTTCGTCAA‑GAAAACGAATTTCCGG‑TTACTACATCCACGAACACGCACCCGACAGCCGCTGGTTAAAAACCAAT‑CGAAGTTTA‑T‑AGGAGGGGGGGGGGGCGATACTTTT‑TAAGCCCAACCCCTCGCTCCTGAGGGGGGGTTTTTTGTTCTGCGCCC‑GACCCCGTAAGCCATACACACCC‑CAGaaaaaacg                                    >  1:588673/1‑248 (MQ=255)
                                                                                                                                                                                             aTGCTGTTGATGTGCAGCCATAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTaa                                 <  2:530555/251‑1 (MQ=255)
                                                                                                                                                                                                 tGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTAGA‑GAAACAGAATGTACGG‑GTTATTCCTCAGCGAACAGGGCCGCGCCAGCCGCTGTTTTAAACCTAAT‑GGGCGTTCA‑C‑CGGGGGCGGGGGGGGGGGGACTTTT‑TAGGCCGGACCCCCCGCCCCTGCGGGGGGGTTTTTTGTTCTGCGCCCCGCCCG‑GTAAACCCTACACACCC‑GAGaaaaaaagtaaaatt                             >  2:31171/1‑241 (MQ=255)
                                                                                                                                                                                                    ggATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCGCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTg                           <  2:359987/249‑1 (MQ=255)
                                                                                                                                                                                                              ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTGTTCGACCGTCa                <  1:21419/251‑1 (MQ=255)
                                                                                                                                                                                                                aaaaTTCACGGGGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGCCGA‑GAACCAGAATGTCCGG‑TTCATGAAGCAGCGAACAGGCCAGCGCCAGCCGCTGGTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGTATATTTTTTCGGTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAATAAGCGTAAATTCTGTTCGACCTTCAgg              <  2:374571/251‑1 (MQ=255)
                                                                                                                                                                                                                      cACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGGTCGTCGA‑GAAACAGAATGCCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGTCCGAACGCCTCGCTCATGCGGGAGATTTTGTCGTTCTGCGCCC‑GACCGCTTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTAcc         <  2:527442/250‑1 (MQ=255)
                                                                                                                                                                                                                            gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGc  >  1:2874/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAGGGGGTCAACGACGGCAGTCGGTTAGTCGA‑GAAACCGAATGTCCGG‑TTCATGCATCCGCGCACAAGCCAGCGCCAGACGCTGTTTAAAACCTAAT‑GGCAGTTAC‑T‑CGGTGGTGTGGGCGGCGAAACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACAC‑CCACGAGAAAAAGCGtcacatatgtt                         >  2:552487/1‑214 (MQ=255)
                                                                                                                                                                                                                                gAGGGGGCCAACGCCGTTAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGGGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCTTAAGCCGAACACACCA‑GAGAAAAAGCGTAAATTCTGtt                         <  1:552487/224‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGAttt                                                                                       <  2:405717/163‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGAttt                                                                                       >  1:405717/1‑163 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGA‑GAAACAGAATGTCCGG‑TTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT‑GGCAGTTCA‑T‑CGGTGGCGTGGGAGGCGATACTTTT‑CAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCC‑GACCGCGTAAGCCATACACACCA‑GAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTAGAGC  >  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: ‑