Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 619,840 T→C 56.7% L383L (TTG→CTG)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE619,8400TC56.7% 16.2 / 42.9 60L383L (TTG→CTG) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base T (4/22);  new base C (27/7);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.16e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAG  >  minE/619774‑619897
                                                                  |                                                         
cacaTCAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1177320/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1070879/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:2650201/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:2950149/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:2477733/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:3417303/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:345100/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:401389/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1744861/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1659960/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1616238/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:159545/69‑1 (MQ=255)
cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                         <  1:1297371/69‑1 (MQ=255)
    acaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAcg                                                     >  1:1417243/1‑69 (MQ=255)
    acaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAcg                                                     >  1:270531/1‑69 (MQ=255)
    acaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAcg                                                     >  1:1012005/1‑69 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:2163479/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:2381204/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:2299335/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:1537235/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:3389611/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:684389/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:1294811/69‑1 (MQ=255)
                          tAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                               <  1:1376759/69‑1 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGAcggcgg                                                 >  1:1200298/1‑48 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTTGTTc                                >  1:2606207/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGtt                                 >  1:41443/1‑64 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:3092776/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:3033876/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:555009/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:624546/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:86794/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2797181/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:271563/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:652556/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2699819/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:950731/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2563367/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1135363/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1275725/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1628580/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1705632/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:183880/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1893395/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1929684/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:1941866/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2072717/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2488167/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2522291/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTc                                >  1:2610896/1‑65 (MQ=255)
                             aaTCAAGTCTATCACTGACTTCGGTATCTTCATCGGCCTGACGGCGGCATCGACGGCCTGGTTc                                >  1:2522317/1‑64 (MQ=255)
                                cAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGt                          >  1:2958611/1‑68 (MQ=255)
                                    tCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGAc                     <  1:894805/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:2177191/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:1958406/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:2247552/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:2854269/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:271132/69‑1 (MQ=255)
                                                       tCTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:1138748/69‑1 (MQ=255)
                                                        cTTCATCGGCCTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAg  <  1:849121/68‑1 (MQ=255)
                                                                  |                                                         
CACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAG  >  minE/619774‑619897

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: