Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 997,102 G→T 100% A61E (GCA→GAA)  ynfL ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE997,1020GT87.3% 95.2 / 7.0 47A61E (GCA→GAA) ynfLpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (41/0);  total (42/5)
Fisher's exact test for biased strand distribution p-value = 3.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCGTCATCCACCATAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCC  >  minE/997044‑997134
                                                          |                                
aGCGTCATCCAACATAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCaa                        <  1:2335957/69‑1 (MQ=255)
         caccaTAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCg               <  1:1400696/69‑1 (MQ=255)
         caccaTAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCg               <  1:475808/69‑1 (MQ=255)
             aTAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTg           <  1:3278007/69‑1 (MQ=255)
                gACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCaa                        <  1:1255813/53‑1 (MQ=255)
                  cAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCg      >  1:1494005/1‑69 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGt     >  1:81095/1‑49 (MQ=39)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:694196/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3308746/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3320991/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3390576/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3390791/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3394090/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3415818/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:358180/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:475383/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:640417/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2998706/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:753367/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:77919/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:81209/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:854096/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:954233/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:968974/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:969102/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3195406/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1210424/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1232353/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1531776/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1549594/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1654333/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2139250/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2192229/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2214871/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:223124/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:228037/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2344975/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2431441/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2460049/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2615073/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2754481/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:2911496/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:1103950/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:3213693/1‑52 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTACGTGc   >  1:2653388/1‑51 (MQ=255)
                                       aaGAAACTGTTTTCCTGCTTCCGTCAGCAATACACTGAGATTGGTCCGTGcc  >  1:290937/1‑52 (MQ=255)
                                       aaGAAAATGTTTTCCTGCTTCCGTCAGCAATACACTGCGATTGGTCCGTGcc  >  1:606003/1‑52 (MQ=38)
                                                          |                                
AGCGTCATCCACCATAGACAGGATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCC  >  minE/997044‑997134

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: