| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 972,949 | G→A | 13.3% | R413H (CGC→CAC) | secD → | protein translocase subunit SecD |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 972,949 | 0 | G | A | 13.3% | 66.4 / 4.2 | 30 | R413H (CGC→CAC) | secD | protein translocase subunit SecD |
| Reads supporting (aligned to +/- strand): ref base G (13/13); new base A (2/2); total (15/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01 | |||||||||||
TTCATCGAGCAGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCTGGGTGCTGACAACATCAC > NC_002947/972807‑973081 | ttCATCGAGCAGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGAACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTAGCGACGATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCAcc > 1:6015/1‑150 (MQ=255) aGCAGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGAATCACCGGCTTGa > 1:563588/1‑150 (MQ=255) gCAGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCCACGAGTAAAAAAAGATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGaa > 1:379256/1‑150 (MQ=255) aGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATg > 1:475290/1‑150 (MQ=255) aGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGGCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCCTCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATg < 2:611133/150‑1 (MQ=255) aCCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGaa < 1:42835/150‑1 (MQ=255) ccGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATc < 1:315789/150‑1 (MQ=255) ccGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATc > 2:55273/1‑150 (MQ=255) gCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGACGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGGTTGAATGGCCAGGGCGAATCGt > 1:150250/1‑150 (MQ=255) cgTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGc < 1:55273/150‑1 (MQ=255) agaAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCggtg > 2:379600/1‑150 (MQ=255) gaAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCggtgg > 2:530949/1‑150 (MQ=255) gACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTg > 2:456728/1‑150 (MQ=255) tGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGc < 1:192937/150‑1 (MQ=255) aaaGCTTCGTCGCAGAGAAGACAATGATCGGAAAGGCGACCATGCAGAAGCTGCTAAGTAGCCAGTTCCGCATCACCGGCTAGAATGGCTTGGGCGAATATTCCGAGCTGGCCCTGCTGCTGCGTGGCGGTGGTTTGGGGCCACCGACGt > 2:181018/1‑150 (MQ=255) aaGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTgaag > 1:478556/1‑150 (MQ=255) aaaaaaaTCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGt > 1:58025/1‑150 (MQ=255) aaaaaaTCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTa < 1:530949/150‑1 (MQ=255) caGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGccc < 2:475290/150‑1 (MQ=255) caGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGccc < 2:478556/150‑1 (MQ=255) gCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCACATCACCGGCTTGAATGGCCa > 1:552372/1‑61 (MQ=255) gCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCACATCACCGGCTTGAATGGCCa > 1:552178/1‑61 (MQ=255) gCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCACATCACCGGCTTGAATGGCCa < 2:552178/61‑1 (MQ=255) gCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCACATCACCGGCTTGAATGGCCa < 2:552372/61‑1 (MQ=255) ggCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCt < 2:379256/150‑1 (MQ=255) gACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCTggg < 1:379600/150‑1 (MQ=255) tCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCTGGGTGCTGacaac < 1:456728/150‑1 (MQ=255) gccgcTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGt > 2:1518/1‑118 (MQ=255) gccgcTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGt < 1:1518/118‑1 (MQ=255) gcTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCTGGGTGCTGACAACATCAc < 2:58025/150‑1 (MQ=255) | TTCATCGAGCAGAAGCCGGTTACCCGCTACGTCAAGCAAGCCGTCGGCGGCGTCGAGAAGGACGTTGCCGTGCAAAGCTTCGTCGAAGAGAAAAAAATCATCAGCCTGGCGACCATCCAGTCGCCGCTTGGTAGCCAGTTCCGCATCACCGGCTTGAATGGCCAGGGCGAATCGTCCGAGCTGGCCCTGCTGCTGCGTGCCGGTGGTCTGGCTGCACCGATGTACTTCGCTGAAGAACGTACCATTGGCCCAAGCCTGGGTGCTGACAACATCAC > NC_002947/972807‑973081 |
| 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: ‑ |