Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,328,417 | G→A | R188C (CGC→TGC) | 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,417 | 0 | G | A | 98.1% | 191.2 / ‑3.9 | 54 | R188C (CGC→TGC) | mreB | cell wall structural complex MreBCD, actin‑like component MreB |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (18/34); total (18/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01 |
TCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATC > REL606/3328344‑3328486 | tCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGcaca > 2:142266/1‑76 (MQ=255) cACCGATCAGAGAACCGTAATTACGACGCACATAGTTGCTGATAGCTTCGTCGAAATGGTCACCACCAATGCTcag > 2:1986696/1‑75 (MQ=255) ccGATCAGAGAACCGTAATTACGACGCAGATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACag < 1:306423/76‑1 (MQ=255) ccGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACag < 1:909823/76‑1 (MQ=255) cGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACaga < 1:680273/76‑1 (MQ=255) agagAACTGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg < 2:821498/76‑1 (MQ=255) agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg < 1:1071013/76‑1 (MQ=255) agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg > 2:423653/1‑76 (MQ=255) gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGt < 2:1516035/76‑1 (MQ=255) aaCCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTa > 2:1379402/1‑76 (MQ=255) ccGTAATTACGNCGNACATAGTTGATCATAGCTTCGTCGAAACGGTCACCACCAATGCTCACAGAAGAGGAGTaaa < 2:306970/76‑1 (MQ=255) ccGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTaaa < 1:1378567/76‑1 (MQ=255) ccGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTaaa > 2:1567492/1‑76 (MQ=255) gTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAAcc > 1:1379880/1‑76 (MQ=255) gTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAAcc < 1:634158/76‑1 (MQ=255) tAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACca < 2:1400054/76‑1 (MQ=255) aaTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcac < 2:1854378/76‑1 (MQ=255) aaTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcac < 1:489151/76‑1 (MQ=255) aTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcaca < 1:945512/76‑1 (MQ=255) tACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACAcc > 1:1336626/1‑76 (MQ=255) acgCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTc > 1:376541/1‑76 (MQ=255) gCACATAGTTGATGNTAGCTTCGTCGAAACGGTCACCACCAATGCACTCAGAAGAGGAGTAAACCAAACCNTTCaa < 2:2043244/76‑1 (MQ=255) acaTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAgg > 1:765540/1‑76 (MQ=255) aTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGgag < 1:773418/76‑1 (MQ=255) tAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGgaga > 1:576525/1‑76 (MQ=255) aGTTGATGCTAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt > 1:865994/1‑76 (MQ=255) aGTTGATGATAGNTTNGTCGAAACGGTCACCACCNATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt < 1:1668368/76‑1 (MQ=255) aGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACNCAGAAGAGGAGTAAACCANACNGTTCAAGGAGAt > 1:1658290/1‑76 (MQ=255) aGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt < 1:904758/76‑1 (MQ=255) gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGANATa < 1:1823185/76‑1 (MQ=255) gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATa > 1:875508/1‑76 (MQ=255) gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATa > 1:562527/1‑76 (MQ=255) ttGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa < 1:1263478/76‑1 (MQ=255) tgatAGCTTCGGCGAAACGGTCACCACCAATGCACACTGAAGAGGAGTAAACCACACCGTTCAAGTAGATAATAGc > 1:53298/1‑76 (MQ=255) gatAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCa < 1:1567493/76‑1 (MQ=255) tAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAAc > 2:1297560/1‑76 (MQ=255) tAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAAc > 2:1586534/1‑76 (MQ=255) cTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc < 2:1555972/76‑1 (MQ=255) ttCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCa < 1:280848/76‑1 (MQ=255) cgtcgAAACGGTCANCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGACAACAGCAACTTCAGt < 1:1748405/76‑1 (MQ=255) gtcgAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTg < 1:1263724/76‑1 (MQ=255) gtcgAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTg > 1:1014468/1‑76 (MQ=255) cgAAACGGTCACCACCAATGCACACAGNAGAGGAGTAAAACACACCGTTCAAGGAGATAACAGCAACTTCAGTGGt > 1:1784116/1‑76 (MQ=255) gAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTa < 2:1265877/76‑1 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac < 2:425661/76‑1 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac > 2:551135/1‑76 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac < 1:1140040/76‑1 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac < 1:80503/76‑1 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac < 1:442121/76‑1 (MQ=255) aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac < 2:1337320/76‑1 (MQ=255) aaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTacc < 1:973218/76‑1 (MQ=255) aCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTacca < 1:686626/76‑1 (MQ=255) ccaccaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGa < 2:652510/76‑1 (MQ=255) accaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGata < 2:1482079/76‑1 (MQ=255) ccaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGatat < 2:905720/76‑1 (MQ=255) caATGCACACAGACGAGGAGNAAACCACACCGTTCANGGAGATAACAGCAGCGTCAGTCGTACCACCACCGATATc < 2:2018476/76‑1 (MQ=37) | TCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATC > REL606/3328344‑3328486 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |