Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,307,446 | C→G | 16.7% | G60R (GGC→CGC) | livK ← | branched‑chain amino acids ABC transporter substrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,307,446 | 0 | C | G | 16.7% | 67.3 / 9.4 | 30 | G60R (GGC→CGC) | livK | branched‑chain amino acids ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (11/14); new base G (3/2); total (14/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01 |
GACGCTGGCTGGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAG > NC_002947/1307298‑1307587 | gACGCTGGCTGGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGccgccc < 1:845790/150‑1 (MQ=255) aCGCTGGCTGGGTGGAGCTGGAGCACANGTGGCCGACCACNAACTTGNCGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGccgcccg > 1:777187/1‑150 (MQ=255) gctggctgGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCcgcg < 2:507659/150‑1 (MQ=255) gctggGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGt < 1:63329/150‑1 (MQ=255) ctggGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGtt < 2:135706/150‑1 (MQ=255) tggagcACAGGTGGCCGATCACGAACTTGACGCCGTCGTTGACCACTTTGTTCGCGACTGCCACGGCCTGTTTATGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGa < 2:271606/150‑1 (MQ=255) tggagcACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGa > 1:282825/1‑150 (MQ=255) cGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGtg > 2:740959/1‑123 (MQ=255) cGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGtg < 1:740959/124‑2 (MQ=255) tGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCa < 2:278018/150‑1 (MQ=255) gtcgtTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGcc < 2:640714/150‑1 (MQ=255) cACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGt > 2:790006/1‑150 (MQ=255) tttGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCAc < 2:187694/150‑1 (MQ=255) gCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTg < 1:790006/150‑1 (MQ=255) tttAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTgg < 1:501745/121‑1 (MQ=255) tttAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTgg > 2:501745/1‑121 (MQ=255) tAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGc > 1:522569/1‑150 (MQ=255) tCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGtg < 1:757079/62‑2 (MQ=255) tCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCGGGCGGCGtg > 2:757079/1‑61 (MQ=255) tCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATc > 1:585196/1‑150 (MQ=255) aGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGa > 2:480783/1‑150 (MQ=255) gCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATg < 1:433664/150‑1 (MQ=255) gCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATg < 1:394950/150‑1 (MQ=255) ttGGCTTCGACCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGGCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAAc < 1:82485/150‑1 (MQ=255) cTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAt > 1:423351/1‑150 (MQ=255) cATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTAt > 1:570101/1‑128 (MQ=255) cATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTAt < 2:570101/128‑1 (MQ=255) cATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTAGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAAtggctggc > 2:294446/1‑150 (MQ=255) aTCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGc < 2:585196/150‑1 (MQ=255) aTCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGc > 2:387947/1‑150 (MQ=255) caaCGCCGCGGGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCa > 2:835559/3‑150 (MQ=255) cACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAg > 1:28592/1‑150 (MQ=255) | GACGCTGGCTGGGTGGAGCTGGAGCACAGGTGGCCGACCACGAACTTGACGCCGTCGTTGACCACTTTGTTGGCGACTGCCACGGCCTGTTTAGGGTCGCAGGCGTCGTCGTATTCCTTGGCTTCGAGCATCTTGCCATCCACGCCGCCCGCGGCGTTGATATCCTTGATGGCTTGCTTGGCGCCGATGAACTGCATGTCGCCGTACTGGGTCACGGGACCCGTCTTGGGACCGGCGATGCCGATCTTGATGGTATCGGCGGCAAACGAATGGCTGGCTACCCCGGCCAG > NC_002947/1307298‑1307587 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |