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

ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445425‑445693
                                                                                                                                     |                                                                                                                                       
aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                   <  2:614550/140‑1 (MQ=255)
                tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTTCGCCCGCGAGAAATACCCCCTGACCGTCAAACATGcc                                                                                                                   >  2:1010282/1‑140 (MQ=255)
                          cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                         <  2:713752/140‑1 (MQ=255)
                                          ttCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACCGCGGCCCGCACTTCGCCCGCGAGAAAAAACCCCCCACCCTCAAAAAAGCCGTTAATCCAACCCCCCCGCGAACCGc                                                                                         >  2:290103/1‑140 (MQ=255)
                                                 cgcAAACTGCTGACATACGGCGGTTTTTTTTTGGTACTGGCGGCTGTCAGCCCCACTGCGGCCCGCATTTCGCCCGCGAGAAAAACCCCCCCACCCTCCAACATACCGGTTATCCAGCCCCCCACCGAACcccccaaccc                                                                                  >  2:581374/1‑134 (MQ=255)
                                                   aaacctgctgCCATCGGGCGTTTTTTTCATGGTCCTGGCGCCTGTCAGCCCCCCTGCGGCCAGCATTGCCCCGGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                                <  2:266217/136‑1 (MQ=255)
                                                                  cggcggTTCTTTTATGGTACTGGCGACTGTCAGCCACACCGCGGCCAGCATTGCGCCAGCGAGGAATACCCCCCCACCCTCAAAAATGCCGTTAATCCAGCCCCCCCCCGAACCGCCCATCCCCCCCCCCAGAAAatggc                                                                 >  2:191405/1‑140 (MQ=255)
                                                                        ttCTTTTATGGTACTGGCGACTGTCAGCCACACCGCGGCCCGCATTTCGCCAGCGAGAAAAAACCCCCGACCGTCAAAAATGCCGGTTATTCAGCCGCCCCGCGAAACGCCCATCGCCCCGCCAAAAAAAAGGCtggggg                                                           >  2:619117/1‑140 (MQ=255)
                                                                        ttCTTTCATGGTCCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtgg                                                           <  1:619117/140‑1 (MQ=255)
                                                                              cATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                     <  1:4909/140‑1 (MQ=255)
                                                                                   tACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccc                                                <  1:179617/140‑1 (MQ=255)
                                                                                    aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                               >  1:355273/1‑140 (MQ=255)
                                                                                      tGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                             >  2:613484/1‑140 (MQ=255)
                                                                                          gcttGTCAGCCACCCTGGGGCCAGCATTGCCCCGGCGAAAAATCCCCCCTGCCGGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                         <  2:702098/137‑1 (MQ=255)
                                                                                            cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                       <  1:1010282/140‑1 (MQ=255)
                                                                                            cTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGACTAAACACCCATCGCCGTa                                       >  1:815977/1‑140 (MQ=255)
                                                                                                      cacTGCGGCCGGTATTGCGCCAGCGAGAAATGCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCCCCCAGCAAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAc                             <  2:1058441/140‑1 (MQ=255)
                                                                                                         tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  2:445382/1‑140 (MQ=255)
                                                                                                         tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                          >  1:477575/1‑140 (MQ=255)
                                                                                                           cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                        >  2:430622/1‑140 (MQ=255)
                                                                                                             gCCAGCTTGGCGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      <  2:948641/140‑1 (MQ=255)
                                                                                                             gCCAGCATTGCGCCAGCGAGAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      <  1:518404/140‑1 (MQ=255)
                                                                                                             cccAGCTTTGCCCCAGCGAAAAATCCCCCCGGACCGCAAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                      <  2:477575/139‑1 (MQ=255)
                                                                                                                aGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                   <  2:355273/140‑1 (MQ=255)
                                                                                                                  cATTGCGCCAGCGAGAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  1:191405/140‑1 (MQ=255)
                                                                                                                  cATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  1:477885/140‑1 (MQ=255)
                                                                                                                  cATTGCGCCAGCGAGAAATCCCCCCGGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATTGCCACGCCAAGAAACTGGCTGGTGGACTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 <  2:815977/140‑1 (MQ=255)
                                                                                                                  cATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                 >  1:518804/1‑140 (MQ=255)
                                                                                                                       cgcCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa            <  1:613484/140‑1 (MQ=255)
                                                                                                                       cgcCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCAGGTGGAGTAAACCCCCATCGCCGTACCTTTGTAACCTGCTgg                      <  1:618889/130‑1 (MQ=255)
                                                                                                                       cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCAGGTGGAGTAAACCCCCATCGCCGTACCTTTGTAACCTGCTgg                      >  2:618889/1‑130 (MQ=255)
                                                                                                                           aGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa        <  1:187264/140‑1 (MQ=255)
                                                                                                                           aGCGAAAAATACCCCCGGCCCTCCAAACATGCCGTTAATCCACCCGCCAAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa        <  2:621662/140‑1 (MQ=255)
                                                                                                                                 aaaTACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTg  >  1:621662/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                       
ACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTG  >  NC_000913/445425‑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: ‑