Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,911,972 | (T)7→6 | intergenic (‑173/‑150) | yifB ← / → ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,911,966 | 0 | T | . | 96.7% | 136.5 / ‑1.8 | 30 | intergenic (‑167/‑156) | yifB/ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base . (9/20); total (9/21) | |||||||||||
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.80e-01 |
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAA > REL606/3911891‑3912036 | gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt < 1:729151/76‑1 (MQ=255) aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt > 1:1000490/1‑76 (MQ=255) tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc < 2:758143/76‑2 (MQ=255) aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTcactcact < 1:663109/76‑1 (MQ=255) aaaTGCAGCGGACAAAGGATGAACTGCGNGGAAGGGAAAGACATTCAGACTGAAATTGAA‑TTTTTTCACTCACtat < 1:122418/76‑1 (MQ=255) aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatt < 2:1060047/76‑1 (MQ=255) aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtatt < 1:879770/76‑1 (MQ=255) gCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtattttattt < 2:494579/76‑1 (MQ=255) ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTAttttt < 2:836200/76‑1 (MQ=255) aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaa < 1:864117/76‑1 (MQ=255) cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaa > 2:610292/1‑76 (MQ=255) cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaa < 1:240615/76‑1 (MQ=255) aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaaaa < 2:29603/76‑1 (MQ=255) aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTGTTaaaaaa < 2:164394/76‑1 (MQ=255) ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAaac < 1:576058/76‑1 (MQ=255) tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAaacaac > 1:1185566/1‑76 (MQ=255) aaCTAGAAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa < 2:1335950/76‑1 (MQ=255) aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa < 2:1237887/76‑1 (MQ=255) aGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat < 1:401522/76‑1 (MQ=255) gAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg < 1:465843/76‑1 (MQ=255) aaGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa > 1:1102267/1‑76 (MQ=255) aGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa > 1:903414/1‑76 (MQ=255) gggAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa < 1:1103656/76‑1 (MQ=255) ggAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat < 1:254151/76‑1 (MQ=255) gAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt < 1:416694/76‑1 (MQ=255) aacaacATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta > 2:818497/1‑76 (MQ=255) acaacaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat < 1:899607/76‑1 (MQ=255) aTTCATACTGAAATTGAA‑NTTTTTCACTCACTATTTTATGTTTAAAGAACAACAATTTATATTGAAATTATTAAAc > 1:1339573/1‑76 (MQ=255) aTACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat > 2:749577/1‑76 (MQ=255) aCTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca < 1:1165797/76‑1 (MQ=255) cTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat > 1:657707/1‑76 (MQ=255) aaaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa > 2:997751/1‑76 (MQ=255) aTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa < 1:383648/76‑1 (MQ=255) | GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAA > REL606/3911891‑3912036 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |