Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,070,720 | G→T | 9.5% | P72Q (CCG→CAG) | PP_4466 ← | TauD/TfdA family dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,070,720 | 0 | G | T | 9.5% | 130.7 / 2.6 | 42 | P72Q (CCG→CAG) | PP_4466 | TauD/TfdA family dioxygenase |
Reads supporting (aligned to +/- strand): ref base G (19/19); new base T (2/2); total (21/21) | |||||||||||
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.66e-01 |
GCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAGCATCGTAGTCGTTGAA > NC_002947/5070573‑5070862 | gCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCggg > 2:329680/1‑150 (MQ=255) gtggACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGaa > 1:134944/1‑150 (MQ=255) ggACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATa > 2:211661/1‑150 (MQ=255) gCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAg > 2:474190/1‑150 (MQ=255) cAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAg > 1:427802/1‑150 (MQ=255) gTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGtga < 2:294872/150‑1 (MQ=255) aTCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGtgatg < 2:499484/150‑1 (MQ=255) tCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGtgatga < 1:474190/150‑1 (MQ=255) tgCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAAc < 1:329680/150‑1 (MQ=255) tcgcttcgctGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGAtgt < 1:542683/150‑1 (MQ=255) tGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCACGAAGTGATGAGAGAACTTGATGTGCTGCTCAt > 1:65317/1‑150 (MQ=255) gTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCTGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGa < 1:442737/137‑1 (MQ=255) gTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCTGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGa < 1:442751/137‑1 (MQ=255) gTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCTGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGa > 2:442751/1‑137 (MQ=255) gTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCTGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGa > 2:442737/1‑137 (MQ=255) gTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCACGAAGTGATGAGAGAACTTGATGTGCTGCTCATc < 2:65317/150‑1 (MQ=255) ttACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCg < 2:130940/150‑1 (MQ=255) cAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAAt < 2:534069/150‑1 (MQ=255) aaaGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATg > 1:562794/1‑150 (MQ=255) aaGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTCCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGt > 1:154848/1‑150 (MQ=255) aGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTc < 2:134944/150‑1 (MQ=255) aGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTc > 2:410093/1‑150 (MQ=255) aGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTc < 2:562794/150‑1 (MQ=255) ttACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGtgctg > 2:287058/1‑122 (MQ=255) ttACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGtgctg < 1:287058/122‑1 (MQ=255) tCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTaca > 1:57079/1‑150 (MQ=255) ttAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTgg > 2:8703/1‑150 (MQ=255) ttAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTgg > 2:293006/1‑150 (MQ=255) aGCCATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCACTTGCAGCTGTTTCGGGTGAATAACAGCTTTCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAg < 2:255773/150‑1 (MQ=255) aaaTTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTc < 1:556920/150‑1 (MQ=255) tttttGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGtgatga > 2:247847/1‑67 (MQ=255) tttttGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGtgatga < 1:247847/67‑1 (MQ=255) gcgcCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAg > 2:41244/1‑148 (MQ=255) gcgcCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAg < 1:41244/148‑1 (MQ=255) tctTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAGCATCGTat > 1:476327/1‑149 (MQ=255) ttGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTc < 1:512893/119‑1 (MQ=255) ttGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTc > 2:512893/1‑119 (MQ=255) ttGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGc < 1:80906/129‑1 (MQ=255) ttGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGc > 2:80906/1‑129 (MQ=255) ttGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAGCATCGTGGTc > 2:234028/1‑150 (MQ=255) cTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAGCATCGTAGTCGTTGaa < 1:393518/150‑1 (MQ=255) tttCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTg > 1:240753/1‑96 (MQ=255) tttCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTg < 2:240753/96‑1 (MQ=255) | GCTGGTGGACGCTCGTAGAACCAGGTATCGGTGTGCCAGGTCGCTTCGCTGTTACCCAAAGCACCCGATGGCTTACCGTCTTTAGTAGTGTTAGCAATAACCAAAATTTCTTTTTGCGTTGGATGAGCGCCCTCTTGCAGCTGTTTCGGGTGAATAACAGCTTCCCCGAAGTGATGAGAGAACTTGATGTGCTGCTCATCGGTAATGTCGATATCGCGGAAACGTACAACACCGTACTGGAGCCAGATTTTTCGTACTTCGGCGACATCCTGAGCATCGTAGTCGTTGAA > NC_002947/5070573‑5070862 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |