Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 2,921,460 | G→A | G3494R (GGG→AGG) | PP_2561 → | hemolysin‑type calcium‑binding bacteriocin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,921,460 | 0 | G | A | 81.8% | 23.0 / ‑5.3 | 13 | G3494R (GGG→AGG) | PP_2561 | hemolysin‑type calcium‑binding bacteriocin |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (8/1); minor base C (0/2); total (8/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.45e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCAT > NC_002947/2921329‑2921599 | aGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGa > 5:147617/1‑148 (MQ=255) gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACg > 8:11539/1‑148 (MQ=255) cacgccggCGCGGGCGAGGACAGCCTCGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGGAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg < 1:106784/142‑1 (MQ=255) ggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGCCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATC‑AGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGa < 3:43170/148‑1 (MQ=255) tGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTACAGATACGGTGGTGGGTGGGGc > 4:7202/1‑148 (MQ=255) cggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTaa > 2:30966/1‑147 (MQ=255) aacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGaca > 4:35047/1‑148 (MQ=255) gcgAGGGTCACGACACTTTGTTCGGCGACGGGGGTACCGACACCCTGATCTGGGGGGCCGGCAACGACTTCCTCTAGGGTGCGGGAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGa < 2:61286/148‑1 (MQ=255) aCGACCCCTTGTTCGGCGGCGGGGGTACCGACACCCTGATCCGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa < 3:144389/148‑1 (MQ=255) cacaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTc > 6:172194/1‑148 (MQ=255) ccGTGGGTCCCGACACCCTCATCCGGGTGGCCGGCAACGACTTCCTCAACGGGGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTgg < 5:135121/147‑1 (MQ=255) cGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTACCGCCCCCATGATCGCCACCTATGGAACCGATGTGGTCCCGTGCGCTGCAGGCTGTGGCAACGccc > 8:158372/1‑146 (MQ=255) ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat > 1:72276/1‑147 (MQ=255) | AGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCAT > NC_002947/2921329‑2921599 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |