Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,507,655 | G→C | 10.9% | A290P (GCG→CCG) | cycA → | D‑alanine/beta‑alanine/D‑serine/glycine permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,507,655 | 0 | G | C | 10.9% | 129.7 / 7.4 | 46 | A290P (GCG→CCG) | cycA | D‑alanine/beta‑alanine/D‑serine/glycine permease |
Reads supporting (aligned to +/- strand): ref base G (22/19); new base C (2/3); total (24/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01 |
GCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCG > NC_002947/5507507‑5507798 | gCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCcggcggc < 1:1205398/150‑1 (MQ=255) gcgGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCACGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCggcggcggc < 1:537377/150‑1 (MQ=255) gCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGc < 2:73489/150‑1 (MQ=255) aCTCGATCCCGATCCGTATCATCGTGTTCGATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGAAACCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGATGTTAGATCTGAAGGGCGTGCAGGAGGCGGCAAGCAtcatca > 2:472889/1‑150 (MQ=255) cTCGACCCCGATCCGTATCATCGTGTTCTATGTGCT‑GCGTGGAATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCaa < 1:374376/150‑1 (MQ=255) atcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTc > 2:1273460/1‑150 (MQ=255) tcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCg > 1:1015427/1‑150 (MQ=255) ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg < 2:250318/150‑1 (MQ=255) atccATATCAGCGTGTGCTATGTGCGGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtggtg > 2:1090154/1‑150 (MQ=255) cctcGTATTCAATGTGCTGGCGCTG‑ATCGCGCTCATAGACATGAACCCGTGGCGCGATGTGGTGCCCGGAAAAAGCAAGTTCGTTGCGTTGTTCGTGCTGGCAGGCTTGAAGACGGCGGCAAACATAATCAACTTCGTGGTGCTGACCTc > 2:1360803/3‑150 (MQ=255) atcGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGc < 1:1176639/142‑1 (MQ=255) atcGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGc > 2:1176639/1‑142 (MQ=255) cTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGcc > 1:821495/1‑150 (MQ=255) gcgcTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc < 1:423718/143‑1 (MQ=255) gcgcTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc > 2:423718/1‑143 (MQ=255) gcgcTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAAc > 1:853230/1‑150 (MQ=255) cgcgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgt < 1:1273460/150‑1 (MQ=255) gATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGtt < 2:817646/150‑1 (MQ=255) cATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTtctc < 1:591290/150‑1 (MQ=255) tGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCTTCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTcca < 1:982851/150‑1 (MQ=255) tGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTcca < 1:983031/150‑1 (MQ=255) cccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCAt > 1:1144575/1‑150 (MQ=255) ccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg > 1:1227179/1‑93 (MQ=255) ccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg < 2:1227179/93‑1 (MQ=255) ggtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTc < 2:1144575/150‑1 (MQ=255) ggAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATg < 1:1044054/129‑1 (MQ=255) ggAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATg > 2:1044054/1‑129 (MQ=255) gAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGct < 2:59434/130‑1 (MQ=255) gAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGct > 1:59434/1‑130 (MQ=255) aaGCCCGTCCGTCGAGATCTACGTGCGTGCATGCTTTCCTTCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGa < 2:1015427/150‑1 (MQ=255) cccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcggc < 1:1403712/78‑1 (MQ=255) cccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcggc > 2:1403712/1‑78 (MQ=255) ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGGTTCGTCGGCCAACAGCGGCGTGTTCTCCACCCGCCGCATGCTCTACGGCCGGTCGCAGGAAGGTGAcgcg > 2:379692/1‑150 (MQ=255) ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc < 2:284523/86‑1 (MQ=255) ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc > 1:284523/1‑86 (MQ=255) ttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGccc > 2:1376117/1‑150 (MQ=255) gttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACAGCCTGTCGCAGGAAGGTTACGCGCCCAAg > 2:249552/1‑150 (MQ=255) taGTTGTTCCTGCTGGCAGCCTGGCCGGTGGCGGCAAGAATCATCAACTTCGTTGTGCTGACCTCGGTGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCGATACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCt < 1:472889/149‑1 (MQ=255) ttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTc < 1:1326135/150‑1 (MQ=255) ttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTc > 2:1333980/1‑150 (MQ=255) ttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCgag > 1:676245/1‑150 (MQ=255) tGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGt > 1:549588/1‑150 (MQ=255) tgatgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGcc > 1:796000/5‑106 (MQ=255) tgatgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGcc < 2:796000/102‑1 (MQ=255) tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGt > 2:806229/2‑106 (MQ=255) tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGt < 1:806229/105‑1 (MQ=255) gCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGcc > 1:646290/1‑150 (MQ=255) ccGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCg < 2:853230/150‑1 (MQ=255) | GCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTG‑ATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCG > NC_002947/5507507‑5507798 |
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: ‑ |