Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,515,325 | G→A | H197Y (CAC→TAC) | dam ← | DNA adenine methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,515,325 | 0 | G | A | 96.2% | 70.7 / ‑4.1 | 26 | H197Y (CAC→TAC) | dam | DNA adenine methyltransferase |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (13/12); total (14/12) | |||||||||||
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.97e-01 |
GGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACG > NC_000913/3515090‑3515544 | ngTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTa > 1:312653/2‑246 (MQ=255) gNCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAgg > 1:395957/1‑251 (MQ=255) gcagTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCg < 1:153882/213‑1 (MQ=255) agTTCGTCCACCTTTTTACGTGTGCTGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATc < 1:240873/251‑1 (MQ=255) cGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCa > 1:329485/1‑250 (MQ=255) aCCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTggn < 1:95180/212‑2 (MQ=255) gctTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATGCCACGCACGCGTTAACATCGGATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGAACTGCGGCGATCTCCGCCAGATGCGGTTGTTGTTCAAGCGTAAAGCTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTgcgcgcgc < 1:160022/249‑1 (MQ=255) ntCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTgcgccgcctgatcaatatatg > 1:303113/2‑234 (MQ=255) naCTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGcatcata > 1:298936/2‑224 (MQ=255) nnnnnnttnNNNNGCTNATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTa < 2:312653/144‑1 (MQ=255) ncgcgcTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCAATGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCACGCGTACAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCgg > 1:47645/2‑176 (MQ=255) naTACCACTCAAGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCACTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCt > 1:124938/2‑251 (MQ=255) ngTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGt > 1:132063/2‑251 (MQ=255) gATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCACTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTcc > 1:268604/1‑251 (MQ=255) ggAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGACAGATGAGCTTGTTGTTCACGCGTCAAACTGTTTGTGTAATACGCCGTGAAGTTGGCGGTCGCAGACAGc > 1:453904/1‑121 (MQ=255) gaTCAGCACTGGAATATGGCGCTCAACCAGCCCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGAAAAGTTGGCgg > 1:220998/1‑107 (MQ=255) cAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAgg < 1:265432/251‑1 (MQ=255) atNNNNNTCANCCAGCCCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACTCCTTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCGATAGTCGACGGATGGatcatc < 2:298936/144‑1 (MQ=255) cTCAACCAGACCTGCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGGAAAGTTGGCGTTCGCAGACAGCGGTGCATAAGGCGGAGCGCAATAGACGACGGAGTCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTgcgc < 1:434312/192‑1 (MQ=255) gACCTGCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGTGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTGTTTCAGCGAAGTGATACAACTCTGCGTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACg < 1:390362/251‑1 (MQ=255) ctcctcCAGATGCGCGTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTTGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTAGCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTGCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGcc < 1:86078/228‑1 (MQ=255) cnGCCAGNTGCGNTTGTTGTTCAAGCGTNAAACTGTTTGTGTAATACCNNGNAAAGTNGGCgg < 2:220998/63‑1 (MQ=255) naGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCCGTAAAGTTGGc > 1:48574/2‑57 (MQ=255) cAGNTGCGCTTGTTGTTCAAGCGNAAAACTGTTNGTGTANNNCNCCGTNAAGTTGGc < 2:48574/57‑1 (MQ=255) cgcTTGTTGTTCAAGCGTAAAACTGTTTGTGTAATACGCAGTAAAGTTGGCGGTCGCAGAAAGCGGTGCATACGGCGGATCGCAATAGACGACGGCTGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCt > 1:334464/1‑160 (MQ=255) tgttCCAGNGTAANANNNTTTGTNTAATACGCCNTAAAGTTGGCGGTCGCAGANAGCgat > 2:55936/1‑58 (MQ=255) gtgtaATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCAt < 2:298639/124‑1 (MQ=255) gtgtaATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACg < 1:447616/224‑1 (MQ=255) | GGTGAAACGACTCCTGGTTTGTACAAAGCCAGCAGTTCGTCCACCTTTTTACGTGTGCCGCCGTTGCTGCTTATACTGCGTCGAACTTTGACGACATGCAATTTTGCGCGCTGATACCACTCACGCGTTAACATCGTATCGTGATTGGAGATCAGCACTGGAATATGGCGCTCAACCAGACCTTCGGCGATCTCCGCCAGATGCGCTTGTTGTTCAAGCGTAAAACTGTTTGTGTGATACGCCGTAAAGTTGGCGGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACG > NC_000913/3515090‑3515544 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |