Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 704,335 | A→T | 6.3% | S104T (TCA→ACA) ‡ | ymcA ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 704,335 | 0 | A | T | 6.3% | 177.8 / 3.6 | 63 | S104T (TCA→ACA) ‡ | ymcA | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (30/29); new base T (2/2); total (32/31) | |||||||||||
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.71e-01 |
ACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG > minE/704268‑704399 | aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg < 1:298368/71‑1 (MQ=255) cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa < 1:4845245/71‑1 (MQ=255) tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc > 1:3675203/1‑71 (MQ=255) tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc > 1:1750856/1‑71 (MQ=255) gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg > 1:2796996/1‑71 (MQ=255) cGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg < 1:2704464/70‑1 (MQ=255) cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc < 1:4096640/71‑1 (MQ=255) cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc < 1:3651759/71‑1 (MQ=255) aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa < 1:2746868/71‑1 (MQ=255) aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa < 1:2124752/71‑1 (MQ=255) gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt < 1:3557014/71‑1 (MQ=255) gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCAGAAAt < 1:2269246/71‑1 (MQ=255) tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa > 1:1528892/1‑68 (MQ=255) tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa > 1:130526/1‑68 (MQ=255) aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt > 1:3500166/1‑54 (MQ=255) aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt > 1:1236879/1‑54 (MQ=255) aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt > 1:1163262/1‑71 (MQ=255) aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt > 1:1807894/1‑71 (MQ=255) aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt > 1:3100519/1‑71 (MQ=255) aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt > 1:4132175/1‑70 (MQ=255) cttcttCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc > 1:4450791/1‑71 (MQ=255) ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt > 1:3540141/1‑65 (MQ=255) ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCAGCTGAAATCTTcc > 1:2706530/1‑70 (MQ=255) cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc > 1:4357608/1‑71 (MQ=255) cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc > 1:3531913/1‑71 (MQ=255) acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa < 1:4293144/55‑1 (MQ=255) acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCta > 1:428365/1‑71 (MQ=255) caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGGTGAAATCTTCCCACTGGc < 1:4277656/68‑1 (MQ=255) caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt > 1:875528/1‑56 (MQ=255) caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCtat < 1:1734292/71‑1 (MQ=255) cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt > 1:1814948/1‑47 (MQ=255) cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt > 1:2227201/1‑47 (MQ=255) gAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa > 1:3783593/1‑38 (MQ=255) aaGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt < 1:5651640/50‑1 (MQ=255) aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc > 1:4380159/1‑57 (MQ=255) aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg < 1:3769852/71‑1 (MQ=255) aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg < 1:5130217/71‑1 (MQ=255) gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg < 1:4611536/44‑1 (MQ=255) gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg < 1:4574332/44‑1 (MQ=255) tttAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGc < 1:3889830/68‑1 (MQ=255) tAAAATCGAATGATTTGTCTTTGTATGACTGATcgcc > 1:1557522/1‑37 (MQ=255) tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCa > 1:2202457/1‑52 (MQ=255) tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCgtgt < 1:2930932/70‑1 (MQ=255) aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc > 1:3072299/1‑48 (MQ=255) aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTgcg < 1:4934951/71‑1 (MQ=255) tCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcac < 1:3229673/71‑1 (MQ=255) gAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcaca < 1:4553275/70‑1 (MQ=255) aTGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt > 1:2240798/1‑38 (MQ=255) tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg > 1:5656822/1‑71 (MQ=255) tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg > 1:3307530/1‑71 (MQ=255) gATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg > 1:4534918/1‑70 (MQ=255) aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta > 1:877800/1‑71 (MQ=255) aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta > 1:3034772/1‑71 (MQ=255) tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat < 1:4540870/71‑1 (MQ=255) tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat < 1:1405437/71‑1 (MQ=255) tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat < 1:3776153/71‑1 (MQ=255) tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat < 1:1142498/71‑1 (MQ=255) tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg < 1:3424822/71‑1 (MQ=255) tCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACg < 1:3245202/71‑1 (MQ=255) tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGaa < 1:4487875/69‑1 (MQ=255) tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGa < 1:5520817/68‑1 (MQ=255) tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGa < 1:843602/68‑1 (MQ=255) tGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATg > 1:184676/1‑71 (MQ=255) aCTGATGGCCGCTGAAATCTTCCCACTGGCTATATTTGCg > 1:715/1‑40 (MQ=255) aCTGATCGCCGCTGAAATCTTCCCACTGGCTATAttt > 1:996110/1‑37 (MQ=255) | ACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTATAACGAATG > minE/704268‑704399 |
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: ‑ |