Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5640AC24.2% 65.8 / 17.1 33*559E (TAA→GAA) ybaLputative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (22/3);  new base C (0/8);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 1.19e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGG  >  NC_000913/501428‑501701
                                                                                                                                        |                                                                                                                                         
aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAccc                                                                                                                                        <  1:115568/140‑1 (MQ=255)
 aGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcccc                                                                                                                                       <  1:1385862/140‑1 (MQ=255)
    tATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                    <  2:698035/140‑1 (MQ=255)
      tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                  >  2:647196/1‑140 (MQ=255)
      tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCACGTcac                                                                                                                                  >  2:373698/1‑140 (MQ=255)
                              gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                          >  1:601453/1‑140 (MQ=255)
                              gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                          >  2:1288954/1‑140 (MQ=255)
                                 ttCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                       <  2:768400/140‑1 (MQ=25)
                                  tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                      >  1:370013/1‑140 (MQ=255)
                                   cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                     >  2:521065/1‑140 (MQ=255)
                                   cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                     >  1:1010562/1‑140 (MQ=255)
                                   cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                     >  2:1034273/1‑140 (MQ=255)
                                     ccccgcCTCTGCCATTTGCCTCTCCATGCTTGAGCCCCCGCGGGCGTTTCTTACCATGCCTACATTTTTTTCATTTTTTCCACCCGCCACCCGCCGTTTCCCCCGCCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                   <  2:838193/138‑1 (MQ=255)
                                                                gCCTGATGCGCCGCGGGCGCTTTTTATCATGCCTACTTTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                        <  2:665826/140‑1 (MQ=255)
                                                                   tGATGCGACGCGGGCGCTTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                     <  2:1112837/140‑1 (MQ=255)
                                                                    gATGCGACGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                    <  2:601453/140‑1 (MQ=255)
                                                                    gATGCGACGCGGACGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCGTTTTCCCCCGCCCCCACCACGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                    <  2:1010562/140‑1 (MQ=255)
                                                                           cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                             >  1:1112837/1‑140 (MQ=255)
                                                                                       ttATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                 >  2:1271384/1‑140 (MQ=255)
                                                                                             tGCCTACATTTTTTTCATATTTTACACCCGGCACCCACCGTTTCCCACGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                           <  1:373698/140‑1 (MQ=255)
                                                                                               ccTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                         <  2:441202/140‑1 (MQ=255)
                                                                                               ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                         >  1:441202/1‑140 (MQ=255)
                                                                                               ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTCCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                         >  1:160313/1‑140 (MQ=255)
                                                                                                cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                        >  1:989532/1‑140 (MQ=255)
                                                                                                cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                        >  1:958492/1‑140 (MQ=255)
                                                                                                  caatttttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTCCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                      <  2:160313/135‑1 (MQ=255)
                                                                                                  aCATTTTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                      <  2:184544/140‑1 (MQ=255)
                                                                                                    atatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                    >  1:130161/1‑140 (MQ=255)
                                                                                                       tttttCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                 <  1:1083316/140‑1 (MQ=255)
                                                                                                         tttCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                               <  2:958492/140‑1 (MQ=255)
                                                                                                           tCATATTTTACATCCGGCACCCACCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                             <  1:793460/140‑1 (MQ=255)
                                                                                                                 tttCCACCCGCCAACCACGTTTTACCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                       <  1:1034273/140‑1 (MQ=255)
                                                                                                                 tttACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                       <  1:521065/140‑1 (MQ=255)
                                                                                                                   tACACCCGGCACCCCCCGTTTCCCCCGTCACCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                     <  1:1041585/140‑1 (MQ=255)
                                                                                                                          ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt              >  1:748350/1‑140 (MQ=255)
                                                                                                                          ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt              >  1:365889/1‑140 (MQ=255)
                                                                                                                          ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt              >  2:955305/1‑140 (MQ=255)
                                                                                                                           gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa             >  2:758332/1‑140 (MQ=255)
                                                                                                                                  ccttttCCCCCCTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg      <  2:373725/137‑1 (MQ=255)
                                                                                                                                    ttttCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg    <  2:292605/139‑1 (MQ=255)
                                                                                                                                     tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg   >  1:373725/1‑140 (MQ=255)
                                                                                                                                      ttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg  >  2:301033/1‑140 (MQ=255)
                                                                                                                                      ttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg  >  1:701521/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGG  >  NC_000913/501428‑501701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑