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

TAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGA  >  NC_000913/445419‑445687
                                                                                                                                             |                                                                                                                                 
tAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATa                                                                                                                                   <  1:655742/140‑1 (MQ=255)
    aaaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCCCACCGCGGCCAGCATT‑TCGCCAGCGAGAAATAcccc                                                                                                                               >  2:634490/1‑140 (MQ=255)
      aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTG‑TCAGCCCCACCGCGGCCCGCATT‑GCGCCAGCGAGAAAAAcccccc                                                                                                                             >  2:655742/1‑139 (MQ=255)
       cGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTg                                                                                                                            <  1:532870/140‑1 (MQ=255)
                  cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                 <  1:70799/140‑1 (MQ=255)
                   gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                <  1:634490/140‑1 (MQ=255)
                   gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                <  2:372871/140‑1 (MQ=255)
                   gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                <  2:779750/140‑1 (MQ=255)
                      tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTCACCCTCAAACATGcc                                                                                                             >  1:649826/1‑140 (MQ=255)
                        gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTTTTTTATGGTACTGGCGGCTG‑TTAGCCACACTTCGGCCAGCATT‑TCGCCAGCGGGAAATACCCCCCGACCGGCAAACAAGCCGt                                                                                                           >  2:840669/1‑140 (MQ=255)
                                          cGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGGACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAAAACCCCCCGACCGTCAAACATTCCGTTAATCCAGCCCCCCcgcg                                                                                         >  2:91505/1‑140 (MQ=255)
                                                         cAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                          <  1:291079/140‑1 (MQ=255)
                                                         aaacctgctgACATACGGGGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCTTT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                          <  2:868342/136‑1 (MQ=255)
                                                                       ccggcggTTCTTTCGTGGTACCGGGCACTGATCA‑CCCCCCTGGGGCTAGCTGG‑GAGCCGGGAAAAAATACCCCCCGACCTCCAAACATGCCGTTAATCCAGCCGCCCAGGGAACCGCCAATCGCCACGCCAAGAAActgg                                                            <  2:363456/139‑1 (MQ=255)
                                                                       aCGGCGGTTCTTTTATGGTACTGGCGACTG‑TCAGCCCCACTGCGGCCAGCATT‑TCGCCAGCGAGAAAAACCCCCCGACCGCCAAACATGCCGTTAATCCAGCCGCCCCGCGAACCCCCAAACGCCACGCCCAAAAAatgg                                                            >  2:803403/1‑140 (MQ=255)
                                                                        gggcggTTTTTTCAGGGTACTGGCGACTG‑TCAGCCACACGGCGGCCAGCATT‑GCGCCAGCGAAAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                           <  2:649826/139‑1 (MQ=255)
                                                                          gcggTTCTTTCATGGTACTGGCGCCTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggctg                                                         <  1:713624/140‑1 (MQ=255)
                                                                                   tCATGGTACGGGCGACTG‑TCACCCACCCGGCGGCCAGCATT‑GCGCCGGCGAGAAATCCCCCCTGACCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaa                                                <  2:42459/140‑1 (MQ=255)
                                                                                         tccgggCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccc                                          <  1:162676/136‑1 (MQ=255)
                                                                                           cTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                        >  1:891160/1‑140 (MQ=255)
                                                                                           cTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                        >  2:454971/1‑140 (MQ=255)
                                                                                                gccTG‑TCAGCCACCCTGCGGCCAGCATT‑GCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                   <  1:91505/138‑1 (MQ=255)
                                                                                                gACTG‑CCAGCAACGCGGCGGCC‑GCATTGGCGCAAGCGAAAAATACCCCCTTACCCTAAAACATGCTGTTAATCCAGCCGCCAAAAAAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACGCCCATCGCCg                                   <  2:924270/140‑1 (MQ=255)
                                                                                                 aCTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                  >  1:924270/1‑140 (MQ=255)
                                                                                                 aCTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                  <  1:112825/140‑1 (MQ=255)
                                                                                                  cTG‑TCGGCCACACTGCGGCAAGCATG‑GCGCCAGCGAAAAATCCCCCCTCCCCGTAAAACATGCCGTAAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                 <  2:145269/140‑1 (MQ=255)
                                                                                                  cTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                 <  2:857945/140‑1 (MQ=255)
                                                                                                   tG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                   >  2:354709/1‑137 (MQ=255)
                                                                                                   tG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                   <  1:354709/137‑1 (MQ=255)
                                                                                                   tG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                >  2:194841/1‑140 (MQ=255)
                                                                                                   tG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                >  1:921890/1‑140 (MQ=255)
                                                                                                         gCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                           >  2:535303/1‑140 (MQ=255)
                                                                                                          ccACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCtttat                          >  1:857945/1‑138 (MQ=255)
                                                                                                               cTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                 >  2:124491/1‑128 (MQ=255)
                                                                                                               cTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                 <  1:124491/128‑1 (MQ=255)
                                                                                                                tGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                    >  1:457802/1‑140 (MQ=255)
                                                                                                                       aGCATT‑GCCCCAGCGAAAAATCCCCCCTGCCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa             <  2:891160/140‑1 (MQ=255)
                                                                                                                         cttt‑GCCCCAGGGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct           <  1:205633/138‑1 (MQ=255)
                                                                                                                               cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa      >  1:352946/1‑140 (MQ=255)
                                                                                                                                   aGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa  <  1:881703/140‑1 (MQ=255)
                                                                                                                                    gCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa             <  1:544998/128‑1 (MQ=255)
                                                                                                                                    gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa             >  2:544998/1‑128 (MQ=255)
                                                                                                                                             |                                                                                                                                 
TAGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTG‑TCAGCCACACTGCGGCCAGCATT‑GCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGA  >  NC_000913/445419‑445687

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