Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,353,663 G→T A251E (GCA→GAA)  yihF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,353,6630GT95.4% 197.0 / 2.4 65A251E (GCA→GAA) yihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (62/0);  total (65/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGTAGGTGAGATTGACATTGATATAATCTTTGTCTTTAACGCGATCCAGCCTGGTTTTTGCAACGAAAGCATCGATTTGGGCAACGTCTTCCCC  >  minE/2353604‑2353697
                                                           |                                  
cGTAGGTGAGATTGACATTGATATAATCTTTGTCTTTAACGCGATCCAGCCTGGTTTTTGCa                                  >  1:479672/1‑62 (MQ=255)
cGTAGGTGAGATTGACATTGATATAATCTTTGTCTTTAACGCGATCCAGCCTGGTTTTTGCa                                  >  1:419499/1‑62 (MQ=255)
cGTAGGTGAGATTGACATTGATATAATCTTTGTCTTTAACGCGATCCAGCCTGGTTTTTGCa                                  >  1:514054/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGc                          >  1:41468/1‑38 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATc                       >  1:329596/1‑41 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:729529/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:730583/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:717386/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:699519/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:681654/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:675852/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:579942/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:572050/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:545965/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:520964/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:510979/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:48273/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:452628/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:739661/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:759722/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:762570/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:773932/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:784452/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:814481/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:820795/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:874640/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:891790/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:895530/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:899282/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:915608/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:920820/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:943907/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:994928/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1148267/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:102725/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1032156/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1051238/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1056313/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:110264/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1102687/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1102866/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1111319/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1112605/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1113212/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1115721/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1119665/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1130733/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:452652/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:239398/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1023815/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:371830/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:34478/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:310909/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:284157/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:267403/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:247421/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1169529/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:219785/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:216984/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:205839/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:125488/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1186910/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTcccc  >  1:1183915/1‑62 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTccc   >  1:780073/1‑61 (MQ=255)
                                tCTTTAACGCGATCCAGCCTGGTTTTTTCAACGAAAGCATCGATTTGGGCAACGTCTTaccc  >  1:937316/1‑62 (MQ=255)
                                                           |                                  
CGTAGGTGAGATTGACATTGATATAATCTTTGTCTTTAACGCGATCCAGCCTGGTTTTTGCAACGAAAGCATCGATTTGGGCAACGTCTTCCCC  >  minE/2353604‑2353697

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