Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,069 | G→A | A578V (GCG→GTG) | mrdA ← | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 649,069 | 0 | G | A | 97.8% | 172.2 / ‑3.6 | 46 | A578V (GCG→GTG) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (24/21); total (24/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01 |
CCAACCGCCGGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGC > REL606/649000‑649139 | ccACCCGCGGGCCCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAg < 2:1466564/75‑1 (MQ=255) cccgccgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCg < 2:337915/74‑1 (MQ=255) aCCGCCGGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCg > 2:592602/1‑75 (MQ=255) gccgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGtc > 2:1869079/1‑75 (MQ=255) cccgGCCCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGtc < 1:1291793/74‑1 (MQ=255) cgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTCTGGCACAAAGGCGGccat > 2:1895125/1‑75 (MQ=255) ggccccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGtcatc < 1:878807/72‑1 (MQ=255) gACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGtcatca > 1:1912719/1‑75 (MQ=255) aCCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAg > 1:654755/1‑75 (MQ=255) cccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGt < 2:2075171/75‑1 (MQ=255) ccccgccgcCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGtt < 2:1578712/72‑1 (MQ=255) cgcgCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGttt > 1:1715794/1‑75 (MQ=255) cgcgCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGttt > 2:899116/1‑75 (MQ=255) cgcgCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGttt > 1:731674/1‑75 (MQ=255) ccgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTgg > 2:1907735/1‑75 (MQ=255) gccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTc < 1:1392345/75‑1 (MQ=255) cgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCAc > 2:851663/1‑75 (MQ=255) ttCTCCAGAATCATGGCGACAGCCACTTGTGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGt > 2:488256/1‑75 (MQ=255) ttCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGt < 2:1828337/75‑1 (MQ=255) tctcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTa > 1:2003472/1‑75 (MQ=255) tctcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTa < 1:2019901/75‑1 (MQ=255) caccAGAATCAGGGCGGCAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTaa < 2:1856784/73‑1 (MQ=255) tcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTaaa < 2:2026830/75‑1 (MQ=255) gAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCt < 1:1343749/75‑1 (MQ=255) aaTCATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTc > 2:818321/1‑75 (MQ=255) cagggCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGc < 1:102081/72‑1 (MQ=255) cATGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGc < 1:816672/75‑1 (MQ=255) aTGGCGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCa < 2:858943/75‑1 (MQ=255) cggcAGCCACTTGCGGATTGTTGTATGGGACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAAttt < 2:1028047/72‑1 (MQ=255) cGACAGCCACTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAAttt < 1:1516490/75‑1 (MQ=255) aCTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCa > 1:1743095/1‑75 (MQ=255) aCTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCa > 2:95774/1‑75 (MQ=255) aCTTGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCa > 1:1329988/1‑75 (MQ=255) ttGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCAtt > 2:363440/1‑75 (MQ=255) tGCGGATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta < 1:1887427/75‑1 (MQ=255) gCGGATTGTTGTATGGCGCAAAGGCAGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATtat < 2:1171128/75‑1 (MQ=255) ggATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAg > 2:1513912/1‑75 (MQ=255) ggATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAg < 2:504041/75‑1 (MQ=255) ggATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAg < 2:504025/75‑1 (MQ=255) gATTGTTGTATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAgg > 1:1773754/1‑75 (MQ=255) tgttgtATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGttt > 1:2070711/1‑75 (MQ=255) gttgtATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttc < 1:1357567/75‑1 (MQ=255) tgtATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgt < 1:1125527/75‑1 (MQ=255) tATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGGAATTTTGTGCGCATTATAGGTttcgttc > 1:1579709/1‑75 (MQ=255) tATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgttc > 2:1853219/1‑75 (MQ=255) tATGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgttc < 1:1590061/75‑1 (MQ=255) aTGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgttcg > 2:1075697/1‑75 (MQ=255) tGGCACAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGc > 1:416307/1‑75 (MQ=255) | CCAACCGCCGGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGC > REL606/649000‑649139 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |