Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 395,849 | Δ1 bp | 11.1% | intergenic (+28/+176) | fdhA → / ← PP_0329 | formaldehyde dehydrogenase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 395,849 | 0 | C | . | 11.1% | 97.5 / 8.9 | 36 | intergenic (+28/+176) | fdhA/PP_0329 | formaldehyde dehydrogenase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (16/15); new base . (2/2); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01 |
TGAAGTGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAA > NC_002947/395701‑395988 | tGAAGTGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCa > 2:1279582/1‑150 (MQ=255) tGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAg < 2:1269126/150‑1 (MQ=255) aGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGTTCCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGGGTCGAGTAAAAAGAGCCCCTCGATG‑CGAGGGGc < 2:550394/150‑1 (MQ=255) ggTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTcc < 2:1298627/150‑1 (MQ=255) gATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCAtt < 1:1279582/150‑1 (MQ=255) cTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCG‑ATGCGAGGGGCtc < 1:340384/133‑2 (MQ=255) cTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCG‑ATGCGAGGGGCtc > 2:340384/1‑132 (MQ=255) gATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTCTCGACCCGCACAAAATGTGCGGTGCGCCGTAATCAGCGTCGAGTAAACAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGCTAAGAc > 2:438811/1‑150 (MQ=255) ggAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCgg > 2:821841/1‑150 (MQ=255) aaGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCCCCCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGGCCTTCTCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAg < 1:1327232/150‑1 (MQ=255) tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTCAAGAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAg > 1:318036/1‑150 (MQ=255) tttGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAg < 1:821841/150‑1 (MQ=255) gtgtGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt < 2:463143/132‑1 (MQ=255) gtgtGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt > 1:463143/1‑132 (MQ=255) gCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCttt < 2:318036/150‑1 (MQ=255) ccGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCtttt > 1:1418984/1‑150 (MQ=255) cGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt < 2:40239/126‑1 (MQ=255) cGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCtt > 1:40239/1‑126 (MQ=255) cGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTg > 2:591816/1‑150 (MQ=255) gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCttt < 2:426592/87‑1 (MQ=255) gaagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCttt > 1:426592/1‑87 (MQ=255) aagaaATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTTaa < 1:1422346/150‑1 (MQ=255) ttATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGcgcc > 1:932748/1‑150 (MQ=255) gACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGccgcccg > 2:1267358/1‑150 (MQ=255) ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAATAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGTCTTTACTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGcac > 1:1156069/1‑150 (MQ=255) gCACCAAATGTGGGGTGCGGGGTAATGAGTGGCGAGTAAACAGGGCCCCTCGCAT‑CGAGTGGCTTTTTGTCCATTCTGATAAGACATCTAGGCCGGAGGAGCCTGCCATAAGCCCTTCACCGCTTTTCACGTCGCGCCGCGCCCGCACAt < 2:1156069/150‑1 (MQ=255) gCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGcaca < 2:219920/149‑1 (MQ=255) gCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGcaca > 1:219920/1‑149 (MQ=255) cgTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGcc > 2:947833/1‑150 (MQ=255) gTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTg < 1:591816/150‑1 (MQ=255) tAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTgg > 2:205744/1‑150 (MQ=255) agagCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTg < 1:939052/87‑1 (MQ=255) agagCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTg > 2:939052/1‑87 (MQ=255) aaagCCCCTCG‑ATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTg > 1:638686/3‑131 (MQ=255) aaagCCCCTCG‑ATGCGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTg < 2:638686/129‑1 (MQ=255) gagCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCGTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGaa < 1:639732/150‑1 (MQ=255) | TGAAGTGGTAGGTGTGCAGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAA > NC_002947/395701‑395988 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |