| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 395,861 | T→C | 14.8% | intergenic (+40/+164) | 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,861 | 0 | T | C | 14.8% | 60.5 / 4.2 | 27 | intergenic (+40/+164) | fdhA/PP_0329 | formaldehyde dehydrogenase/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (11/12); new base C (2/2); total (13/14) | |||||||||||
| 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.29e-01 | |||||||||||
AGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTGCGCAGTCATTTT > NC_002947/395718‑396009 | aGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAAAGCCCCTCGATG‑CGAGGGGCtcttttta < 1:155240/150‑4 (MQ=255)aGGTGATCAACCTGGATCAGGCGCAGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCACCGTCGAGTCAAAACAGCCCCTCGATG‑CGAGGAGCtcttttta > 2:155240/1‑147 (MQ=255) aTCAACCTGGATCAGGCGCCGGAAGGGTATGGAGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTc < 2:228950/150‑1 (MQ=255) tGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGATTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAg < 1:188267/150‑1 (MQ=255) ccGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGCGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGcc < 2:521749/150‑1 (MQ=255) gAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGa < 1:601875/150‑1 (MQ=255) ggCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGt < 1:4452/150‑1 (MQ=255) gCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTGTTGTCCATTCTGATAAGACATCTCGGCAAGACATCCTTGTCATAAGCCTTTc > 1:38860/1‑150 (MQ=255) aTTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATc > 2:37218/1‑136 (MQ=255) aTTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATc < 1:37218/136‑1 (MQ=255) ttCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAg > 1:262743/1‑150 (MQ=255) tCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAgc < 2:262743/150‑1 (MQ=255) ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCt > 1:364623/1‑117 (MQ=255) ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCt < 2:364623/117‑1 (MQ=255) ccGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGcac > 1:710392/1‑150 (MQ=255) gTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCa < 1:752559/94‑1 (MQ=255) gTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCa > 2:752559/1‑94 (MQ=255) aGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGTCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCt > 2:13459/1‑150 (MQ=255) gTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTg < 2:710392/150‑1 (MQ=255) aaagCCCCTCG‑ATGCGAGGGGCTCTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGaa > 1:680544/3‑89 (MQ=37) aaagCCCCTCG‑ATGCGAGGGGCTCTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGaa < 2:680544/87‑1 (MQ=37) aagCCCCTCG‑ATGCGAGGGGCTCTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCAt > 1:399725/2‑78 (MQ=37) aagCCCCTCG‑ATGCGAGGGGCTCTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCAt < 2:399725/77‑1 (MQ=37) ccccTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGa < 1:451215/131‑1 (MQ=255) ccccTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGa > 2:451215/1‑131 (MQ=255) ccccTCGCAT‑CGAGGGGCTTTTTGACATCGCCGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAggg > 1:209314/1‑150 (MQ=255) gCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCAt < 2:244757/69‑1 (MQ=255) gCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCAt > 1:244757/1‑69 (MQ=255) t‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTgcgc > 1:38046/1‑150 (MQ=255) tttttGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTGCGCAGTCAtttt > 1:685536/1‑150 (MQ=255) | AGGTGATCAACCTGGATCAGGCGCCGGAAGGGTATGGCGAGTTTGATGCAGGTGTGCCGAAGAAATTCGTTATCGACCCGCACAAAATGTGGGGTGCGGCGTAATCAGCGTCGAGTAAAAAGAGCCCCTCGCAT‑CGAGGGGCTTTTTGTCCATTCTGATAAGACATCTCGGCCGGAGGTGCTTGTCATAAGCCTTTCATCGCTTTTGAAGTCGAGCGCCGCCCGCACATACCTCATCAGTTCCGAAGCCTGGTGAGTGGCCCTGGAGCGAAGGGCCACTGCGCAGTCATTTT > NC_002947/395718‑396009 |
| 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: ‑ |