Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 450,008 | A→T | I226N (ATC→AAC) | ylaB ← | conserved inner membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 450,008 | 0 | A | T | 91.2% | 197.4 / 10.7 | 57 | I226N (ATC→AAC) | ylaB | conserved inner membrane protein |
Reads supporting (aligned to +/- strand): ref base A (1/4); new base T (30/22); total (31/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01 |
AGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTG > REL606/449939‑450080 | aGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGaaa > 2:17879/1‑76 (MQ=255) gTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATAtt > 2:1385321/1‑76 (MQ=255) gTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATAtt > 2:703130/1‑76 (MQ=255) gTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATAtt > 2:643708/1‑76 (MQ=255) tAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTc > 2:275555/1‑76 (MQ=255) cAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCAtct > 2:1566997/1‑76 (MQ=255) aaaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCAtctc > 2:1373602/1‑76 (MQ=255) aaaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCAtctc < 2:575107/76‑1 (MQ=255) aTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCg > 2:477917/1‑76 (MQ=255) aTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCg > 2:477909/1‑76 (MQ=255) tAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCgg > 1:1409552/1‑76 (MQ=255) aaaGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCggg < 1:1106972/76‑1 (MQ=255) aaaGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCggg < 1:1511304/76‑1 (MQ=255) aGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGaa < 1:681906/76‑1 (MQ=255) aGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGaa < 1:681922/76‑1 (MQ=255) aGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCGGGaa < 1:1504862/76‑1 (MQ=255) cTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAggg < 2:1625492/76‑1 (MQ=255) ttGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGc < 2:1019861/76‑1 (MQ=255) ttGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGc < 1:952461/76‑1 (MQ=255) ttGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGc < 1:512300/76‑1 (MQ=255) ttGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGc < 1:797747/76‑1 (MQ=255) tGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCa > 2:350843/1‑76 (MQ=255) gCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAg < 1:777863/76‑1 (MQ=255) gATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGat < 1:882935/76‑1 (MQ=255) gATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCGGGAAGGGCAGGat < 1:211463/76‑1 (MQ=255) gCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATc < 2:445162/76‑1 (MQ=255) ccAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCg > 1:712328/1‑76 (MQ=255) aGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTa < 1:560636/76‑1 (MQ=255) gCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTaa < 1:1481370/76‑1 (MQ=255) gCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTaa > 1:1608751/1‑76 (MQ=255) ttatGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGa > 2:1131730/4‑76 (MQ=255) tCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGAt > 1:347327/1‑76 (MQ=255) tCTGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGAt > 1:930633/1‑76 (MQ=255) tGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTc > 1:1317190/1‑76 (MQ=255) tGCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTc > 2:1316266/1‑76 (MQ=255) gCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTcc > 2:1439503/1‑76 (MQ=255) gCAACATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTcc > 2:484026/1‑76 (MQ=255) gCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTcc > 2:363136/1‑76 (MQ=255) gCAACATTTTCGTTGAAGCCCAGGAGATGTTCGAAATATTCATCTCCGGGAGGGGCAGGATATCGTAAATGATTcc > 2:1439492/1‑76 (MQ=255) cATTTTCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCatta > 2:1015713/1‑76 (MQ=255) ttttCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCattatt < 1:826794/76‑1 (MQ=255) ttttCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCattatt < 1:1486980/76‑1 (MQ=255) ttCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTAtttt < 2:1197629/76‑1 (MQ=255) ttCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTAtttt > 1:258471/1‑76 (MQ=255) tCGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTc > 1:597915/1‑76 (MQ=255) cGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCg > 1:408742/1‑76 (MQ=255) cGTTGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCg > 1:256130/1‑76 (MQ=255) ttGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGat < 1:697730/76‑1 (MQ=255) tGAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGata < 1:505621/76‑1 (MQ=255) gAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGatat > 2:1474534/1‑76 (MQ=255) gAAGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGatat > 2:1085089/1‑76 (MQ=255) aaGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATg < 1:157450/76‑1 (MQ=255) aGCCCAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGc > 1:1113102/1‑76 (MQ=255) ccAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCa > 1:1510188/1‑76 (MQ=255) ccAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCa < 2:9031/76‑1 (MQ=255) cAGGTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCAc < 1:777893/76‑1 (MQ=255) ggTGATGTTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCg < 1:267064/76‑1 (MQ=255) tgatgtTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCggg > 2:1589351/1‑76 (MQ=255) tgatgtTCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCggg > 2:83188/1‑76 (MQ=255) atgttCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGgtg > 2:832628/1‑76 (MQ=255) | AGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAGGTGATGATCGAAATATTCATCTCCGGGAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTG > REL606/449939‑450080 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |