Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 413,107 | C→G | 15.0% | A686A (GCC→GCG) | glnE → | glutamate‑ammonia‑ligase adenylyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 413,107 | 0 | C | G | 15.0% | 67.2 / 5.8 | 28 | A686A (GCC→GCG) | glnE | glutamate‑ammonia‑ligase adenylyltransferase |
Reads supporting (aligned to +/- strand): ref base C (13/10); new base G (2/2); total (16/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01 |
GCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC > NC_002947/412969‑413256 | gCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtct < 1:458309/150‑4 (MQ=255) cTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCaa > 2:344262/1‑150 (MQ=255) tGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAAc > 1:477246/1‑150 (MQ=255) ggAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCTACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCTCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGa > 1:192118/1‑150 (MQ=255) ggCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtc < 1:84435/109‑3 (MQ=255) ggCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCGGGCCACggtc > 2:84435/1‑107 (MQ=255) ggCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtc > 1:387300/1‑150 (MQ=255) aaCCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtcatc < 1:595417/150‑1 (MQ=255) tgatgaAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGCCCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGtt < 1:344262/150‑1 (MQ=255) aaGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCa > 1:352927/1‑150 (MQ=255) gACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGggc > 2:444624/1‑150 (MQ=255) gacctgGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGa < 2:163699/150‑1 (MQ=255) ctggctggCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCt > 2:297207/1‑150 (MQ=255) ggctggttGACGCCACCCTCGACCAGGTGCTGGCCCTGGCCGGGCGCCAGACCGTTGCCCTCCACGGCCAACCCAAGCGCAGCGACGGCCGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCg < 2:192118/150‑1 (MQ=255) tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCCTTATCGGTTACGGCAAGATGGGCGGCCTGGAGCCCGGCCACgg > 2:294956/1‑150 (MQ=255) tgAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACgg > 1:94226/1‑150 (MQ=255) aTCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggac > 2:3073/1‑150 (MQ=255) cTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCATCTTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCggacctg < 2:298890/150‑1 (MQ=255) cAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg < 1:297207/150‑1 (MQ=255) tggccctggacAGGTGCCATACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGCCAGCAGGTTCTCCGACCCTGTTTTCATCATTATCGGTTCCGCCACGATGAGCGGCCTGGCGCTCGGTCACGGTTCGGACCTGGACCTGGTCTTCa < 1:444624/150‑1 (MQ=255) gccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtca > 1:532872/1‑79 (MQ=255) gccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTtca < 2:532872/79‑1 (MQ=255) ccctggccTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATcc > 2:601160/1‑150 (MQ=255) ctggcGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACgg < 2:387300/150‑1 (MQ=255) gcgcCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCAGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTGCATCCATGACGGCGa > 1:216758/1‑150 (MQ=255) ttggccGTGGCGGGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCg > 1:321924/5‑150 (MQ=255) aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg < 1:540116/125‑1 (MQ=255) aGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGgacctggacctg > 2:540116/1‑125 (MQ=255) ccGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGCCTTCATCCATGACGGCGATCCGCAGGCGGAAAc > 1:272494/1‑150 (MQ=255) | GCTGGCCCACAGCCTGCGCGTGGCCGCCTCGGAAATCAGCGGCAACCTGCCGTTGATGAAAGTCAGCGACTACCTGACCTGGCTGGCTGAAGCCATCCTCGACCAGGTGCTGGCCCTGGCCTGGCGCCAGACCGTGGCCCGCCACGGCCAACCCAAGCGCAGCGACGGCAGCCTGTGCGACCCGGGTTTCATCATTATCGGTTACGGCAAGATGGGCGGCCTGGAGCTCGGCCACGGTTCGGACCTGGACCTGGTCTTCATCCATGACGGCGATCCGCAGGCGGAAAC > NC_002947/412969‑413256 |
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: ‑ |