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

AAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGA  >  NC_000913/445423‑445694
                                                                                                                                       |                                                                                                                                        
aaaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATAcccc                                                                                                                                      <  1:872098/140‑1 (MQ=255)
  aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                    <  2:34381/140‑1 (MQ=255)
  aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                    <  1:802358/140‑1 (MQ=255)
  aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                    >  2:170053/1‑140 (MQ=255)
        ttttAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTTCGCCCGCGAGGAAAACCCCCCGACCgc                                                                                                                              >  2:221095/1‑139 (MQ=255)
         tttAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTc                                                                                                                             <  2:709924/140‑1 (MQ=255)
               gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                       <  1:170053/140‑1 (MQ=255)
               gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                       <  1:333596/140‑1 (MQ=255)
                  tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGcc                                                                                                                    >  2:585883/1‑140 (MQ=255)
                    gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGt                                                                                                                  >  1:779709/1‑140 (MQ=255)
                              aTGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAgcc                                                                                                        <  2:302911/140‑1 (MQ=255)
                                           tttCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGAATTGCGACCGTCAACCACACTCCGGCCCGCATTGCGCCCGCGAGAAAAAACCCCTCACCCTCAAAAATGCCGGTAAACCAGCCCCCCAGCGaaacc                                                                                           >  2:403922/1‑137 (MQ=255)
                                                                    ggggggTTTTTTCAGGGTACGGGCGACTGTCAGCCCCCCTGCGGCCAGCATTGCGCCGGGGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                  <  2:62404/136‑1 (MQ=255)
                                                                    cggcggTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                  <  1:403922/140‑1 (MQ=255)
                                                                         gTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                             <  1:585883/140‑1 (MQ=255)
                                                                          ttttttCATGGGACTGGCGACCGTCCCCCACACTGCGGCCAGCATTTCGCCAGCGAGAAAAAACCCCCGACCGCCAAACATGCCGCTAATCCACCCGCCCCGCGAACCGCCCATCGCCCCGCCAAAAAAAAGGCtggggg                                                            >  2:542559/4‑140 (MQ=255)
                                                                           tCTTTCATGGGACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAAACCCCCCGACCGTCAAAAATGCCGTTAATCCACCCGCCCCGCGAAACGCCAAACGCCCCGCCAAGAAACCGGCtgggggg                                                           >  2:405805/1‑139 (MQ=255)
                                                                              ttCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTa                                                        <  1:795807/140‑1 (MQ=255)
                                                                               tCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCGGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaa                                                       <  1:405805/140‑1 (MQ=255)
                                                                                cATGGTCCTGGCGCCTGCCAGCCACACTGGGGCCAGCATGGCCCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                      <  2:779709/140‑1 (MQ=255)
                                                                                cATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                      <  1:221095/140‑1 (MQ=255)
                                                                                    gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                  >  2:199430/1‑140 (MQ=255)
                                                                                     tACTGGCGCCTGCCGGCCACCCTGGGGCCAGCATGGCGCCGGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccc                                                 <  2:585441/140‑1 (MQ=255)
                                                                                       cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                               >  2:393855/1‑140 (MQ=255)
                                                                                            cacggTCACCCACCCGGGGGCAAGCTTTGCGCCAGCGAAAAATCCCCCCTGCCCTTAAAACATGCCGTTAATCCAGCCGCCCAGCGAGCCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                          <  2:559924/136‑1 (MQ=255)
                                                                                             aCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                         >  2:267839/1‑140 (MQ=255)
                                                                                              cTGTCAGCCCCCCTGCGGCCAGGTTGGCGCCGGCGAAAAATCCCCCATGCCGTTAAAACATGCCGTAAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                        <  2:128271/140‑1 (MQ=255)
                                                                                              cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                        <  1:875072/140‑1 (MQ=255)
                                                                                               tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                       >  1:723821/1‑140 (MQ=255)
                                                                                                    gCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                                  >  1:102491/1‑140 (MQ=255)
                                                                                                     ccccccTGCGGCCAGCATTGCCCCAGGGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCACCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGt                                 <  1:542559/135‑1 (MQ=255)
                                                                                                         acTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                             >  2:90007/1‑140 (MQ=255)
                                                                                                           tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                           >  2:221866/1‑140 (MQ=255)
                                                                                                           tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                           >  1:865633/1‑140 (MQ=255)
                                                                                                           tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                           >  1:11260/1‑140 (MQ=255)
                                                                                                               gCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                         <  1:875892/122‑1 (MQ=255)
                                                                                                               gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                         >  2:875892/1‑122 (MQ=255)
                                                                                                                         cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCAAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa             >  1:314182/1‑140 (MQ=255)
                                                                                                                         cgcCAGCGAAAAATCCCCCCGGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa             <  2:102491/140‑1 (MQ=255)
                                                                                                                          cccAGCGAGAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTAc            <  2:865633/139‑1 (MQ=255)
                                                                                                                              gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    <  1:817778/128‑1 (MQ=255)
                                                                                                                              gCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    >  2:817778/1‑128 (MQ=255)
                                                                                                                                  aaaaTCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt    <  1:199430/139‑1 (MQ=255)
                                                                                                                                    aaTCCCCCCTGCCCTTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGa  <  2:737474/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
AAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGA  >  NC_000913/445423‑445694

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