Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,517,282 | G→C | 11.1% | V158L (GTG→CTG) ‡ | murE → | UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,517,282 | 0 | G | C | 11.1% | 122.4 / 8.1 | 45 | V158L (GTG→CTG) ‡ | murE | UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase |
Reads supporting (aligned to +/- strand): ref base G (19/21); new base C (2/3); total (21/24) | |||||||||||
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.35e-01 |
CACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTC > NC_002947/1517137‑1517425 | cACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg > 3:288285/1‑148 (MQ=255) ccAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg < 8:218976/146‑1 (MQ=255) ccAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg > 7:218976/1‑146 (MQ=255) cAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAAt > 7:58117/1‑149 (MQ=255) aaCGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATc < 2:224316/149‑1 (MQ=255) cgccAAGACAAGTGTTACCACACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGATTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgg < 8:236694/146‑2 (MQ=255) aaGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgggtcg < 8:152938/149‑6 (MQ=255) ttACCCAACTGGTCGCCAAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTc < 8:45100/149‑1 (MQ=255) ccAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCATAAGGCGAGTTGAGATAGGCCGCCGCCTGACCACCCCCGACCAGATTGGCGTGCAATCGACCCTCAACCACCTc > 1:256411/1‑148 (MQ=255) gACCTGCTTGGCCTGCGCTGCGGGCTGATTGGCACCCTGGGTACCGTTTTC‑TACTGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCGCGACCGGATTGCCGTGGTATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGc < 2:259998/149‑1 (MQ=255) gACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTgg < 4:288285/148‑1 (MQ=255) ccTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCa > 4:54407/1‑149 (MQ=255) cTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg < 1:145916/99‑1 (MQ=255) cTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg > 2:145916/1‑99 (MQ=255) cTTGGCCAGCGCTGCGGCCTGATTGGCAACCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATgg > 8:146203/1‑148 (MQ=255) cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTc > 1:232599/1‑149 (MQ=255) cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATc < 8:258672/95‑1 (MQ=255) cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATc > 7:258672/1‑95 (MQ=255) gcTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACTACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTcc < 8:58117/149‑1 (MQ=255) gCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGccac < 2:299146/124‑3 (MQ=255) gCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGccac > 1:299146/1‑122 (MQ=255) aTTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgtgt > 3:122646/1‑133 (MQ=255) aTTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgtgt < 4:122646/133‑1 (MQ=255) gCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg < 5:327354/71‑1 (MQ=255) gCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg > 6:327354/1‑71 (MQ=255) cACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAg < 2:256411/149‑1 (MQ=255) aCCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCGAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGc > 3:177628/1‑149 (MQ=255) gTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCa > 8:6985/1‑149 (MQ=255) tC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACt < 3:54407/148‑1 (MQ=255) c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg < 3:216347/148‑1 (MQ=255) c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg < 1:119567/148‑1 (MQ=255) c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg < 2:232599/148‑1 (MQ=255) aCGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCa > 5:223379/1‑127 (MQ=255) aCGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCa < 6:223379/127‑1 (MQ=255) cGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTgg < 7:146203/146‑1 (MQ=255) ttG‑CAGA‑GCGGCGGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACAt > 4:197429/1‑149 (MQ=255) cAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATcgc > 1:169108/1‑149 (MQ=255) a‑gCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGt > 5:254166/1‑149 (MQ=255) a‑gCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGt < 7:6985/149‑1 (MQ=255) aCCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCGAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCt < 4:177628/149‑1 (MQ=255) ccgacccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGAGATGACAAACCTGTCCcgc > 8:236248/1‑149 (MQ=255) tcgaTTGCACGGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACCTCGCGGTGATGACCAACCTGTCCCGCGccca > 2:230862/2‑149 (MQ=255) ccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCaa < 6:225257/134‑1 (MQ=255) ccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATCCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCaa > 5:225257/1‑134 (MQ=255) ttGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTc < 6:254166/149‑1 (MQ=255) | CACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTC > NC_002947/1517137‑1517425 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |