Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913445,5580AC30.6% 70.0 / 11.5 36V370G (GTA→GGA) yajRputative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (20/5);  new base C (0/11);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 7.27e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445424‑445693
                                                                                                                                       |                                                                                                                                         
aaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATAccccc                                                                                                                                      <  1:1243066/140‑1 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTTTTTCATGGGACTGGCGACCGTCAACCACACTGCGGCCCGCATTTCGCCAGCGGGAAAAAcccccc                                                                                                                                     >  2:147310/1‑139 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGGCTGTCAGCCCCACCGCGGGCCGCCTTTCGCCCGCGGGGAAAAAccccc                                                                                                                                      >  2:30080/1‑140 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                     <  1:140739/140‑1 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                     <  2:1039347/140‑1 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGGTCTTTCATGGTACTTGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGGGAGAAATAcccccc                                                                                                                                     >  2:661862/1‑139 (MQ=255)
              gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                       <  1:852007/140‑1 (MQ=255)
                 tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGAGAAAAACCCCCCGACCGTCAAACATGcc                                                                                                                    >  2:335243/1‑140 (MQ=255)
                           cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                          <  2:1068074/140‑1 (MQ=255)
                           aaaTTTTCCCCGGGATTTCAATGCGTAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                          <  2:553972/139‑1 (MQ=255)
                             aTGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTTCGCCCGCGAGAAAAACCCCCCGACCGTCAAACAAGCCGTTAATCCAgcc                                                                                                       >  2:996971/1‑140 (MQ=255)
                                        gATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAAc                                                                                            <  2:1246182/140‑1 (MQ=255)
                                                  cgcAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGACCGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAAAACCCCCGACCGTCAAACATGCCGTTAATCCCGCCCCCCCGCGAACCGCCAATcgc                                                                                  >  2:1210019/1‑140 (MQ=255)
                                                            tgACATACGGCGGTTCTTTTCTGGGACGTGGGACTGGGAAGCACACCGTGGCCCGCACTGCTGCCGCGGGAAATACCCCCTGTCCACGAAACATCACGTTAATCCAGCCTCCCCGCGAACCCCCCATCGCCCCTCCCAGa                                                                        >  2:256636/1‑140 (MQ=255)
                                                            tgACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAAAACCCCCCGACCGTCAAACATGCCGTTAATCCAGCCCCCCAGCGAACCGCCAAACGCCCCCCCAAGa                                                                        >  2:353885/1‑140 (MQ=255)
                                                             gACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCCGTAATCCAGCCGCCCCGCGAACCGCCAATCGCCACGCCAAGaa                                                                       >  1:730740/1‑140 (MQ=255)
                                                              aCATACGGCGGTTTTTTCAGGGTACGGGCGCCTGTCACCCACACTGCGGCCAGCATTGCGCCAGAGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGaaa                                                                      <  2:731168/140‑1 (MQ=255)
                                                                                   gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                 <  1:798966/140‑1 (MQ=255)
                                                                                   gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                 >  1:577747/1‑140 (MQ=255)
                                                                                   gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                 >  2:47962/1‑140 (MQ=255)
                                                                                      cTGGCGACTGTCAGCCACACTGCGGCAAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                              <  1:27310/140‑1 (MQ=255)
                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  1:345591/1‑140 (MQ=255)
                                                                                        ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                            >  2:426339/1‑140 (MQ=255)
                                                                                        ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                            >  2:1097671/1‑140 (MQ=255)
                                                                                         gtgACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           <  1:66115/138‑1 (MQ=255)
                                                                                         gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  2:1036260/1‑140 (MQ=255)
                                                                                         gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  2:109948/1‑140 (MQ=255)
                                                                                         gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                           >  2:180942/1‑140 (MQ=255)
                                                                                           gccTGCCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCGCCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                         <  2:199863/138‑1 (MQ=255)
                                                                                             cTGTCAGCCACCCGGCGGCCAGCATGGCGCCAGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  2:577747/140‑1 (MQ=255)
                                                                                             cTGTCACCCACACGGCGGCCAGCATGGCGCCAGCGAAAAATCCCCCCGGCCCGTCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  2:726419/140‑1 (MQ=255)
                                                                                              tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:31436/1‑140 (MQ=255)
                                                                                              tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:619317/1‑140 (MQ=255)
                                                                                              tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  2:94972/1‑140 (MQ=255)
                                                                                                     cacaCTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTa                               <  1:180942/140‑1 (MQ=255)
                                                                                                        acTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                            >  1:770250/1‑140 (MQ=255)
                                                                                                         cTGCGGCCAGCATGGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACct                           <  1:661862/140‑1 (MQ=255)
                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  1:1166595/1‑140 (MQ=255)
                                                                                                           gCGGCCAGCATGGCCCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgc                         <  2:696557/140‑1 (MQ=255)
                                                                                                            cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  2:523951/1‑140 (MQ=255)
                                                                                                              gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      >  1:454615/1‑140 (MQ=255)
                                                                                                                 agcttggcgcAAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  1:1278217/134‑1 (MQ=255)
                                                                                                                 acaaTTCCGCCAGCAAGAAATCCCCCCGGACCGTAAAACACGCCGTTAACCCACCCGCCCAGCAAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  1:256636/137‑1 (MQ=255)
                                                                                                                 aGCATTCCGCCGCCAAAAATCCCCCCTGGCCCTAAAAACATGCGTTAAACCCGCCCCCCCACCAAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  2:639348/140‑1 (MQ=255)
                                                                                                                  gctttGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAc                  <  1:336261/137‑1 (MQ=255)
                                                                                                                   cATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  1:996971/140‑1 (MQ=255)
                                                                                                                     ttccgcCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct               <  2:454615/137‑1 (MQ=255)
                                                                                                                       gcgcCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  2:213895/140‑1 (MQ=255)
                                                                                                                       gcgcCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAACCACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  1:353885/140‑1 (MQ=255)
                                                                                                                        cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa            >  2:1009189/1‑140 (MQ=255)
                                                                                                                            aGCGAGAAATCCCCCGTGACCGCCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa        <  2:105368/140‑1 (MQ=255)
                                                                                                                             gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt       >  1:121370/1‑140 (MQ=255)
                                                                                                                             gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGGAACCTGCTGGCGACTCTTTACTGAt       >  1:1105618/1‑140 (MQ=255)
                                                                                                                                  aaaaTCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt   <  1:47962/139‑1 (MQ=255)
                                                                                                                                   aaaTACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTg  >  2:189643/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                         
AACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445424‑445693

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