Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,429,067 | C→T | 12.9% | S244S (AGC→AGT) | PP_3044 → | lambda family portal protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,429,067 | 0 | C | T | 12.9% | 72.6 / 4.0 | 31 | S244S (AGC→AGT) | PP_3044 | lambda family portal protein |
Reads supporting (aligned to +/- strand): ref base C (12/15); new base T (2/2); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01 |
CCCCAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCTT > NC_002947/3428920‑3429216 | ccccAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCt < 2:82658/150‑1 (MQ=255) ttCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGAGAACTacgacga > 1:57558/1‑150 (MQ=255) ttCGCTGAACGCCGGTTACAACCCGCTNGTGCGTGTCCCGNCCNATCNGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTacgacgg > 1:457261/1‑149 (MQ=255) gCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGTCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGAcgcgagacat < 2:279872/150‑7 (MQ=255) ggTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCTCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTc < 1:54052/150‑1 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc < 1:279225/145‑1 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc < 1:279301/145‑1 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc > 2:279301/1‑145 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGTCTGGACAACTACGACGACGCGGTGCTCTTcc > 2:279225/1‑145 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAg > 1:343537/1‑150 (MQ=255) aaCCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAg > 1:343580/1‑150 (MQ=255) gtCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTg < 2:448255/150‑1 (MQ=255) cGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCg < 1:249343/150‑1 (MQ=255) aTCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCt < 2:343537/150‑1 (MQ=255) aTCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCt < 2:343580/150‑1 (MQ=255) aTGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCAAACTTGTTCGCCGGCTTTGtt > 2:124025/1‑150 (MQ=255) tgctgcACATCTTCGAGCCGTTGGAGCCTCGCGAATCCCGTGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTc < 2:40888/150‑1 (MQ=255) tgcACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCa > 2:156416/1‑150 (MQ=255) cGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATg < 1:124025/150‑1 (MQ=255) ttGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCgg > 1:330846/1‑122 (MQ=255) ttGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTCGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCgg < 2:330846/122‑1 (MQ=255) tGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAACTGGCCACCTTGTTAGACGGCTTTCTATACAGACCTGCGCCGGATGGGCTGGGTgg > 2:310303/1‑150 (MQ=255) gcccggTGGCCCCCGGCTTGCCCCGATCCTCTATTGCCTGCCCATCCTGGAAACCTACGAAGCCCCCTATAACTACCGGCAGTAAGTGGCCAACTGGTGCGCCGGCTTTGTTCGCAAACCGGCGCCGGATCGGCTGGGTGGCCCGcccct < 2:57558/146‑4 (MQ=255) ccGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCAc < 1:224159/150‑1 (MQ=255) ggCGCCGATCCTTAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACcggcgg < 2:75489/150‑1 (MQ=255) gcCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGc > 1:320782/1‑150 (MQ=255) gATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGgt > 2:34027/1‑150 (MQ=255) aTCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGgtt < 2:172040/150‑1 (MQ=255) tCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTTTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTc > 2:222983/1‑150 (MQ=255) tCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTc < 1:310303/150‑1 (MQ=255) ccTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCg > 2:382988/1‑150 (MQ=255) ccTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCtt < 1:382988/150‑1 (MQ=255) | CCCCAGCGACAAGGCTTCGCTGAACGCCGGTTACAACCCGCTTGTGCGTGTCCCGGCCGATCAGATGCTGCACATCTTCGAGCCGTTGGAGCCTGGCCAATTGCGAGGTGTCCCCCGGCTGGCGCCGATCCTCAAGCGCCTGCGCAGCCTGGACAACTACGACGACGCGGTGCTCTTCCGGCAGGAAGTGGCCAACTTGTTCGCCGGCTTTGTTCGCAAACCGGCGCCGGATGGGCTGGGTGGCCCGCCAATGGACATGATCACCGGCGGAGCGGTTGTTCATGACCGTGATGCCTT > NC_002947/3428920‑3429216 |
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: ‑ |