Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,097,189 | Δ2 bp | coding (383‑384/846 nt) | glpF ← | glycerol facilitator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,097,187 | 0 | A | . | 93.0% | 178.1 / 8.9 | 43 | coding (386/846 nt) | glpF | glycerol facilitator |
Reads supporting (aligned to +/- strand): ref base A (2/1); new base . (22/18); total (24/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.46e-01 | |||||||||||
* | REL606 | 4,097,188 | 0 | C | . | 93.0% | 180.5 / 10.4 | 43 | coding (385/846 nt) | glpF | glycerol facilitator |
Reads supporting (aligned to +/- strand): ref base C (2/1); new base . (22/18); total (24/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 = 7.54e-01 |
ACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAC > REL606/4097112‑4097263 || aCCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCaa > 2:1840528/1‑76 (MQ=255) cATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAact < 1:542619/76‑2 (MQ=255) tctcAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAacttt < 1:1054588/76‑4 (MQ=255) gCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTg > 1:1179940/1‑76 (MQ=255) gCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTg > 1:1643574/1‑76 (MQ=255) gCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTg > 2:1694721/1‑76 (MQ=255) gCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTACCAGCCAGATCA‑‑ACTTTCAACGCTg > 1:358185/1‑76 (MQ=255) gAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGcc > 1:994967/1‑76 (MQ=255) aaaGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCcg > 1:698827/1‑76 (MQ=255) aaaGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCcg > 2:1076929/1‑76 (MQ=255) aGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCcgcg > 1:1075225/1‑76 (MQ=255) gCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGa > 1:618383/1‑76 (MQ=255) cTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAAc > 2:854730/1‑76 (MQ=255) caAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAAtgt < 2:188097/76‑1 (MQ=255) caAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAAtgt > 1:170819/1‑76 (MQ=255) aaTTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGtgat < 2:1990371/76‑1 (MQ=255) aTTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGtgatg > 1:285860/1‑76 (MQ=255) ttGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGtgatga < 2:855911/76‑1 (MQ=255) tGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAg < 2:1993037/76‑1 (MQ=255) tGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAg > 1:737855/1‑76 (MQ=255) atGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTc > 1:1063113/1‑76 (MQ=255) ttggggTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAg < 2:761608/73‑1 (MQ=255) tAGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGt < 1:1493492/76‑1 (MQ=255) aGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAg > 2:1923818/1‑76 (MQ=255) aGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTc > 1:574921/1‑76 (MQ=255) aGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTc < 2:1770336/76‑1 (MQ=255) aGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTc < 2:1168365/76‑1 (MQ=255) aGGGTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTc > 2:1051218/1‑76 (MQ=255) gggTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCg < 2:46408/76‑1 (MQ=255) ggTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGa > 1:1960273/1‑76 (MQ=255) ggTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGa > 1:743770/1‑76 (MQ=255) ggTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGa < 1:1686523/76‑1 (MQ=255) gTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGaa > 2:1628544/1‑76 (MQ=255) gTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGaa < 2:1949854/76‑1 (MQ=255) gTAAGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGaa < 2:610534/76‑1 (MQ=255) aaGTAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGaaaa < 1:1592946/76‑1 (MQ=255) tAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAaaata < 1:1958748/76‑1 (MQ=255) tAGAGAAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAaaata > 1:2000721/1‑76 (MQ=255) agaAAGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAaaataaat > 1:596057/1‑76 (MQ=255) aaGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgt < 2:566174/76‑1 (MQ=255) aGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgta > 2:1310818/1‑76 (MQ=255) aGTGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgta < 1:176620/76‑1 (MQ=255) tGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgtagt < 2:875859/76‑1 (MQ=255) tGCCAGCCAGATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgtagt > 2:928698/1‑76 (MQ=255) ccagccagATCA‑‑ACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTaa < 2:618383/76‑1 (MQ=255) cAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTa < 2:655042/76‑1 (MQ=255) cacTTTCAACGCTGCCGCGAGCAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAc < 1:1825111/76‑1 (MQ=255) cacTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAc < 1:1825128/76‑1 (MQ=255) || ACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAC > REL606/4097112‑4097263 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |