Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 372,491 T→C 53.0% G1146G (GGT→GGC entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE372,4910TC53.0% 0.8 / 24.4 34G1146G (GGT→GGCentFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (3/13);  new base C (7/11);  total (10/24)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCC  >  minE/372424‑372556
                                                                   |                                                                 
ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGt                                                                   <  1:3197962/68‑1 (MQ=255)
                          taCCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTg                                          <  1:1652562/67‑1 (MQ=255)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCg                                            >  1:1352968/1‑61 (MQ=255)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:2272703/1‑69 (MQ=39)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:2465860/1‑69 (MQ=39)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:2579821/1‑69 (MQ=39)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:1826873/1‑69 (MQ=39)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:3429452/1‑69 (MQ=39)
                              tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc                                    >  1:524889/1‑69 (MQ=39)
                               gctgGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTg                                          <  1:2007268/62‑1 (MQ=255)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:3223445/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:3382497/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:3130289/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:2970658/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:2522928/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:583376/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:92744/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:1521132/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:1915921/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:1812631/68‑1 (MQ=39)
                                     ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc                             <  1:1796389/68‑1 (MQ=39)
                                                    tACGCTGGCGCAGGGTATTGCTGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGc               >  1:889250/1‑68 (MQ=255)
                                                    tACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGc               >  1:1628692/1‑68 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:1035685/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:2247538/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:843068/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:1189004/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:563706/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:1783870/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:3277477/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:3168836/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:2963004/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             >  1:2903532/1‑69 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:1945761/69‑1 (MQ=255)
                                                     aCGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:2066157/69‑1 (MQ=255)
                                                                   tATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTg             <  1:2786992/55‑1 (MQ=255)
                                                                   tATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGgcc  >  1:1792650/1‑66 (MQ=255)
                                                                   tATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGgcc  >  1:2818635/1‑66 (MQ=255)
                                                                   |                                                                 
TTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCC  >  minE/372424‑372556

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: