| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 702,809 | G→C | 7.9% | D612E (GAC→GAG) | ymcA ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 702,809 | 0 | G | C | 7.9% | 183.6 / 5.1 | 63 | D612E (GAC→GAG) | ymcA | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (31/27); new base C (2/3); total (33/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 | |||||||||||
TAGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGC > minE/702745‑702878 | tAGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTc < 1:1132283/67‑1 (MQ=255) aGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATc < 1:1937831/70‑1 (MQ=255) aGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCg > 1:3990515/1‑71 (MQ=255) aCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGt < 1:3184954/61‑1 (MQ=255) aCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGt < 1:3173059/61‑1 (MQ=255) aCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGt < 1:1915187/61‑1 (MQ=255) aCGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGt < 1:2120484/61‑1 (MQ=255) aCTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGcc > 1:4071958/1‑69 (MQ=255) aCTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCtt > 1:1781640/1‑71 (MQ=255) aCTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCt < 1:5609738/70‑1 (MQ=255) aCTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCt < 1:5111406/70‑1 (MQ=255) cTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCtt > 1:3056408/1‑70 (MQ=255) aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATc < 1:5502744/71‑1 (MQ=255) aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTAAAAGCGTTTGGCGCCg > 1:4320808/1‑51 (MQ=39) ccgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCttt < 1:3106262/70‑1 (MQ=255) cgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATc > 1:2368684/1‑52 (MQ=255) gccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGc > 1:5199373/1‑71 (MQ=255) gccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGc > 1:1568474/1‑71 (MQ=255) gccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGc > 1:3342361/1‑71 (MQ=255) cATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTg > 1:5674165/1‑67 (MQ=255) tACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTAt < 1:2749373/59‑1 (MQ=255) aCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGAGCGTTGCGCCTTTATCTTTTGCCAGa > 1:1238703/1‑69 (MQ=255) ccccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtttt > 1:1812543/1‑62 (MQ=255) ccccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGatat > 1:2570779/1‑71 (MQ=255) ccccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGatat > 1:3008816/1‑71 (MQ=255) cTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGata > 1:4318178/1‑67 (MQ=255) tACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtttt > 1:3411183/1‑58 (MQ=255) tACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtttt > 1:5010922/1‑58 (MQ=255) accaccACGCCGCTATCAAAGCGTTTGGCGATCTCCAGCGTGCCGCCtt > 1:891513/1‑49 (MQ=12) accaccACGCCGCTATCAAAGCGTTTGGCGATCTCCAGCGTGCCGCCTTTATCCCCTGCCAGATACTGccc < 1:298773/71‑1 (MQ=255) ccaccACGCCGCTATCAAAGCGTTTGGCGATCTCCAGCGTGCCGCCTTTATCCCCTGCCAGa < 1:5241443/59‑1 (MQ=255) ccaccACGCCGCTATCAAAGCGTTTGGCGATCTCCAGCGTGCCGCCTTTATCCCCTGCCAGa < 1:326507/59‑1 (MQ=255) ccaccACGCCGCTATCAAAGCGTTTGGCGATCTCCAGCGTGCCGCCTTTATCCCCTGCCAGATACTGccc > 1:5143422/4‑70 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCtt > 1:3179493/1‑48 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGc > 1:5126942/1‑58 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGccc > 1:1884513/1‑70 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGccc > 1:3446746/1‑70 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGccc > 1:561246/1‑70 (MQ=255) ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGccc < 1:881730/70‑1 (MQ=255) cGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTa < 1:601751/49‑1 (MQ=255) ccACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCaa < 1:288315/69‑1 (MQ=255) cACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGc < 1:1653045/71‑1 (MQ=255) cACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGc < 1:459706/71‑1 (MQ=255) cACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGc < 1:90810/71‑1 (MQ=255) aCGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGa < 1:3561650/57‑1 (MQ=255) cgccgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGa < 1:1531551/56‑1 (MQ=255) cgccgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTa < 1:528289/71‑1 (MQ=255) gccgcTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTa < 1:521562/70‑1 (MQ=255) gTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtca > 1:1289911/1‑71 (MQ=255) tCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtca < 1:2857177/70‑1 (MQ=255) aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtt > 1:4887726/1‑38 (MQ=255) aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGGCAGATATTGCCCAACGCTAAGTTtcat > 1:418394/1‑69 (MQ=255) aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCaa > 1:2892223/1‑56 (MQ=255) aaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAg < 1:1924502/71‑1 (MQ=255) aGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtt > 1:4304310/1‑36 (MQ=255) aGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCtt > 1:2909569/1‑36 (MQ=255) cGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGa > 1:2043584/1‑42 (MQ=255) cGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCAAGa > 1:1337801/1‑42 (MQ=255) tttGGCGACGTCGATCGTTGCGCCTTTATCTTTTGcc < 1:2275777/37‑1 (MQ=255) gACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTg < 1:2907826/70‑1 (MQ=255) cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGc > 1:1998746/1‑71 (MQ=255) cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGc > 1:2266784/1‑71 (MQ=255) cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGc < 1:4467349/71‑1 (MQ=255) cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGc < 1:730130/71‑1 (MQ=255) gTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTtca > 1:2535448/1‑53 (MQ=255) | TAGATACGTTACTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGC > minE/702745‑702878 |
| 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: ‑ |