Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,948,123 T→C 25.8% T40A (ACC→GCC)  yjjB ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,948,1230TC25.8% 167.3 / 67.0 93T40A (ACC→GCC) yjjBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (26/43);  new base C (20/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGT  >  minE/2948060‑2948174
                                                               |                                                   
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:255853/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1511653/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1535561/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1641932/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1649753/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1678030/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:177446/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1831718/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:191382/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1953312/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1958784/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1985195/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:2048580/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:2110937/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:2119345/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:102045/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:291464/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:346079/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:368554/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:455650/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:557963/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:600747/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:651422/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:761400/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:770232/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:902801/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:907265/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:912230/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:997408/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:100535/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1046010/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1057541/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:108692/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1432029/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1116576/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1217900/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1242494/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1382445/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:124274/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1334731/66‑1 (MQ=255)
ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:1268213/66‑1 (MQ=255)
 cgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:188917/65‑1 (MQ=255)
            ttGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTaa                                                   <  1:2041144/54‑1 (MQ=255)
                        ccGTTATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATAcc                           >  1:295949/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATAcc                           >  1:865778/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATAcc                           >  1:1725247/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:1325415/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:651154/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:353485/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:1490763/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:1208430/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCa                          >  1:2131951/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATgcgc                                   >  1:571211/1‑58 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:1079007/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:1219184/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:997120/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:849110/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:442521/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:148462/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATAcc                           >  1:1527802/1‑66 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1567153/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:715801/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1060118/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:84806/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1931535/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:378141/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1321556/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:981580/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1814554/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:1228554/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:2080341/1‑67 (MQ=255)
                        ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCa                          >  1:2068548/1‑67 (MQ=255)
                                              gATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt     <  1:977938/66‑1 (MQ=255)
                                              gATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt     <  1:1848953/66‑1 (MQ=255)
                                              gATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt     <  1:1465097/66‑1 (MQ=255)
                                               aTAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt     <  1:2005113/65‑1 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCTGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:112498/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATgg   >  1:1922585/1‑66 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:665563/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1656402/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:584742/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:514782/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:472202/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1152794/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:882745/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1163582/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1170552/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1022188/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:927606/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:330342/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:29510/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1401104/1‑67 (MQ=255)
                                                tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt  >  1:1773977/1‑67 (MQ=255)
                                                               |                                                   
CCGCCGAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGT  >  minE/2948060‑2948174

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