Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 780,563 | G→T | 15.6% | intergenic (‑90/+32) | PP_0670 ← / ← glyA‑II | bile acid/Na+ symporter family transporter/serine hydroxymethyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 780,563 | 0 | G | T | 15.6% | 81.7 / 5.9 | 32 | intergenic (‑90/+32) | PP_0670/glyA‑II | bile acid/Na+ symporter family transporter/serine hydroxymethyltransferase |
Reads supporting (aligned to +/- strand): ref base G (12/15); new base T (3/2); total (15/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.45e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01 |
GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGG > NC_002947/780414‑780704 | ggCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGCCTGCGTGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTgg > 1:439497/1‑150 (MQ=255) cggcggCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGct < 1:544311/150‑1 (MQ=255) tAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGAGCCTGGCCCAG‑ATAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTa < 2:222798/150‑1 (MQ=255) aTGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGc < 1:234878/150‑1 (MQ=255) ggATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAgc > 1:712838/1‑150 (MQ=255) ccGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATc < 2:182365/150‑1 (MQ=255) tttCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCggg < 1:796180/117‑1 (MQ=255) tttCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg > 1:10395/1‑150 (MQ=255) tttCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTg > 2:757748/1‑150 (MQ=255) tttAGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCgg > 2:796180/1‑117 (MQ=255) gAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTc < 1:11416/121‑1 (MQ=255) gAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCGGGCTCTTTc > 2:11416/1‑121 (MQ=255) tGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAAc < 2:69183/150‑1 (MQ=255) ctctctGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCg < 1:745965/150‑1 (MQ=255) gTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGctt > 1:102204/1‑150 (MQ=255) aaGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGcttctt > 1:130631/1‑150 (MQ=255) tG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAg > 2:329088/1‑150 (MQ=255) aaaGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGAcc < 1:455113/150‑1 (MQ=255) aaGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCCCAGTGCCggg > 1:195111/1‑81 (MQ=38) aaGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCt > 1:751004/1‑150 (MQ=255) gTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTtc > 2:101447/1‑150 (MQ=255) tCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCAcc < 1:101447/150‑1 (MQ=255) tACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCg < 2:712838/150‑1 (MQ=255) gCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg < 1:329088/150‑1 (MQ=255) gCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAg < 2:10395/150‑1 (MQ=255) gcCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGcc < 2:439497/150‑1 (MQ=255) tAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTc < 2:164381/150‑1 (MQ=255) aaaCGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCg > 2:190850/1‑150 (MQ=255) gagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCCGCAGCTTCTTGCAGATg < 1:253455/79‑1 (MQ=37) gagCCCGGCCCAGTGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCCGCAGCTTCTTGCAGATg > 2:253455/1‑79 (MQ=37) gCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGc < 2:751004/150‑1 (MQ=255) cccGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGcc > 2:42997/1‑150 (MQ=255) cGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGAGCACCGCTTCGTTGTTCAGGTCGGCCAGGATTTCGGAGATCCAGCCgg > 2:65099/1‑150 (MQ=255) | GGCCAGTACGGCGGCGAAGTAGATGGGGATCTGCCGGGTTTCGAGTTGCTCTCTGGTCAAGGTG‑CGGCCTGTGCAAAGTGGAGGTATTGGGGGCGCAGTCTTTGTACGAGGGGCGGTGCGCCTGTATTAAACGAAAGAGCCCGGCACTGGGCCGGGCTCTTTCACTTTAGGCAAGCGCCAATCAGTTGCCGTAAACCGGCAGCTTCTTGCAGATGGCCTTGACCTTCTCACGGACGGCGTCGATCACCGCTTCGTTGTTCAGGTCGGCCAGGATGTCGCAGATCCAGCCGG > NC_002947/780414‑780704 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |