Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913782,1940TG29.0% 90.8 / 10.4 62L37R (CTA→CGA) nadAquinolinate synthase
Reads supporting (aligned to +/- strand):  ref base T (18/26);  new base G (0/18);  total (18/44)
Fisher's exact test for biased strand distribution p-value = 6.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCATCCCGCTTCTACTTT  >  NC_000913/781895‑782332
                                                                                                                                                                                                                                                                                                           |                                                                                                                                          
aCTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCt                                                                                                                                            >  2:421712/1‑300 (MQ=255)
   gtgCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGANTATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTaaa                                                                                                                                         >  2:525188/1‑300 (MQ=255)
   gtgCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGANTATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTATGCTaaa                                                                                                                                         >  2:525182/1‑300 (MQ=255)
   gtgCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTaaa                                                                                                                                         >  1:414458/1‑300 (MQ=255)
   gtgCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTaaa                                                                                                                                         >  2:424899/1‑300 (MQ=255)
            aCTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTaa                                                                                                                                >  1:635100/1‑300 (MQ=255)
                tATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGc                                                                                                                             <  2:695832/299‑1 (MQ=255)
                      ggTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATg                                                                                                                      <  1:424899/300‑1 (MQ=255)
                        tGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGt                                                                                                                    <  1:726963/300‑1 (MQ=255)
                          gAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGtt                                                                                                                   >  2:647720/1‑299 (MQ=255)
                           aaTTAGGCTTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTGAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCGAATGTCGGTTGTAACGTTAAACATGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGAGTTATCCTTTCCCCCCGAAGCCGACGCCGGTAAGCATTGATGAAAAAGCGTAGTACCGCGAGGAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATgggggc                                                                                                                 <  1:360258/300‑5 (MQ=255)
                           aaTTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGGAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGc                                                                                                                 <  1:360268/300‑1 (MQ=255)
                                   gTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACTAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCactac                                                                                                          <  2:709250/299‑1 (MQ=255)
                                   gTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCactac                                                                                                          <  1:60798/299‑1 (MQ=255)
                                   gTAAAACGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTCCCCCACCATTCAACTATCTTGTTTAGCATATAAAACAAACTACACCGACAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCAATTTATCCTTTCCCCCCGAAGCTGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCTAGAAGATAAAACGTCTGCTAAAAGAACGTAATCCGGTGATGGTTGCCCactac                                                                                                          <  2:484622/299‑1 (MQ=255)
                                    tAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACTAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCTAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAGAACGTCTGCTAAAAGAACGTAATGCGTTGATGGTTGCCCactat                                                                                                          <  2:341662/298‑2 (MQ=255)
                                     aaaaTGACGCATCCTGCACATTAGGCGTAATTCGAGTTACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTANGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCactac                                                                                                          <  2:218787/297‑1 (MQ=255)
                                       aagaaCGCACCCTGCACATGAGGCGGAATTTGAGTGACTTTTCCCCACCATTCGACTATCTTTTTTAGCATATAAAACAAATTACACCGAGAACAGCGAACATTACGCTAAGGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGTGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGTTGATGGTTGCCCACTACTATAc                                                                                                     <  2:84331/296‑1 (MQ=255)
                                       aaTGATGCATCCTGCACATGATGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCGAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  1:269595/300‑1 (MQ=255)
                                       aaTGACGCATCCTGCACATTAGGCGTAATTCGATTGACTTTTCCCCACCATTCTACTACCTTGTTGAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCGAATGTCGGTTTGAACGTTAAGCCTGTAAAACGAGATGGGAAGATGAGTGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGAGGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  2:508279/300‑1 (MQ=255)
                                       aaTGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTTTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  1:701240/300‑1 (MQ=255)
                                       aaTGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGGAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  2:610552/300‑1 (MQ=255)
                                       aaTGACGCATCCTGCACATGAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATGACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGGAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGAGGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  1:694395/300‑1 (MQ=255)
                                       aaTGACGCATCCTGCACATGAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATAGAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGGAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  1:647720/300‑1 (MQ=255)
                                       aaTGACGCATCCGGCACATGAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCAAATATTACGCTAATGTCGGTTTTAACGTGAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  1:701274/300‑1 (MQ=255)
                                       aaTGACGCACCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATAGAAAACAAATTACACCGACAACAGCGAATATGACACGAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGGAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAc                                                                                                     <  2:421898/300‑1 (MQ=255)
                                          gACGCATCCTGCACATTATGCGTAGTTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACTCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGTGGAATGTTTGATCCAGACACGGCGATTGATCCTTTCCCCCCGAAGCCGACGCCGTGAAGCATTGATGAAAAAGCGTATTACTGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGa                                                                                                  <  1:543643/300‑1 (MQ=255)
                                               aTCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAAAGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAAAGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATAcc                                                                                                    >  2:94072/1‑293 (MQ=255)
                                                       caTTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTcaa                                                                                     >  1:415745/1‑300 (MQ=255)
                                                         ttAGGCGTAGTTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTGAGCATAGAAAACAAATTACACCGATAACAGCGAACATTACTCGAATGTCGGTTTTAACGTTAAACCTGTAAAACGAGATGGTAAGAGGAGCGTAAGGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCAACGCCGTGAAGCATTGAGGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGGGAGGGTTGCCCACTACTATAc                                                                                                     <  2:189697/282‑1 (MQ=255)
                                                            ggCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCaacaac                                                                                 >  2:17056/1‑299 (MQ=255)
                                                                         ttatttttCCCCCCCCTTTGACTACCTTTTTTAGCGTGGAAAACAAATTACACCGGTAACAGCGAAGATGACACTAATGTAGTTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTGGATCCAGACACGGCGGTTTATCCTTTCCCCCCGAAGCCGACGTCGTGAAGCATTGAGGAAAAAGCGTATTACCGTGAGAAGGTAAAACGTCTGCTAAAAGAACGTAATGTGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAAcc                                                                   <  2:418058/296‑1 (MQ=255)
                                                                         tGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAACATTACGCTAATGTCGGTTTTAACGTTAAGCCTGGAAAACGAGATGGGAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTGATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGTGGTGGTGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAAcc                                                                   <  1:491862/300‑1 (MQ=255)
                                                                         tGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATAGAAAACAAATTACACCGATAACATCGAACATTACGCTAGTGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGGAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATGGAGGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGTTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAAcc                                                                   <  1:701900/300‑1 (MQ=255)
                                                                           atttttCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAAGATTACACCAATGTCGTTTTTAACGTGAAGCCTGTAAAACGAGGTGGTAAGATGAGCGGAATGTTGGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTGAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCgg                                                                 <  1:296399/298‑1 (MQ=255)
                                                                              tttCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCACCAACTGGCAGAAGAAACCGGt                                                                >  1:310216/1‑298 (MQ=255)
                                                                                     accgttCGACTACCTTGTTTATCATATAAACAAAATTACACCAAGAACAGCAAATATTACTCTAAGGTCGGTTTTAACTTGAAGCCTGGAAAATGAGATGGTAAGATAGGCGTAATGTTTGATCCAGACATGGCGGTTTGTCCTTTCCCCCAGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGTGAAAAGAACGTAATGCGGGGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTAtt                                                       <  2:603554/296‑1 (MQ=255)
                                                                                     accaTTCGACTACTTTGTTTAGCATAGAAAAAAAATGACACCAAGAACAGCAAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTGAAGCATTGGTGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTAtt                                                       <  2:603659/300‑1 (MQ=255)
                                                                                          tCGACTATCTTGTTTAGCATAGAAAACAAATTACACCGAGAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGGAAAACGAGATGGTAAGAGGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGGTGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTAtt                                                       <  2:603533/295‑1 (MQ=255)
                                                                                                aTCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTctct                                            >  1:718855/1‑300 (MQ=255)
                                                                                                     gTTTAGCATATAAAACAAATTACACCGATAACAGTGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACAGCGATTTATCCTTTCCCCAAGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCTGTGGCTGTCTTTCTGATTCTCTGGaaa                                       >  2:153630/1‑300 (MQ=255)
                                                                                                       ttAGCGTAGAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGCGAAAAGAACGTAATGCGGTGAGGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAAcc                                                                   <  2:418075/270‑1 (MQ=255)
                                                                                                          gCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGgc                                   >  1:593253/1‑299 (MQ=255)
                                                                                                            atatAAAGCCAATTACACCGATAACAGCGCATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTACAACGAGATGGTAAGCTGGGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGGACGTAATGCGGTGCTGGTTGACCACTCCTATACCGATCCCAAAATTCAACAACTGACAGAAGAAACCGGTGACTGTATTTCTGATTctct                                            >  2:36859/1‑288 (MQ=255)
                                                                                                            atatAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGgcgc                                 >  2:36948/1‑299 (MQ=255)
                                                                                                               taAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGGAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGGGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCt                              <  1:490202/299‑1 (MQ=255)
                                                                                                                aaaaCAAATTACACCGATAACAGCGAAGATTACGCTAATGTCGGTTTTAACGTTAAGCCTGGAAAACGAGATGGTAAGATGAGCGGAATGTTTGATCCAGACATGGCGGTTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGCTAAAAGAACGTAATGCGTTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTAtt                                                       <  2:603656/273‑1 (MQ=255)
                                                                                                                 aaaCAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTc                            >  2:683937/1‑299 (MQ=255)
                                                                                                                     aaggtaaatagATAACAGCGAGAATTCCTCTAATGTCTTTTCTAATTTTAAGCCGGGAAAAAGAGATGGTAAGATGAGCGGAATGTTTGATCCAGACGTGGCGATTTATCCTTTCCCCCCGAATCCGAAGCTTTTGAGCGTGGAGGAAAAAGCGGAGTACTGTGAGAAAAGAAACCGTTGGCGAAAAAAATGGAAGGTGGGGAGGGTGGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATTGCACGCTTCggggc                       <  2:283420/290‑4 (MQ=255)
                                                                                                                     aaaTTACACCGATAACATCGAACATTACGCTAATGTCGTTTTTAACGTTAATCCTGGAAAACGAGGTGGTAAGATGAGCTTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACTCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGGGAAAACGTCTGCTAAAAGAACGTAATGTGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCATAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGc                       <  1:83825/300‑1 (MQ=255)
                                                                                                                     aaaTTACACCGATAACANCGAATATTACGCTAATGTNGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTGAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGCGAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCAGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAANTGGCGCGCTTCGGTGc                       <  1:614654/300‑1 (MQ=255)
                                                                                                                          acacCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAACCTGGAAAACGAGAGGGTAAGATGAGCGTAGGGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGAGAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGc                  <  1:213904/300‑1 (MQ=255)
                                                                                                                            acCGATAACATCGAATATTACGCTAATGTCGGTTTTAACGTGAAGCCTGTAAAACGAGATGGGAAGATGAGCGTAAGGTTGGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCTTTAAGCATTGAGGAAAAAGCGTATTACCGCGAGAAGATAAAACGCCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCac                <  1:300713/300‑2 (MQ=255)
                                                                                                                                   acaacGAATATTACGCTAATGTCGTTTTTAACGTTAAGCCTGTAAAACGAGATGGGAAGATGGGTGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGTAATGCGGTGAGGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaa         <  1:298661/296‑9 (MQ=255)
                                                                                                                                      gCGAATATTACTCTAATGTCGGTTTTAACGTGAAGCCTGTAAAACGAGATGGTAAGAGGAGCGGAATGTTGGATCCAGACACGGCGATTTACCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAATGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaatcc      <  1:382683/300‑12 (MQ=255)
                                                                                                                                       cGAATATTACGCGAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGAGGAGTGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCTGTTAAGCATTGATGAAAAAGTGTATTACCGCGAGAAGAGAAAACGTCTGCGAAAAGAACGTAATGCGGTGGTGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaa         <  1:650162/296‑9 (MQ=255)
                                                                                                                                          atatTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaatccccta  >  1:269650/1‑285 (MQ=255)
                                                                                                                                          atatTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTACCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaatccccta  >  2:166438/1‑285 (MQ=255)
                                                                                                                                                             tttAACGTTAAGCCTGTAAAACGAGAGGGTAAGATGAGCGGAATGTGTGATCCAGACACGGCGATTTATCCTTCCCCCCCGAAGCCGACGCCTTTAAGCATTGAGGAAAAAGCGGATTACCGCGAGAAGATAAAACGTCTGCGAAAAGAACGAAATTCGGTGATGGTTTCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGc                  <  2:414458/265‑1 (MQ=255)
                                                                                                                                                                              aaaccGAGATGGTAAGAGGAGCGGAAGGTTTGATCCAGACAAGTCGATTTATCCTTTCCCCCCGAAGCCGACTCCTTTAGGCGTTGGTGAAAAAGCGGATTACCGCGGGAAGATAAAACGTCGGCGAAAAGAACGTAATGCGGTGGGGGTTGCCCACTACTATACCGATCCCGAAATTcaacaa                                                                                  <  1:602272/180‑1 (MQ=255)
                                                                                                                                                                                                        tGTTTGATCCAGAAATGGCGATTTATCCTTTCCCCCCGAAGCTGACGCCGTTAAGCATTGGGGAAAAAGCGGATTACCGCGAGAAGAGAAAACGTCGGCGAAAAGAACGCAATGCGGGGAGGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatg           <  2:47591/229‑7 (MQ=255)
                                                                                                                                                                                                             ggtCCAGACACGGCGGTTTATCTTTTCCCCCCTAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAGGATAAAACGTCTGCTAAAGGAACGTAATGTGTTGATGGTTCCCCACTACTATAc                                                                                                     <  2:224980/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tgatgaAAAAGCGTATTACCGCGAGAAGAGAAAACGTTGGCGAAAAGAATGGAGTGTGGGGATTGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTTGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCactgatgaatcc      <  2:254827/176‑12 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                          
ACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCCTGCACATTAGGCGTAATTCGAGTGACTTTTCCCCACCATTCGACTATCTTGTTTAGCATATAAAACAAATTACACCGATAACAGCGAATATTACGCTAATGTCGGTTTTAACGTTAAGCCTGTAAAACGAGATGGTAAGATGAGCGTAATGTTTGATCCAGACACGGCGATTTATCCTTTCCCCCCGAAGCCGACGCCGTTAAGCATTGATGAAAAAGCGTATTACCGCGAGAAGATAAAACGTCTGCTAAAAGAACGTAATGCGGTGATGGTTGCCCACTACTATACCGATCCCGAAATTCAACAACTGGCAGAAGAAACCGGTGGCTGTATTTCTGATTCTCTGGAAATGGCGCGCTTCGGTGCAAAGCATCCCGCTTCTACTTT  >  NC_000913/781895‑782332

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