| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,645,948 | A→T | 15.5% | L31Q (CTG→CAG) | PP_3212 ← | Rieske 2Fe‑2S family protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,645,948 | 0 | A | T | 15.5% | 63.4 / 2.8 | 25 | L31Q (CTG→CAG) | PP_3212 | Rieske 2Fe‑2S family protein |
| Reads supporting (aligned to +/- strand): ref base A (9/12); new base T (2/2); total (11/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTG > NC_002947/3645806‑3646095 | ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaa > 1:166314/1‑150 (MQ=255) cAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAAt > 1:386107/1‑150 (MQ=255) ccGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaagcgc < 2:166314/150‑5 (MQ=255) ccGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGc > 2:132823/1‑150 (MQ=255) tACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTgggg < 1:11841/150‑1 (MQ=255) ccGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCTCGCATTATCACCTGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCg > 1:296998/1‑150 (MQ=255) gCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGt < 2:386107/150‑1 (MQ=255) ccTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCa > 2:628477/1‑150 (MQ=255) tGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCa > 2:732518/1‑150 (MQ=255) cGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAggg < 1:132823/150‑1 (MQ=255) cTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATAcc < 2:3933/150‑1 (MQ=255) gACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGaaa < 2:296998/150‑1 (MQ=255) gTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgc < 1:732518/150‑1 (MQ=255) tttGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgcg > 1:116046/1‑150 (MQ=255) gcggtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCGTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAgg < 2:761491/150‑1 (MQ=255) ggCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgttgtt < 2:736758/150‑1 (MQ=255) gcagcaGCGGTCCCGCATTGCCACCGGGGGACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt > 1:415947/1‑150 (MQ=255) tGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGTCGCAATACAGGTTGTTTGGTAACCAACATAATCAg < 2:347564/150‑1 (MQ=255) ccATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCg > 2:480612/1‑66 (MQ=255) ccATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCg < 1:480612/66‑1 (MQ=255) cATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCg > 2:547379/1‑150 (MQ=255) acaacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATgg < 1:628477/150‑1 (MQ=255) ccAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCa < 2:116046/150‑1 (MQ=255) cactCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg > 2:599449/4‑93 (MQ=255) cactCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg < 1:599449/90‑1 (MQ=255) gcagGGTGAATGGCTTGGGTCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGagcag < 2:779568/150‑5 (MQ=255) | CCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTG > NC_002947/3645806‑3646095 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |