Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT28.6% 102.8 / 17.3 56V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (4/36);  new base T (15/1);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 3.47e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCT  >  NC_000913/803622‑804116
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
tATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGTGGCCCAGCGTTATGGCGGTATTTTTTTGc                                                                                                                                                                                                                                                            >  1:680265/1‑251 (MQ=255)
          tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCat                                                                                                                                                                                                                                                    <  1:96567/249‑1 (MQ=255)
          tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCat                                                                                                                                                                                                                                                    <  1:436575/249‑1 (MQ=255)
                       ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATGGCGGTcatccat                                                                                                                                                                                                                                    <  1:549163/251‑1 (MQ=255)
                        cGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccatc                                                                                                                                                                                                                                   <  1:222697/251‑1 (MQ=255)
                           aaaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAat                                                                                                                                                                                                                                <  1:741334/251‑1 (MQ=255)
                                  gATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCCCATATTGCGGTCATCCCTCAATATATGCg                                                                                                                                                                                                                         >  1:7083/1‑251 (MQ=255)
                                   aTAAAATTAAAATTAAGTTTGATTCTTTAGCCTCTCACGATATTACCTTTGTTGGTATGGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGt                                                                                                                                                                                                                        <  2:253346/251‑1 (MQ=255)
                                     aaaattAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCCATTGCGGTCATCCATCAATATATTCGTgg                                                                                                                                                                                                                      >  1:476161/5‑250 (MQ=255)
                                        attAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCATATTGCGGGCCTCCCTCAatatattggtgag                                                                                                                                                                                                                     >  1:300209/2‑242 (MQ=255)
                                          aaaaaTTAATTTTGATTAAGTAGCCTCTCACGATAGTACATTTGTTGGTAGTGTCCAGTACGCTAAAGCGTCCGGTATGCGACGTTGCCCGCTGCCGAATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGGACTATTAACGGTGATGGCCATTTTTTTGGGTTATCGGGGGCCCCGCGTTATGGCGGTATTTTGGTGCCTCCGCATATTGCGGGCAGCCATCAATATATGCGTGAgatga                                                                                                                                                                                                                 <  2:315761/250‑1 (MQ=255)
                                                          tCATTAGCCTTCCACGATATTACCTTTGTGGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGGTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCa                                                                                                                                                                                                 <  2:421212/250‑1 (MQ=255)
                                                               aGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCTCATATTGCTGTCATCCATCa                                                                                                                                                                                                                                  <  1:307920/213‑1 (MQ=255)
                                                                    ctcACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTGTCTCCCCCTATTTTGGTCCTCCCACAATATATTCGTGAGATGATGGcgcggggcggcgaaaaatgtctc                                                                                                                                                                                       >  2:134481/1‑229 (MQ=255)
                                                                      cACGATACTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCGCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTc                                                                                                                                                                                      <  1:784171/249‑1 (MQ=255)
                                                                        cGATATTGCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGt                                                                                                                                                                                 <  1:140847/251‑1 (MQ=255)
                                                                            atTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAg                                                                                                                                                                              <  1:170886/250‑1 (MQ=255)
                                                                                ccTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTGCGGGCATCCATCAATATATTCGTGAGAAGATGGCAGGGGGCGGCAAAAtaaatcccggggcagacagg                                                                                                                                                                          >  1:396114/1‑232 (MQ=255)
                                                                                    tgttgGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTTCCTCCCCATATTGCGGGTATCCATCACTATATGCGTGGGGTGATGGCCGGCGGCGGCAAAAAGATCCTCGGGGCAGACAGcccc                                                                                                                                                                     >  1:209494/1‑248 (MQ=255)
                                                                                         gTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCTCCCCCTTTTGCGGGGCTCCATCAATATAATCGTGCGGTGATGGCGGGGGGCGGGAAAAAGAttccccggtgaggacgcacccacccc                                                                                                                                                                 >  2:127788/1‑226 (MQ=255)
                                                                                           aTTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGtt                                                                                                                                                              <  1:239955/250‑1 (MQ=255)
                                                                                                 catcacGCTAAAGCGTCCGGTATGGAACTTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGTCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGGGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGc                                                                                                                                                       <  2:442590/246‑1 (MQ=255)
                                                                                                   gACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATCCCTCAATATATTCGTGGGATGATGGCGGGGGGCGGCAAAAAGATCCTCGGGGCAGAAAGCCCCACCCGTTACGGtaga                                                                                                                                                      >  1:168820/1‑247 (MQ=255)
                                                                                                   gACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGt                                                                                                                                                         <  2:168820/247‑1 (MQ=255)
                                                                                                      cGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAg                                                                                                                                                 <  2:158947/251‑1 (MQ=255)
                                                                                                       gCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTCTTGCGGTCATCCATCAATATATGCGTGGGATGAAGGGAGGGGGCGGGCAAATGATCCTCGGGGCAGACAGCCACACCCCTTACGGGGCAtttgg                                                                                                                                                 >  1:642454/1‑248 (MQ=255)
                                                                                                          aaaGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCAtt                                                                                                                                                    <  1:166682/245‑1 (MQ=255)
                                                                                                          aaaGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCTCCGCATTTTTGGGGCATCCCACAATATATGTGTGGGGTGATGGTCGGGGGGGGCAAAATTATTCTCGGGGTCGAACACCACACCCGTTACgggggtctatatcg                                                                                                                                              >  2:166682/1‑239 (MQ=255)
                                                                                                                   ggTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCTCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGGCAATGGCAGt                                                                                                                                   <  2:736407/251‑1 (MQ=255)
                                                                                                                    gTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTc                                                                                                                                  <  2:300209/251‑1 (MQ=255)
                                                                                                                          aaCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTGCTCCGCATATTGCGGGCATCCATCAATATATTCGGGAGATGATGGCAGGGGGGGGCAAAATGATTCTCGGGGCAGAAAGCCACACCCGTTACGGGGGATTTGGGGCAATGGCCGTCgggggg                                                                                                                            >  1:51605/1‑247 (MQ=255)
                                                                                                                             gTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTTCCCCTATTGTCGGCAATCATCAATATTTGGGTGGGGTGGTGGCAAGCGGCGGGAAAAAGATTCTCGGGGCAGACAGGCCAACCCGATACGGGGCATTTAGGGAAATgggacgtggggagggg                                                                                                                         >  1:174002/1‑237 (MQ=255)
                                                                                                                                       gCCGTATGTGCGGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGt                                                                                                               <  2:327897/251‑1 (MQ=255)
                                                                                                                                         cGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGtt                                                                                                              <  2:184917/250‑1 (MQ=255)
                                                                                                                                         cGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCATATTGCGGGCCTCCCTCCATATATTAGTGAGATGATGGCAGGCGGCGGCGAAAATATCCTCGGGGCAGACAGCCACACCCCTTACGGGGGATTTGGGACAATGGCAGGCgggggggggggggggggggtt                                                                                                             >  1:184917/1‑232 (MQ=255)
                                                                                                                                                   gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTGCTCCGCCAATTTCGGGCATCCCTCAATATTTTCGGGGGGTGAAGGCAGGGGGGGGGAAAAAGATCCTCGGGGCCGGAAGCCACACCCCCTTCGGTGTATTTGGGGCAATGGCAGTGGGGGAGGGgggcgggggggtggggaaaca                                                                                                     >  2:665060/1‑234 (MQ=255)
                                                                                                                                                   gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgc                                                                                                   <  2:381519/251‑1 (MQ=255)
                                                                                                                                                             cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCCGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGa                                                                                         <  1:649212/251‑1 (MQ=255)
                                                                                                                                                              aTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTAGCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGac                                                                                        <  2:21739/251‑1 (MQ=255)
                                                                                                                                                                    cactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTgg                                                                                  <  1:500176/251‑1 (MQ=255)
                                                                                                                                                                     actcTGCGCCGTCGGCGGCCCTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCCGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGTCACCTggg                                                                                 <  2:492751/251‑1 (MQ=255)
                                                                                                                                                                          gcgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATc                                                                            <  1:121102/251‑1 (MQ=255)
                                                                                                                                                                          gcgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTGTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCGCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAGGGCAGGCGGCGGCAAAATGATCCGCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGatat                                                                             <  2:291229/250‑1 (MQ=255)
                                                                                                                                                                             ccGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGAc                                                                         <  2:7083/251‑1 (MQ=255)
                                                                                                                                                                                        aCTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCCTATTGCGGGCATCCCTCAATATATGCGTGAGGTGATGGCCGGCGGCGGCAAAAAGATCCTCGGGGCAGACAGCCACACCCGTTACGGGGCATTAGGGAACAAGGGAGTCGGGGAGGGGGGGGGGGGGGTGGTAAAAAAGCTGCTTAATGACACCTGGGGTATCGACTATTCGGGCgg                                                              >  1:637695/1‑250 (MQ=255)
                                                                                                                                                                                             tAACGGTGATGACTATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGGGGCGGCAAAATGATGCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGCGCCAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGtt                                                          <  2:222082/250‑1 (MQ=255)
                                                                                                                                                                                              aaCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTCTTTTTTTGCCTCCGCATATTGCGGGCCTCCATCAATATATGCGTGAGATGATGGCAGGGGGCGGCAAAAAGATCCTCGGGGCAGACAGCCACACCCGTTTCGGTGCCTTAGGGACAAAGGCAGTCGGGGAGGGCGGCGGGGAGTTGGTAAAACAGATGCTTAATGACACCTGGGATATCGACTATACGGGCgggggtgg                                                        >  1:271198/1‑244 (MQ=255)
                                                                                                                                                                                                 ggTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCTTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGTGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGCCTATCCGGGCGTGGTTGCGGt                                                     <  2:251999/251‑1 (MQ=255)
                                                                                                                                                                                                               ttttGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGa                                        <  1:701696/250‑1 (MQ=255)
                                                                                                                                                                                                                    gTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCt                                  <  1:127788/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTTTATCGGCGGCCCAGCGTTAGGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCt                                  <  2:586800/251‑1 (MQ=255)
                                                                                                                                                                                                                          cggcggCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGcc                             <  1:568363/250‑1 (MQ=255)
                                                                                                                                                                                                                             cggcCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGt                           <  1:628525/249‑1 (MQ=255)
                                                                                                                                                                                                                                        tATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCcac               <  1:311683/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ggTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtgt        >  1:8271/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    tttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCt  >  1:399797/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                             
TATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCT  >  NC_000913/803622‑804116

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