Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,954,513 G→A 63.5% P13P (CCG→CCA rimI → acetylase for 30S ribosomal subunit protein S18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,954,5130GA63.5% 34.1 / 37.9 63P13P (CCG→CCArimIacetylase for 30S ribosomal subunit protein S18
Reads supporting (aligned to +/- strand):  ref base G (16/7);  new base A (17/23);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACG  >  minE/2954451‑2954567
                                                              |                                                      
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAgctgct                                                  <  1:3297299/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:971409/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2365717/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2406522/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2078817/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2506051/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1750366/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:257447/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1657533/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2677809/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1649783/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:2705997/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1559800/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:3115654/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1528990/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:418529/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1298758/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:515671/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1250289/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:555147/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:598658/69‑1 (MQ=255)
cgTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCt                                                  <  1:1157511/69‑1 (MQ=255)
     tGGCAACAAATTTTCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGTCt                                                  <  1:332150/64‑1 (MQ=255)
      ggCAACAAATTTGCACATATGAACACGATTTCTTCCCTCTAAACGACTGATTTACCGGCGGCTTACca                                             <  1:2308998/68‑1 (MQ=255)
      ggCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACca                                             <  1:1613751/68‑1 (MQ=255)
      ggCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACca                                             <  1:481949/68‑1 (MQ=255)
            aaaTTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGa                                       >  1:1070444/1‑68 (MQ=255)
            aaaTTTGCACATATGAACACGATTTCTTCCCTCGAAACGAATGATTTACCGGCGGCTTACCACATTGa                                       >  1:1467647/1‑68 (MQ=255)
                     caTATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             <  1:2366870/69‑1 (MQ=255)
                     caTATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             <  1:2372810/69‑1 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGGGcc                             >  1:3063500/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:2290443/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:1002957/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:537672/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:1278550/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:490430/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:432022/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:1336430/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:2201821/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:2180656/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:2482225/1‑67 (MQ=255)
                       tatGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGcc                             >  1:2172391/1‑67 (MQ=255)
                        atGAATACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:487371/1‑69 (MQ=255)
                        atGAACCCGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGcc                             >  1:673597/1‑66 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1654222/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:76451/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:701419/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:624668/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1213617/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:241438/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:398746/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1482652/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:3279365/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:155596/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1932683/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:2838685/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1598173/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGATTTACCAGCGGCTTACCACATTGAACAACGCGCCCAc                          >  1:1687329/1‑69 (MQ=255)
                        atGAACACGATTTCTTCCCTCGAAACGACTGAATTACCAGCGGCATACCCCATAGAACAAcg                                 >  1:2894708/1‑62 (MQ=255)
                         tGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTAGCACATTGAACAACGCGCCCACg                         >  1:2404198/1‑69 (MQ=255)
                         tGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACg                         >  1:224556/1‑69 (MQ=255)
                                                aCGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGTAGTGAAAAAACg  <  1:1232725/69‑1 (MQ=255)
                                                aCGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACg  <  1:1310226/69‑1 (MQ=255)
                                                              |                                                      
CGTTATGGCAACAAATTTGCACATATGAACACGATTTCTTCCCTCGAAACGACTGATTTACCGGCGGCTTACCACATTGAACAACGCGCCCACGCCTTTCCGTGGAGTGAAAAAACG  >  minE/2954451‑2954567

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: