Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,052,930 | A→C | 14.7% | L198R (CTG→CGG) ‡ | PP_5305 ← | epimerase/dehydratase‑like protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,052,930 | 0 | A | C | 14.7% | 89.1 / 12.1 | 41 | L198R (CTG→CGG) ‡ | PP_5305 | epimerase/dehydratase‑like protein |
Reads supporting (aligned to +/- strand): ref base A (21/14); new base C (3/3); total (24/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.79e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01 |
ATCGGACCATTCGGTAACCCCCATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGACCACGCTCGCCGGGATGC > NC_002947/6052791‑6053074 | aTCGGACCATTCGGTAACCCCCATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTggcggc > 2:257531/1‑150 (MQ=255) ccATTCGGTAACCCCCATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTc < 2:121580/150‑1 (MQ=255) aTTCGGTAACCCCCATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCgg > 2:284386/1‑150 (MQ=255) ccccATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCCTGCTggcgtttctg < 2:717752/136‑7 (MQ=255) ccccATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCg > 2:124609/1‑150 (MQ=255) ccccATGTAGGCCCGCAGCCAGGCAACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCCTGCTggcgtttctg > 1:717752/1‑130 (MQ=255) cccATGTAGGCCCGCAGCCAGGCCACCACATCGGAAAGCCGCGCCGAGTTGTCAGCAACGCCAATGTAGCAGGCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCgg > 1:15111/1‑150 (MQ=255) ccATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAgg < 2:616767/122‑1 (MQ=255) ccATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAgg > 1:616767/1‑122 (MQ=255) aTGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTg < 2:383322/150‑1 (MQ=255) cGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAgc > 1:703447/1‑150 (MQ=255) gCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAgcg > 2:258680/1‑150 (MQ=255) aGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGCAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAgcggc < 1:258680/150‑1 (MQ=255) ccAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCCTGCTggcgtttc > 2:587130/1‑112 (MQ=255) ccAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCCTGCTggcgtttc < 1:587130/116‑5 (MQ=255) aGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGct > 2:309899/1‑150 (MQ=255) caccacATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATTGGCCTGCAGCAGAAACGCCAGCCTGCTggcgttt > 1:518383/1‑106 (MQ=255) caccacATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATTGGCCTGCAGCAGAAACGCCAGCCTGCTggcgttt < 2:518383/109‑4 (MQ=255) caccacATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGctcctc < 1:268824/150‑1 (MQ=255) gcCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGa < 1:309899/150‑1 (MQ=255) gcCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGa < 1:284386/150‑1 (MQ=255) ccGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGAc < 1:124609/150‑1 (MQ=255) tcatcaACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTgg > 2:770876/1‑150 (MQ=255) tcaACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTc > 2:331370/1‑150 (MQ=255) aCGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGc < 2:539720/150‑1 (MQ=255) gCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCaa > 2:129350/1‑150 (MQ=255) gCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCaa > 1:632808/1‑150 (MQ=255) aGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAAc < 1:506661/143‑1 (MQ=255) aGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAAc > 2:506661/1‑143 (MQ=255) ccATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCAATAGATACCg > 1:331114/1‑150 (MQ=255) cATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCgg > 1:94806/1‑150 (MQ=255) cAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCaa > 2:262227/1‑87 (MQ=255) cAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCaa < 1:262227/87‑1 (MQ=255) aTCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATAc < 1:702469/140‑1 (MQ=255) aTCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATAc > 2:702469/1‑140 (MQ=255) gCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGCCGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGAcc < 1:257531/150‑1 (MQ=255) cagcagAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGACCACGCt > 1:112705/1‑150 (MQ=255) gCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGACCACGCTCGCCGGGATg > 1:65878/1‑150 (MQ=255) ccAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTaa > 1:600343/1‑71 (MQ=255) ccAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTaa < 2:600343/71‑1 (MQ=255) ccAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGACCACGCTCGCCGGGATGc > 2:338901/1‑150 (MQ=255) | ATCGGACCATTCGGTAACCCCCATGTAGGCCCGCAGCCAGGCCACCACATCGGCCAGCGGCGCCGGGTCGTCATCAACGCCAATGTAGCAGTCGTCCAGAGCCTTGCCATCAGCATCGGCCTGCAGCAGAAACGCCAGCAGGCTGGCGGCATCTTCGGCATGAATGCGGTTGCCATACAGCGGCGGCTCCTCGGCAACGCGGTAACCCTGACGCACCTGGCTCAGCAACCACTCGCGCCCTGGGCCATAGATACCGGTCAGGCGGACCACGCTCGCCGGGATGC > NC_002947/6052791‑6053074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |