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

TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCG  >  NC_000913/501426‑501697
                                                                                                                                          |                                                                                                                                     
ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGG‑‑CTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAccc                                                                                                                                    >  1:167966/1‑140 (MQ=255)
      tATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                <  2:956183/140‑1 (MQ=255)
        tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                              >  1:78850/1‑140 (MQ=255)
             ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCt                                                                                                                         >  1:3225/1‑140 (MQ=255)
                      gCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgccg                                                                                                                >  1:103796/1‑140 (MQ=255)
                        cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                              >  2:856764/1‑140 (MQ=255)
                           aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                           >  2:224104/1‑140 (MQ=255)
                                gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                      >  2:787970/1‑140 (MQ=255)
                                gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                      >  1:251902/1‑140 (MQ=255)
                                   tccACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                                                   <  2:1201847/138‑1 (MQ=25)
                                    tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                  >  1:279679/1‑140 (MQ=255)
                                     cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATCTTTTTCATTTTTTACATCCGGCAACCACCGTTTACCCCGTCACCCCCTCACCCCCCGGTGGCGTTTCcagca                                                                                                 >  2:355928/1‑140 (MQ=38)
                                     cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                 >  1:51395/1‑140 (MQ=255)
                                     cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                 >  2:1068600/1‑140 (MQ=255)
                                       cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                               >  2:1115817/1‑140 (MQ=255)
                                         ccacATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCACCGTTTACCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                             <  1:659019/137‑1 (MQ=255)
                                            cATCGGGCTTTTGCCTCTCGATCCCTGATGCGCCGCGGGCCCGTTTTATCATCCCTACATTTTTTTCATTTTTCACATCCGCCACCCCCCGTTTCCCCCGCCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                          <  2:580554/140‑1 (MQ=255)
                                                   tttttGGCTCTCGATCCATGATGCGCCGCGGGCGCGTCTTATCTTGCCTACATTTTTTCCATATTTTACACCCGGAACCCACCTTTTCCCCCGCCCCCCCCCCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa                                                                                   <  2:103796/138‑1 (MQ=255)
                                                       tGGCCCTCGATCCCTGATGCGCCGCGGGCGCGTCTTATCATCCCTACTTTTTTTTAATATTTTACACCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                               <  2:1054115/140‑1 (MQ=255)
                                                            ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                          >  2:1153036/1‑140 (MQ=255)
                                                                  gCCTGATGCGACGCTGGGGCGTCTTATCATGCCTACTTTTTTTTCATTTTTAACACCCGCCACCCCCCGTTCCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                    <  2:167966/140‑1 (MQ=255)
                                                                      gATGCGCCGCGGGCGCTTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                <  2:354999/140‑1 (MQ=255)
                                                                       gtgggccgcgGGCGTGTCTTTTCATCCCTATATTTTTTTCTTTTTTTCCACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                               <  2:922247/134‑1 (MQ=255)
                                                                                gggCGCTTTTTATCATGCCTATATTTTTTTTATATTTAACATCCGGAAACCACCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGTAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                      <  2:813593/139‑1 (MQ=255)
                                                                                gggCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                      <  2:3225/139‑1 (MQ=255)
                                                                                       tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                               >  2:247154/1‑140 (MQ=255)
                                                                                            tcagccttACATTTTTTTCATATTTCACACCCGCCAACCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                          <  2:251902/133‑1 (MQ=255)
                                                                                            tCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                          <  1:1229963/140‑1 (MQ=255)
                                                                                             cATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                         <  2:51395/140‑1 (MQ=255)
                                                                                                 ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                     >  1:1228404/1‑140 (MQ=255)
                                                                                                 ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                     >  1:1268546/1‑140 (MQ=255)
                                                                                                  cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAATTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                    >  2:155101/1‑140 (MQ=255)
                                                                                                   tACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                   <  2:279679/140‑1 (MQ=255)
                                                                                                   tACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                   <  1:856764/140‑1 (MQ=255)
                                                                                                     cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                 >  1:411238/1‑140 (MQ=255)
                                                                                                            ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                          >  2:221523/1‑140 (MQ=255)
                                                                                                             tCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                         <  1:273700/140‑1 (MQ=255)
                                                                                                                    ttccacccGGCACCCACCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                  <  2:78850/134‑1 (MQ=255)
                                                                                                                     tACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                 <  1:1115817/140‑1 (MQ=255)
                                                                                                                         cccGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt             <  2:205612/139‑1 (MQ=255)
                                                                                                                            ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt          >  1:1107308/1‑140 (MQ=255)
                                                                                                                             gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa         >  1:931347/1‑140 (MQ=255)
                                                                                                                                    ccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg  >  2:731430/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                     
TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCG  >  NC_000913/501426‑501697

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