Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,167,337 | A→G | K18E (AAG→GAG) | secE → | preprotein translocase subunit SecE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,167,337 | 0 | A | G | 94.7% | 57.1 / ‑3.2 | 19 | K18E (AAG→GAG) | secE | preprotein translocase subunit SecE |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (6/12); total (6/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01 |
AAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCG > NZ_CP009273/4167203‑4167462 | aaaaGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTg < 1:515470/139‑1 (MQ=255) cTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCg < 1:16829/139‑1 (MQ=255) tttCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGc < 2:562051/139‑1 (MQ=255) tttGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATcg < 1:326286/139‑1 (MQ=255) tttGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATcg < 2:524219/139‑1 (MQ=255) aCAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCg < 2:207178/138‑1 (MQ=255) ggttTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCtgcgtgc > 1:157128/1‑139 (MQ=255) cGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCgtag < 2:354635/138‑1 (MQ=255) aTACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAAt < 2:44028/139‑1 (MQ=255) gAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGa > 1:357450/1‑139 (MQ=255) gAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGa > 1:559808/1‑139 (MQ=255) cTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATtgc > 1:83835/1‑139 (MQ=255) aGGAAGCGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCa < 1:264447/139‑1 (MQ=255) nngngCGNGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCa < 2:83835/132‑1 (MQ=255) cGGGCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGggt > 2:324963/1‑139 (MQ=255) gggCGCGGCCTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGggtg > 1:367148/1‑139 (MQ=255) ccTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGTCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTcgcg < 2:422072/139‑1 (MQ=255) ccTGGAAGCGATGGAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTcgcg < 1:5836/139‑1 (MQ=255) ccTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTcgcg < 1:381090/139‑1 (MQ=255) | AAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGATCGCGGGGTGAAAATGTTTGTAGAAAACTTCTGACAGGTTGGTTTATGAGTGCGAATACCGAAGCTCAAGGAAGCGGGCGCGGCCTGGAAGCGATGAAGTGGGTCGTTGTGGTGGCATTGCTCCTGGTGGCGATTGTCGGCAACTATCTTTATCGCGACATTATGCTGCCGCTGCGTGCGCTGGCCGTAGTAATTCTGATTGCTGCAGCGGGTGGTGTCGCG > NZ_CP009273/4167203‑4167462 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |