Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,328,579 | C→A | G134C (GGT→TGT) | mreB ← | cell wall structural complex MreBCD, actin‑like component MreB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,328,579 | 0 | C | A | 97.4% | 130.4 / ‑5.4 | 38 | G134C (GGT→TGT) | mreB | cell wall structural complex MreBCD, actin‑like component MreB |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (18/19); total (19/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
GCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAC > REL606/3328506‑3328651 | ncTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGGTCTTCAATCAGGAAGACTTCACGGGCACAAg > 1:147318/2‑76 (MQ=255) gCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAg > 1:661197/1‑76 (MQ=255) gCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAg < 1:1323742/76‑1 (MQ=255) gCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAg < 1:1389560/76‑1 (MQ=255) ttCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAgcg < 1:1382004/76‑1 (MQ=255) ttCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAgcg > 2:100191/1‑76 (MQ=255) gAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCt < 1:69348/76‑1 (MQ=255) aaCCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCCATCAGGAAGACTTCACGGGCACAAGCGCCCCgc > 2:120096/1‑76 (MQ=255) ccGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTgcgc < 1:853926/76‑1 (MQ=255) cAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGAt > 2:937416/1‑76 (MQ=255) aGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGAtt < 2:817509/76‑1 (MQ=255) ccAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTg < 2:194131/76‑1 (MQ=255) cGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCgcggcg > 1:693918/1‑76 (MQ=255) cGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCgcggcg > 1:666563/1‑76 (MQ=255) ggCAGCCATCGGTTCTTCAATCAGGCAGACCTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGt > 1:371212/1‑76 (MQ=255) cAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTc > 2:549487/1‑76 (MQ=255) gCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCaa < 2:813359/76‑1 (MQ=255) ccATCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCACTGCGCGGATTCACGAATTACGCGGCGCTCacc > 2:103387/1‑74 (MQ=255) aTCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAAccg > 2:1469979/1‑75 (MQ=255) tCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCGCGAATTGCGCGGCGTTCAACCTg < 2:1247302/76‑1 (MQ=255) tCGGTTCTTCAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTg < 2:1473421/76‑1 (MQ=255) ttcttcAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTc > 1:1316845/1‑76 (MQ=255) cAATCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGcc > 2:956253/1‑76 (MQ=255) aaTCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCa < 1:952934/76‑1 (MQ=255) aTCAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCaa > 2:253988/1‑76 (MQ=255) cAGGAAGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAAcc < 1:879571/76‑1 (MQ=255) aaGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGca > 2:874075/1‑76 (MQ=255) aaGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGca > 2:1257018/1‑76 (MQ=255) aGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGcac < 1:1594490/76‑1 (MQ=255) aGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGcac < 1:1481943/76‑1 (MQ=255) aGACTTCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGcac < 1:1097672/76‑1 (MQ=255) tCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACaaa < 1:1659021/76‑1 (MQ=255) tCACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACaaa > 2:1480188/1‑76 (MQ=255) cACGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAAc < 1:1186285/76‑1 (MQ=255) aCGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAAcc > 2:988074/1‑76 (MQ=255) cGGGCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCa > 2:1204110/1‑76 (MQ=255) gggCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAg < 1:464244/76‑1 (MQ=255) gCACAAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGaa < 1:1588893/76‑1 (MQ=255) cacaAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAc < 1:1441129/76‑1 (MQ=255) | GCTTCAGAAACCGGCAGGCCAGCACCAATTGCGGCAGCCATCGGTTCTTCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAC > REL606/3328506‑3328651 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |