Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5750AC28.1% 47.5 / 14.1 32V555G (GTG→GGG) ybaLputative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (21/2);  new base C (0/9);  total (21/11)
Fisher's exact test for biased strand distribution p-value = 1.96e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501438‑501710
                                                                                                                                         |                                                                                                                                       
cgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCaccacc                                                                                                                                       >  1:491502/1‑140 (MQ=255)
   ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                                    >  2:982295/1‑140 (MQ=255)
               aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                        >  1:360369/1‑140 (MQ=255)
               aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                        >  2:144095/1‑140 (MQ=255)
                    gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                   >  1:72390/1‑140 (MQ=255)
                         cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                              >  2:185156/1‑140 (MQ=255)
                              cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                         >  2:577197/1‑140 (MQ=255)
                                  tCTGGTTTTTGGCTCTCGATGCTTGACGCGCCCCTGGCGCTTCTTATCATCCCTATATATTTTTCATTTTTTCCCCCCGGCAACCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                     <  2:759295/140‑1 (MQ=255)
                                      gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                 >  2:819399/1‑140 (MQ=255)
                                                      gCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTACTTTTTTTTCATTTTTTACACCGGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:871609/140‑1 (MQ=255)
                                                      gCCTGATGCGACGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:550051/140‑1 (MQ=255)
                                                      cccTGGTGCGCCGCTGGCGCCTCTTATCATGCCTACATTTTTTTCATATTTTCCCCCCGGCACCCCCCGTTTCCCCCCTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:642824/139‑1 (MQ=255)
                                                       ccggttGCGCCGCGGGGGCTTTTTATCATGCCTAATTTTTTTTCATTTTTAACACCCCCCACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                <  2:491502/135‑1 (MQ=255)
                                                         tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                              <  1:577197/140‑1 (MQ=39)
                                                          gATGCGACGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGCCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                             <  2:87689/140‑1 (MQ=255)
                                                             gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                          >  1:518580/1‑140 (MQ=255)
                                                             gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                          >  1:173045/1‑140 (MQ=255)
                                                                ccGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                       <  1:982295/139‑1 (MQ=255)
                                                                         cgTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                              <  1:144095/140‑1 (MQ=255)
                                                                          tttttATCGTGCCTACCTTTTTTTCATTTTTCCCACCCGGCACCCCCCGTTTCCCCCCACCCCCCCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                             <  2:665106/137‑1 (MQ=255)
                                                                              tATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                         >  1:42853/1‑140 (MQ=255)
                                                                                tCATCCCTAAATTTTTTTCATATTTCCCATCCGGAACCCACTGTTTCCCCCCTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                       <  2:577375/140‑1 (MQ=255)
                                                                                     ccTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  <  2:42853/140‑1 (MQ=255)
                                                                                     ccTACATTTTTTTCATTTTTTACACCCGCCACCCCCCTTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  <  2:215278/140‑1 (MQ=255)
                                                                                     ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  >  1:167022/1‑140 (MQ=255)
                                                                                     ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  >  1:420260/1‑140 (MQ=255)
                                                                                     ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  >  1:985778/1‑140 (MQ=255)
                                                                                         cATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                              >  2:249250/1‑140 (MQ=255)
                                                                                           tttttttCATTTTTTACACCCGCCACCCCCCGTTTCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                            <  2:72390/138‑1 (MQ=255)
                                                                                             tttttCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                          <  2:420260/140‑1 (MQ=255)
                                                                                                ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                       >  1:495681/1‑140 (MQ=255)
                                                                                                    tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                   >  2:686862/1‑140 (MQ=255)
                                                                                                               cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                        >  1:196057/1‑140 (MQ=255)
                                                                                                                ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  2:618128/1‑140 (MQ=255)
                                                                                                                ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  1:73813/1‑140 (MQ=255)
                                                                                                                  caccccccGTTCACCCCTCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATaa                     <  2:495681/134‑1 (MQ=255)
                                                                                                                        ccttttCCCCCGTTCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg               <  2:452348/137‑1 (MQ=255)
                                                                                                                        ccttttCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg               <  2:285118/137‑1 (MQ=255)
                                                                                                                        cctttccccccGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg               <  2:73813/133‑1 (MQ=255)
                                                                                                                           tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg            >  2:187388/1‑140 (MQ=255)
                                                                                                                              cccccGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCa         <  1:827403/139‑1 (MQ=255)
                                                                                                                                cccGTTACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat       >  1:452348/1‑140 (MQ=255)
                                                                                                                                   gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    <  1:185156/140‑1 (MQ=255)
                                                                                                                                     caccacCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct  >  2:135253/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501438‑501710

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