Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,105,255 G→A 100% A123A (GCG→GCA ydiO → predicted acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,105,2550GA84.4% 199.4 / 26.0 83A123A (GCG→GCAydiOpredicted acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base A (0/70);  total (8/75)
Fisher's exact test for biased strand distribution p-value = 3.26e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTAAA  >  minE/1105198‑1105315
                                                         |                                                            
gcTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGaa                                                     <  1:1699184/67‑1 (MQ=255)
gcTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGaa                                                     <  1:1055082/67‑1 (MQ=255)
gcTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGaa                                                     <  1:863861/67‑1 (MQ=255)
gcTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGaa                                                     <  1:660877/67‑1 (MQ=255)
 cTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGaa                                                     <  1:1070906/66‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACCACAGTGCCACTACCACTTacac         <  1:402949/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACTGTGCCACTACCACTTacac         <  1:2302506/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:228796/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:333839/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:315281/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:30566/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:274860/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:262491/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:254073/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:252523/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:398066/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2281140/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2213889/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2210961/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2188863/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:875799/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2142947/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:639907/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:778482/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:736380/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:700160/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:694273/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:684107/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:681855/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:661761/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:657060/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:337042/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:603191/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:579566/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:779844/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:531874/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:518016/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:40098/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2121630/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:356258/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1225514/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:970973/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1398327/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1364300/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1342175/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1305905/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:993394/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1261065/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1241358/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1418866/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1218300/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1176890/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:115211/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1130590/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1128019/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1107741/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1062707/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:100318/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:14785/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2089905/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2075716/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:20490/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2046042/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2043698/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1982955/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1968156/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1929684/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1907909/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1907092/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1903261/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1749141/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1670280/67‑1 (MQ=255)
                                            ccTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:2135161/67‑1 (MQ=255)
                                             ctttGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:957055/63‑1 (MQ=255)
                                             cTATGCCCTGGCATTGACGGAGCCAGGCGCAGGCTCAGACAACAACAGTGCCACTACCACTTacac         <  1:1830835/66‑1 (MQ=255)
                                                   ccTGGCGTTTACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:1013605/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:566971/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:794295/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:2182703/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:2151351/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:2078905/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:1403226/1‑67 (MQ=255)
                                                   ccTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTaaa  >  1:1302214/1‑67 (MQ=255)
                                                         |                                                            
GCTACGTAAAACGGCAGAGAGCACCCTGGAAACGGGTGATCCCGCCTATGCCCTGGCGTTGACGGAACCAGGTGCTGGCTCAGATAACAACAGTGCCACTACCACTTACACGCGTAAA  >  minE/1105198‑1105315

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: