Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 326,067 | C→A | 22.4% | *556Y (TAG→TAT) | PP_0269 ← | glutamate synthase large subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 326,067 | 0 | C | A | 22.4% | 60.8 / 12.5 | 32 | *556Y (TAG→TAT) | PP_0269 | glutamate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base C (9/15); new base A (3/4); total (12/20) | |||||||||||
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.44e-01 |
TCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAA > NC_002947/325920‑326207 | tCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTa < 1:150161/150‑1 (MQ=255) aaGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCgg < 2:231453/150‑1 (MQ=255) aGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCggg < 2:136993/150‑1 (MQ=255) gTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCCACCAAGTCCGCCCATTGGTCGAGTACCCCATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGtt < 1:189427/150‑1 (MQ=255) aGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTc < 2:67620/150‑1 (MQ=255) tCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCttt > 2:319179/1‑150 (MQ=255) acaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCg > 2:362245/1‑121 (MQ=255) acaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCg < 1:362245/121‑1 (MQ=255) caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCTGttcga > 2:502868/1‑128 (MQ=255) caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCTGttcga < 1:502868/131‑4 (MQ=255) aacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCTTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGt < 2:294183/150‑1 (MQ=255) gCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTg > 1:206676/1‑150 (MQ=255) cagcagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCgg < 1:321629/150‑1 (MQ=255) aaCGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTgg > 2:78817/1‑150 (MQ=255) aaCGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAAc < 1:201669/117‑1 (MQ=255) aaCGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAAc > 2:201669/1‑117 (MQ=255) acgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGt < 1:104221/150‑1 (MQ=255) gAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCt < 1:78817/150‑1 (MQ=255) gATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGtt < 2:220611/86‑1 (MQ=255) gATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGtt > 1:220611/1‑86 (MQ=255) cAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCcgcg < 1:319179/150‑1 (MQ=255) cATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACCTCCGCGCAt > 2:369620/1‑150 (MQ=255) ttctgattctgaTCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaaaa > 2:224029/1‑150 (MQ=255) tctgattctgatCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAAc > 2:327499/1‑150 (MQ=255) tGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGc < 2:242703/150‑1 (MQ=255) aagaaaCAGCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaaaa < 2:275516/126‑1 (MQ=255) aagaaaCAGCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaaaa > 1:275516/5‑130 (MQ=255) aTCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCAc < 1:327499/150‑1 (MQ=255) gAACAGCCCCGGCCTCTGCCGGGGCTGTTTCTTTTCGGCGACCGGTAACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCTATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGa < 1:362576/150‑1 (MQ=255) aGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTaa < 1:224029/150‑1 (MQ=255) ccGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACg > 2:101039/1‑66 (MQ=38) ccGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACg < 1:101039/66‑1 (MQ=38) | TCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAA > NC_002947/325920‑326207 |
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: ‑ |