Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 605,134 | T→C | 13.5% | G148G (GGT→GGC) | thiL → | thiamin monophosphate kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 605,134 | 0 | T | C | 13.5% | 82.7 / 6.2 | 37 | G148G (GGT→GGC) | thiL | thiamin monophosphate kinase |
Reads supporting (aligned to +/- strand): ref base T (21/11); new base C (2/3); total (23/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.46e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01 |
TGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGC > NC_002947/604989‑605276 | tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCaa > 1:232043/1‑150 (MQ=255) tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCaa < 1:529085/150‑1 (MQ=255) tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGACTTGCGCCGCAGCGGTGCaa > 1:135204/1‑150 (MQ=255) ggcggcCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCCGGCCTTGCACCGCTGCGGCGCAAg < 2:140844/150‑1 (MQ=255) cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCCGGCCTTGCACCGCTGCGGCGCAAGGCCtgg < 1:262532/150‑1 (MQ=255) ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg < 1:397422/150‑1 (MQ=255) ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg > 1:423750/1‑150 (MQ=255) gCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtg > 1:351401/1‑150 (MQ=255) cATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg > 2:539901/1‑123 (MQ=255) cATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg < 1:539901/123‑1 (MQ=255) ccACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCt > 2:168849/1‑150 (MQ=255) cACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg > 1:132398/1‑150 (MQ=255) cACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg > 1:43323/1‑150 (MQ=255) ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTAGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGc < 1:784866/150‑1 (MQ=255) ggtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg < 2:43323/150‑1 (MQ=255) acacCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGc > 2:791660/1‑150 (MQ=255) cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtgt < 1:376872/100‑1 (MQ=255) cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtgt > 2:376872/1‑100 (MQ=255) cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg > 1:217639/1‑150 (MQ=255) cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg > 1:355170/1‑150 (MQ=255) cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg > 2:381168/1‑150 (MQ=255) cgGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAgcg < 2:423750/150‑1 (MQ=255) gCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGGGAGcgcg > 1:265856/1‑150 (MQ=255) tGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc < 2:135204/150‑1 (MQ=255) gAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGcc < 1:381168/150‑1 (MQ=255) ccctgtTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCg > 1:797503/4‑150 (MQ=255) cgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGa > 1:793894/1‑150 (MQ=255) gtgtTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAg > 1:734751/1‑150 (MQ=255) atCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGa < 1:304442/118‑1 (MQ=255) atCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGa > 2:304442/2‑119 (MQ=255) gTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCt > 1:652220/1‑150 (MQ=255) gTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCt > 2:397764/1‑150 (MQ=255) ttGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGTGCCTGCGGCGCAGGCCGAGCCGCTGCTGGCTCATTACt < 2:232043/149‑1 (MQ=255) catgCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTAc < 1:631912/141‑1 (MQ=255) cAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTAc > 2:631912/1‑144 (MQ=255) ttGCACCGCTGCGGCGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATg > 2:285561/1‑150 (MQ=255) gcAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGc > 1:74814/1‑150 (MQ=255) | TGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGC > NC_002947/604989‑605276 |
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: ‑ |