Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,141,978 T→C 27.2% P214P (CCA→CCG nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,141,9780TC27.2% 31.0 / 10.0 33P214P (CCA→CCGnusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base T (9/15);  new base C (7/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 5.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAC  >  minE/2141918‑2142035
                                                            |                                                         
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                                    >  1:2128297/1‑68 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:2994724/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:3393471/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:3219240/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:2963781/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:2584059/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:3180463/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:3145929/1‑69 (MQ=255)
cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   >  1:1839870/1‑69 (MQ=255)
 gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:2036770/68‑1 (MQ=255)
 gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:2200000/68‑1 (MQ=255)
 gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:3058335/68‑1 (MQ=255)
 gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:1774195/68‑1 (MQ=255)
 gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:1578400/68‑1 (MQ=255)
   aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                 <  1:319993/68‑1 (MQ=255)
   aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                 <  1:2332024/68‑1 (MQ=255)
   aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                 <  1:2335523/68‑1 (MQ=255)
   aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                 <  1:1587639/68‑1 (MQ=255)
   aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                 <  1:1065023/68‑1 (MQ=255)
       gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACgg                                            <  1:304069/69‑1 (MQ=255)
              gATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc                                                   <  1:2090222/55‑1 (MQ=255)
                                    ttCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTgg                <  1:541627/68‑1 (MQ=255)
                                    ttCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTgg                <  1:211708/68‑1 (MQ=255)
                                          cACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATGAGCATTTCCGGCTTGGAACGAGt         <  1:1232074/69‑1 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCTCTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:532465/1‑69 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:2664093/1‑69 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:1839685/1‑69 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:1821389/1‑69 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:1425094/1‑69 (MQ=255)
                                          cACTTCTTCGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:391456/1‑69 (MQ=255)
                                          cACTTCTTAGCCGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGt         >  1:162996/1‑69 (MQ=255)
                                                    ccGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGTGACGAAc  <  1:149309/66‑1 (MQ=255)
                                                    ccGATTTCCGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAGCGGGTGACGAAc  <  1:3404718/66‑1 (MQ=255)
                                                            |                                                         
CGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAC  >  minE/2141918‑2142035

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