Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 400,660 A→G 100% F292F (TTT→TTC dacA ← D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE400,6600AG81.8% 93.2 / 14.8 44F292F (TTT→TTCdacAD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base A (8/0);  new base G (36/0);  total (44/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 = 8.51e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAAACCAAACCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCA  >  minE/400602‑400705
                                                          |                                             
caaacgaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1367519/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGac                                          >  1:2186084/1‑63 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:616286/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1092106/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1967676/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:2066466/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:2120686/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:2246763/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:2296182/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:265693/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:526956/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1857671/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:629246/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:639885/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:709150/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:734652/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:741178/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:863740/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:878347/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:934797/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1098600/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1262691/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1201852/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1520686/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1525658/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1565659/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1569957/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1573333/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1182053/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1175858/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1629763/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:169961/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1705708/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1706003/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCGAAGAAAc                                        >  1:1797674/1‑66 (MQ=255)
caaaccaaacCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGATAACGGTTTCGAAGAAAc                                        >  1:1806191/1‑66 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTTCTTTCGGCTTCa  >  1:572118/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:1186499/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:516669/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:1284985/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:2262741/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:2260064/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:159890/1‑67 (MQ=255)
                                     tttCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCa  >  1:1843315/1‑67 (MQ=255)
                                                          |                                             
CAAACCAAACCGGTTCAGAGGCGAACTCTTTACCTACTTTCAGTGGGTTAACGGTTTCAAAGAAACGGAAGCCCCAGGTTAGCAGTTTTTTACTTTCGGCTTCA  >  minE/400602‑400705

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: