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

AGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGG  >  NC_000913/445420‑445696
                                                                                                                                          |                                                                                                                                          
aGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTTTTTCATGGGACTGGCGACTGTCTGCCTCACTGCGGCCCGCATTTCGCCAGCGAGAAAAAc                                                                                                                                           >  2:9220/1‑140 (MQ=255)
aGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGGCTGTCAGCCACACTGCGGCCCGCATTGCGCCAGCGAGAAATAc                                                                                                                                           >  1:545956/1‑140 (MQ=255)
aGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATAc                                                                                                                                           >  1:1186621/1‑140 (MQ=255)
   aaaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATAcccc                                                                                                                                        <  2:1186621/140‑1 (MQ=255)
     gacagcTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                      <  1:649989/135‑1 (MQ=255)
     aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                      <  2:700044/140‑1 (MQ=255)
     aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                      <  2:588020/140‑1 (MQ=255)
       gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGa                                                                                                                                    <  1:423907/140‑1 (MQ=255)
                 cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                          <  1:636910/140‑1 (MQ=255)
                 cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                          <  1:9220/140‑1 (MQ=255)
                     tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCCCACTGCGGCCCGCATTGCGCCCGCGAGAAAAACCCCCCGACCGTCAAACATGcc                                                                                                                      >  2:229616/1‑140 (MQ=255)
                     tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCATCAAACATGcc                                                                                                                      >  1:676981/1‑140 (MQ=255)
                     tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGGCTTTCAGCCACACTTCGGCCAGCATTGCGCCAGCGAGAAATACCCCCCGACCGGCCAACCTgtc                                                                                                                      >  2:649989/1‑138 (MQ=255)
                      cGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCg                                                                                                                     <  1:198035/140‑1 (MQ=255)
                       gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTTAGCCACACCGCGGCCAGCATTGCGCCCACGGGGAAAACCCCCTGACCCTCAAACAAGCCgc                                                                                                                    >  2:977831/1‑139 (MQ=255)
                              gCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTTTTTTATGGTACTGGCGGCTGTCCGCCACACTTCGGCCAGCAATGCGCCAGCGAGAAATAACCCCCCACCGCCCAACAAGCCGTTTATccc                                                                                                             >  2:198035/1‑139 (MQ=255)
                               cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCATCAAACATGCCGTTAATCCAg                                                                                                            <  2:676981/140‑1 (MQ=255)
                                            gATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCCGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAAc                                                                                               >  1:1144159/1‑140 (MQ=255)
                                                      cgcAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCCGACCGTCAAACATGCCGTTTATCCAGCCGCCCCGCGAACCGCCAATcgc                                                                                     >  1:342619/1‑140 (MQ=255)
                                                        cAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGGGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                                   <  1:451391/140‑1 (MQ=255)
                                                           aCTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAAAACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgccacgc                                                                                >  1:815591/1‑140 (MQ=255)
                                                                  aCATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAAAACCCCCCGACCCTCAAACATGCCGTTAATCCAGCCGCCCCGCGAACCGCCCAACGCCACGCCCAaaaa                                                                         >  2:658944/1‑140 (MQ=255)
                                                                    aTAGGGCGTTTCTTTCATGGTACGGGCGACTGTCAGCCACCCGGGGGCCAGCATCGCCCCGGGAAGAATTCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAAct                                                                       <  2:250971/140‑1 (MQ=255)
                                                                           ggTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                                <  1:658944/140‑1 (MQ=255)
                                                                            gTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTGCGCCAGCGAGAAAAACCCCCCGACCGTCAAAAATGCCGTTAATCCAGCCGCCCCACGAAACCCCCATCCCCCCCCCAAGAAAatggctgggg                                                               >  2:1236390/1‑138 (MQ=255)
                                                                            gTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCCGCGAAACGCCCAACGCCACGCCAAGAAActggctgggg                                                               >  2:1226656/1‑138 (MQ=255)
                                                                                       gTACTGGCGACTGTCAGCCACCCTGCGGCCGGCTTTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                    <  1:49530/140‑1 (MQ=255)
                                                                                         ccTGGCGCCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                                  <  1:229616/139‑1 (MQ=255)
                                                                                          cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                                 >  1:376818/1‑140 (MQ=255)
                                                                                            ggCGCCTGCCACCCACACTGCGCCCAGCATTGCCCCAGCAAAAAATCCCCCCTGCCCGCCAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                               <  2:810961/140‑1 (MQ=255)
                                                                                            ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                               >  1:810961/1‑140 (MQ=255)
                                                                                            ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                               >  2:137183/1‑140 (MQ=255)
                                                                                                aCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                           >  1:1025289/1‑140 (MQ=255)
                                                                                                 cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                          <  1:691478/140‑1 (MQ=255)
                                                                                                    tCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCt                                       >  1:933583/1‑140 (MQ=255)
                                                                                                       cccACACGGCGCCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                                    <  1:581551/139‑1 (MQ=255)
                                                                                                            acTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                               >  1:488069/1‑140 (MQ=255)
                                                                                                              tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                             >  1:9310/1‑140 (MQ=255)
                                                                                                                     aGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                      <  1:137183/140‑1 (MQ=255)
                                                                                                                        aTTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActc                   >  1:1246273/1‑140 (MQ=255)
                                                                                                                         ttGCGCCAGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct                  <  2:815591/140‑1 (MQ=255)
                                                                                                                           gcgcCGGCAAAAAATCCCCCCTCCCCATAAACAATGCCGTAAAACAGGCCGCCCGGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt                <  2:308125/140‑1 (MQ=255)
                                                                                                                            cgcCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                         <  1:47118/130‑1 (MQ=255)
                                                                                                                            cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                         >  2:47118/1‑130 (MQ=255)
                                                                                                                             gcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActcctgtc              >  1:468675/1‑135 (MQ=255)
                                                                                                                             cccAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTAc              <  1:1226656/139‑1 (MQ=255)
                                                                                                                             cccAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTAc              <  1:1236390/139‑1 (MQ=255)
                                                                                                                                 gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt          >  1:308125/1‑140 (MQ=255)
                                                                                                                                 gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt          >  2:1120003/1‑140 (MQ=255)
                                                                                                                                     aaaaTCCCCCCTGCCCGCAAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt      <  2:9310/139‑1 (MQ=255)
                                                                                                                                         tACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAgg  >  1:183231/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                          
AGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGG  >  NC_000913/445420‑445696

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