| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,639,013 | G→C | 18.2% | A18A (GCG→GCC) | PP_4950 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,639,013 | 0 | G | C | 18.2% | 49.4 / 6.7 | 22 | A18A (GCG→GCC) | PP_4950 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (7/11); new base C (2/2); total (9/13) | |||||||||||
| 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.34e-01 | |||||||||||
CAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATGGGCGCCAGGCGGAACAGC > NC_002947/5638870‑5639155 | cAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCCAg < 1:609542/150‑1 (MQ=255) aaGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAgc < 1:303595‑M2/150‑3 (MQ=255) ccAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCt > 1:215409/1‑150 (MQ=255) aTCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCgg < 1:220382/150‑1 (MQ=255) gtggtgCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACa < 1:270409/150‑1 (MQ=255) gtgCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTg > 2:280669/1‑150 (MQ=255) tCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggttt > 1:125690‑M2/1‑120 (MQ=255) tCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggttt < 2:125690‑M2/141‑22 (MQ=255) ccTGTGTGTTCCGACAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGttt < 1:280669/150‑1 (MQ=255) tgtTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCa > 2:600448/1‑150 (MQ=255) ttCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggtttcgaa < 1:277548‑M2/119‑26 (MQ=255) ttCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCTCGCCCGCCTGAAGGGCGGCGTgcacctccccgaccagggtttcgaa > 2:277548‑M2/1‑94 (MQ=255) ttCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAg < 2:468323/146‑1 (MQ=255) ttCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAg > 1:468323/1‑146 (MQ=255) gaCGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCa < 1:600448/150‑1 (MQ=255) tgaAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAgccaggcca < 2:215409/150‑1 (MQ=255) gaAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAgccaggccag > 1:63414/1‑150 (MQ=255) gCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAgcgc > 2:473797/1‑150 (MQ=255) cccTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTgg < 2:63414/150‑1 (MQ=255) tGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTgc > 2:42946/1‑150 (MQ=255) gggAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATg < 1:413654/150‑1 (MQ=255) ttCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGACATCCGGGCCAGCGCGCAACTGCGGCTGGTGCGCGATGGGCGCCAGGCGGAAcagc < 2:303595/150‑1 (MQ=255) | CAAGGGCGAACCCAATCAGGTGGTGCGGGTCGGCCACGCCTCGCCGGCCTGTGTGTTCCGCCAGATCGATGTATTCGGAGGAGACGCCTGATGAAGCACGCTTTCGAAACCCTGGTCGGGGAGGTGCACGCCGCCCTTCAGGCGGGCGAGCAGTTCACCCTCGGTTACAGTGCCGAGCAATCGCAGTTCGTGCGTTTCAACCATGCCAAGGTGCGCCAGGCCGGCGAAGTCAGCCAGGCCAGCGCGCAACTGCGGCTGGTGCGCGATGGGCGCCAGGCGGAACAGC > NC_002947/5638870‑5639155 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |