Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG34.4% 132.3 / 12.1 61D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (3/37);  new base G (21/0);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 1.66e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTA  >  NC_000913/827273‑827742
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
tCGTGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAAc                                                                                                                                                                                                                                             >  1:21820/1‑239 (MQ=255)
tCGTGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCaa                                                                                                                                                                                                                                              >  1:591317/1‑238 (MQ=255)
                       ttGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaa                                                                                                                                                                                                                      <  2:591317/239‑1 (MQ=255)
                        tGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCCGAAGTCGGtgttccgggaaa                                                                                                                                                                                                                     >  2:649151/1‑228 (MQ=255)
                        tGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGCAAGGCCGGAAGACGGTtcgtcgcgga                                                                                                                                                                                                                      >  2:1145007/1‑234 (MQ=255)
                        tGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAGACGGGTCGTGGGGaaa                                                                                                                                                                                                                     >  2:1369491/1‑239 (MQ=255)
                        tGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCGGGAGGCGGGtctgcgcggaa                                                                                                                                                                                                                     >  2:266503/1‑230 (MQ=255)
                           ttcATTGTGGGATCGGGTTGCTCATCGTTAGCCGACTTTGCTTTCAAAGCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                                                                                                                                                                                                   <  2:1222529/236‑1 (MQ=255)
                                 gtggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                            <  2:1075816/239‑1 (MQ=255)
                                  tggGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                            <  1:158198/238‑1 (MQ=255)
                                   gggATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCg                                                                                                                                                                                                          <  2:1441742/239‑1 (MQ=255)
                                               cTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCGGTTCGTTGAGGAAAAAAATGTTCGGGTCACGccac                                                                                                                                                                                               >  1:930886/1‑235 (MQ=255)
                                               cTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTCCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCAGCCATAGAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTCCCATGCTGTGGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAGCGCCGGAAGTCGGTTCGTCGAGAAACAGAACGTCCGGTACATGCATCa                                                                                                                                                                                              <  2:592460/239‑1 (MQ=255)
                                                    cGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAATCCCGGGAATGGGTTCGTCGGAGAACAGAATATCTCGTgcatgcatgaccggc                                                                                                                                                                                         >  2:809796/1‑232 (MQ=255)
                                                         gCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAACCCGGGAGGCGGTTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCACCggaaaag                                                                                                                                                                                    >  1:1078478/1‑232 (MQ=255)
                                                              cTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGACGGAAGGCGGGTTGTCCAGAAAACGAATTGCCGCGTCATGCATCAACGCACAAgccacc                                                                                                                                                                               >  1:981038/1‑236 (MQ=255)
                                                                gtgCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcc                                                                                                                                                                            <  1:1145007/239‑1 (MQ=255)
                                                                     ttCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCg                                                                                                                                                                      <  1:1046597/239‑1 (MQ=255)
                                                                                                   aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGt                                                                                                                                       <  1:649151/239‑1 (MQ=255)
                                                                                                               tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAgg                                                                                                                           <  1:266503/238‑1 (MQ=255)
                                                                                                                 cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                                                                                                                                                                     <  1:910020/196‑1 (MQ=255)
                                                                                                                           aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCAAGGCGGAAAGCGGGTTGTCGAGAAAAAGAAAGTCCGGGTTATGCATCAGCGAAAAGAGCAGCGCCCGCCGCCGGTTTAAAAATAAAGGGCGGTCTCCGGGGGCGTGGGGGGGGAGAAttttttg                                                                                                              >  2:584847/1‑236 (MQ=255)
                                                                                                                           aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTTGTCGAGAAACAGAAAGTCCGGGTCATGCATACGCGAACAAGCCCGCGCCAGCCGCCGTTTAAAAc                                                                                                                                                          >  2:910020/1‑196 (MQ=255)
                                                                                                                            tGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGCAAGCGGGTCGTCGAGAAAAAAAATGTCCGGGTCATTCATCCACGAAACAGCCAGCCCCCGCCCCCTGTTTAAAACCAAAAGCCAGTTCATCGGGGGGGGGGGGGGGGGAACTTTCCAg                                                                                                              >  2:1423494/1‑239 (MQ=255)
                                                                                                                               ggTCGAAATTTTCCGCTTCCTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCg                                                                                                          <  2:929714/239‑1 (MQ=255)
                                                                                                                                 tCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGa                                                                                                         <  2:272081/238‑1 (MQ=255)
                                                                                                                                   gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACg                                                                                                      <  2:21820/239‑1 (MQ=255)
                                                                                                                                      atcttCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCGTGCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCt                                                                                                   <  2:443394/236‑1 (MQ=255)
                                                                                                                                               ttCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGGGAAACAGAAAGTCCGGGTCCTTGATCACCCGAAAAGCCACGCCCCACCCCCGTTTTAAAACCAATAGCAGTTAATCGGGGGGCGGGGGGGGGGTACTTTTTAGGCCGAAGGccccccctctcg                                                                                          >  2:796252/1‑229 (MQ=255)
                                                                                                                                                    ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGTGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgccagc                                                                                                                                                                        <  1:34099/158‑1 (MQ=255)
                                                                                                                                                    ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                            <  1:653249/124‑1 (MQ=255)
                                                                                                                                                    ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATg                                                                                                                                                    <  1:623741/177‑1 (MQ=255)
                                                                                                                                                              ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTGAACGCCGGGAGGCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAAAAGGGCCGCGCCAGCCGCTGTTTTAAAACTAAAGGCAGTTCATc                                                                                                                                          >  2:623741/1‑177 (MQ=255)
                                                                                                                                                              ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCTTGc                                                                                                                                                                                                  >  2:653249/1‑124 (MQ=255)
                                                                                                                                                              ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTTGTCGGGAAAAAGAAAGTCCGGGTCATTCCTCCCCGAAAAGGCCCACGCCCGCCCCCGTTTAAAACCTAATTGCCGGTCCTCTGTGGGGGTGGGGGGGGGGATTTTTTGGCCGCACGCCCCGCCCCTTCGGGGGGTTTTTTTGt                                                                           >  2:665822/1‑239 (MQ=255)
                                                                                                                                                                   tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCg                                                                                                                                                                      <  1:462708/145‑1 (MQ=255)
                                                                                                                                                                     aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGt                                                                                                                                       <  1:1602168/173‑1 (MQ=255)
                                                                                                                                                                     aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGAttttt                                                                               <  1:894486/228‑1 (MQ=255)
                                                                                                                                                                      ccATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCCACCCCGGGAGGCGGGTTGTCGAGAAAACAAAAGTCCGGGTCATGCATCAACGGAAAAGCCAGCGCCACCCGCGGTGTTAAAACCAATAGCAGTTCATCGGGGGCGGGGGGGGCGGTCGTTTTTTGGGCGCACGCCCCCCTCATCCGGGGGGGTTTTTTGTTCTtcgc                                                                    >  2:1565001/1‑238 (MQ=255)
                                                                                                                                                                         tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGGCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGa                                                                <  2:782900/239‑1 (MQ=255)
                                                                                                                                                                             cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACcgc                                                            <  2:1584118/239‑1 (MQ=255)
                                                                                                                                                                               tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGGGACAGGTCTGCGGGAAACAGAAATTTCTGTTACTATATCAGCGCACCGGCCCGCGCCAGCCGCTGTGTTAAAACTAAAGACAGTTTCTAGGGGGGGAGGGGGGGGGTACTTTTTATGCCCGAAGCCCCCCTCCCTCGGGGGATGTTTTTGTTCtttg                                                                     >  2:894486/1‑225 (MQ=255)
                                                                                                                                                                                  cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCa                                                                                                                                                                                                     <  1:1276416/101‑1 (MQ=255)
                                                                                                                                                                                  cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAg                                                       <  1:665822/239‑1 (MQ=255)
                                                                                                                                                                                       ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGt                                                                                                                                                                                                        <  1:871507/93‑1 (MQ=255)
                                                                                                                                                                                         ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATaca                                                <  1:1496804/239‑1 (MQ=255)
                                                                                                                                                                                         ccctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATaca                                                <  2:1390922/237‑1 (MQ=255)
                                                                                                                                                                                            tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                                                                                                            <  1:1172572/84‑1 (MQ=255)
                                                                                                                                                                                            tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCCGGAGGCGGGTCGTCGGGAAAAAGAATGTCCGGGTTATTCATCcgcg                                                                                                                                                                                          >  2:1276416/1‑101 (MQ=255)
                                                                                                                                                                                                 tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGCGAAACAGAATTTTCGGTGCGTGCATcc                                                                                                                                                                                              >  2:871507/1‑92 (MQ=255)
                                                                                                                                                                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATg                                                                                                                                                                                                   <  1:356536/86‑1 (MQ=255)
                                                                                                                                                                                                    tGTTGATGTGCAGCGAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTGCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAgg                                                                                                                           <  1:741508/153‑1 (MQ=255)
                                                                                                                                                                                                      ttGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTGACGCCCGAAGGCGGTTCTTCGGGAAACAGAATGACCGGGTCATGc                                                                                                                                                                                                  >  2:1172572/1‑84 (MQ=255)
                                                                                                                                                                                                         aTGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGtt                                                                                                                                                                                                                               <  1:919330/52‑1 (MQ=255)
                                                                                                                                                                                                             gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCCGGAGACGGTTCGTTGAGAAACAGAAATCCCGGTTCATTCAACAACGaa                                                                                                                                                                                        >  2:356536/1‑86 (MQ=255)
                                                                                                                                                                                                                        tCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAg                                                       <  1:1374315/201‑1 (MQ=255)
                                                                                                                                                                                                                        tCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGAc                 <  1:182825/239‑1 (MQ=255)
                                                                                                                                                                                                                           cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGt              <  1:205599/239‑1 (MQ=255)
                                                                                                                                                                                                                              cggGTGAGGGGGTCAACTCCGGAAGGCGGGTCGGCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCTGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCa            <  2:864057/238‑1 (MQ=255)
                                                                                                                                                                                                                              cggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCa            <  1:638593/238‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGGGGGAGGGGGTAATTTTCAGGCCGAACGCCTCGCTCATGCGGGGGATTTTTTCGGTCTGCGCCCGACCGCGGAAGCCATAAACACCAGAGAAAAAGCGGAAATTATGGTCGAACGTCAGGTTAACGTa  >  2:638593/1‑238 (MQ=255)
                                                                                                                                                                                                                                          tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACcgc                                                            >  2:612405/1‑178 (MQ=255)
                                                                                                                                                                                                                                          tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACcgc                                                            <  1:612405/178‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
TCGTGGATCAACTGAATAAAGGCTTGCTCCATGGTGGGATCGGGTTGCTCATCGTTAGCCGACTGTGCTTTCAAATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGGAGATTTTTTCGTTCTGCGCCCGACCGCGTAAGCCATACACACCAGAGAAAAAGCGTAAATTCTGTTCGACCGTCAGGTTACCGTA  >  NC_000913/827273‑827742

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