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

TAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445419‑445693
                                                                                                                                           |                                                                                                                                       
tAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTCGCGACTGTCAGCCACACTCCGGCCAGCATTGCGCCCGGGAGAAATa                                                                                                                                         >  1:1256429/1‑140 (MQ=255)
      aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                   <  2:926315/140‑1 (MQ=255)
      aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                   <  1:977277/140‑1 (MQ=255)
                aaCGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCaaa                                                                                                                         <  2:879612/140‑1 (MQ=255)
                aaCGCCTCGTTTGCGGCAATGTTCGCCGAGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCCTGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAATACCCCCCGACCGTCaaa                                                                                                                         >  1:900738/1‑140 (MQ=255)
                    ccTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATg                                                                                                                     >  2:364167/1‑140 (MQ=255)
                    ccTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATg                                                                                                                     >  2:730744/1‑140 (MQ=255)
                      tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTGCGCCCGCGGGAAATACCCCCTGACCGTCAAACATGcc                                                                                                                   >  2:1182262/1‑140 (MQ=255)
                      tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                      <  2:906956/137‑1 (MQ=255)
                      tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                      >  1:906956/1‑137 (MQ=255)
                        gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCCAACATGCCGt                                                                                                                 >  2:972899/1‑140 (MQ=255)
                                                         aaacctgctgACATACGGCGTTTCTTTCATGGTACGGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                                <  2:878650/136‑1 (MQ=255)
                                                                              ttCTTTTATGGTACTGGCGGCTGTCAGCCACACCGCGGCCCGCATTGCGCCCGCGAGAAATACCCCCCCACCGTCAAACATGCCGCTAATCCAGCCGCCCCGCGAACCGCCAATCGCCACCCCAAGAAAACGGCtggggg                                                           >  2:262722/1‑140 (MQ=255)
                                                                                 tttCATGGTACTGGCGCCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGt                                                        <  1:972899/140‑1 (MQ=255)
                                                                                  ttCATGGTACTGGCGACTGTCGGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAT‑CCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaa                                                      <  1:963975/140‑1 (MQ=255)
                                                                                     aTGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGGCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGGAAac                                                    >  1:1002459/1‑140 (MQ=255)
                                                                                          aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                               >  2:1196627/1‑140 (MQ=255)
                                                                                           cgggCGACTGTCAGCCACACGGCGCCCAGCATGGCCCCAGGGAAAAATCCCCCCTGACCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                              <  2:143239/138‑1 (MQ=255)
                                                                                                gACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcccctgtctc                                         >  1:500465/1‑133 (MQ=255)
                                                                                                 aCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                        >  1:143239/1‑140 (MQ=255)
                                                                                                 aCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCTTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                        <  2:13010/140‑1 (MQ=255)
                                                                                                  cTGTCAGCCACCCTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  1:364167/140‑1 (MQ=255)
                                                                                                  cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  1:730744/140‑1 (MQ=255)
                                                                                                   tGTCAGCCATACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:1142508/1‑140 (MQ=255)
                                                                                                   tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:822491/1‑140 (MQ=255)
                                                                                                   tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                      >  1:1008305/1‑140 (MQ=255)
                                                                                                        gCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCtttt                                 >  2:352967/1‑139 (MQ=255)
                                                                                                            cagggCGGCCAGTTTCGCGCCAGCGAGAAATCCCCCCTCACCGCAAAACAGGCCGTTAACCCAGCCACCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAc                             <  2:1002459/136‑1 (MQ=255)
                                                                                                            cacTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAc                             >  2:1177242/1‑140 (MQ=255)
                                                                                                             acTGCGGCCAGCATTGCGCCAGCGAGAAATACCCACTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                            >  2:620522/1‑140 (MQ=255)
                                                                                                               tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  1:1305225/1‑140 (MQ=255)
                                                                                                               tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  1:221841/1‑140 (MQ=255)
                                                                                                                 cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  1:209766/1‑140 (MQ=255)
                                                                                                                 cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  1:299789/1‑140 (MQ=255)
                                                                                                                   gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      >  1:384913/1‑140 (MQ=255)
                                                                                                                      aGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  1:544314/140‑1 (MQ=255)
                                                                                                                      aGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   >  2:243767/1‑140 (MQ=255)
                                                                                                                      aGCATTGCCCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  1:262722/140‑1 (MQ=255)
                                                                                                                       gctttGCGCCAGGGAAAAATACCCCCTGCCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAc                  <  1:1316555/137‑1 (MQ=255)
                                                                                                                        cATTGCGCCGGCGAAAAACCCCCCCGGACCGACAAACATGCCGTTAATCCAGCCCCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  2:209766/140‑1 (MQ=255)
                                                                                                                        cATTGCGCCAGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  1:1214990/140‑1 (MQ=255)
                                                                                                                          ttGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct               <  1:1182262/140‑1 (MQ=255)
                                                                                                                          ttGCCCCAGGGAAAAATACCCCCTGCCCGTCAAACATGCCTTTAATCCAGCCGCCCAGCGACCCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAATAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct               <  2:1256429/140‑1 (MQ=255)
                                                                                                                            gcgcCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  1:1177242/140‑1 (MQ=255)
                                                                                                                            gcgcCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  1:1098701/140‑1 (MQ=255)
                                                                                                                            cgaccACCGAAAAACCCCCCCCGCCCATCAACCTTGCCTTTAATCCAGCCCCCCACCGATCCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  2:1305225/137‑1 (MQ=255)
                                                                                                                            agacaggCGAAAAATCCCCCCGGCCCGTAAAACAGGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt             <  2:51892/134‑1 (MQ=255)
                                                                                                                             cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa            >  1:1384249/1‑140 (MQ=255)
                                                                                                                              cccggCGAGAAATCCCCCCGGACCGTAAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTAc           <  2:822491/136‑1 (MQ=255)
                                                                                                                                 aGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa        <  2:299789/140‑1 (MQ=255)
                                                                                                                                  gggaaaaaTCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt       <  2:221841/135‑1 (MQ=255)
                                                                                                                                       aaaTCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTg  <  2:548072/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445419‑445693

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: ‑