Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,911,972 | (T)7→8 | 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,965 | 1 | . | T | 91.2% | 114.7 / 4.9 | 34 | intergenic (‑166/‑157) | yifB/ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Reads supporting (aligned to +/- strand): ref base . (1/2); new base T (17/14); total (18/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
CTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑TTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATA > REL606/3911899‑3912032 | cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑tttttt‑ttc > 1:497058/1‑74 (MQ=255) tCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑tttttt‑ttca < 2:875960/76‑4 (MQ=255) tCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑tttttt‑ttca < 1:409435/76‑4 (MQ=255) aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑cactcact < 1:1169522/76‑1 (MQ=255) tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACtatt < 2:4847/76‑1 (MQ=255) cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACtatttt < 2:421676/76‑1 (MQ=255) aaaGGATGAACTATGAGGAAGGGGACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTaa < 1:1136556/76‑1 (MQ=255) aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑TTTTTTCT‑TACTCACTATTTTATTTTTaa < 1:6561/76‑1 (MQ=255) aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAaac > 1:1116687/1‑76 (MQ=255) aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTGATTTTTAAAAaac < 1:88695/76‑1 (MQ=255) gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAaacaa > 2:939252/1‑76 (MQ=255) aCTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAacaaca < 1:972692/76‑1 (MQ=255) cTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAAcaacaa < 1:232646/76‑1 (MQ=255) aCGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAAtt < 2:128283/76‑1 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAAttt < 2:497058/76‑1 (MQ=255) gAGGAAGGGAACAANATTCATACTGAAATTGA‑A‑TNTTTT‑TACANTCACTATTTTATTTTTAAACAACAACAATTta < 2:1091462/76‑1 (MQ=255) ggAAGGGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTtata > 1:1222102/1‑76 (MQ=255) aaGAGAACAACATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt > 2:981314/1‑76 (MQ=255) aGGGAACAACATTCATACTGAAATTGACT‑TTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTAATATTg > 1:1292326/1‑76 (MQ=255) aacaacATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat > 1:589019/1‑76 (MQ=255) aacaacATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat > 1:1005606/1‑76 (MQ=255) aacaacATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat > 1:1143781/1‑76 (MQ=255) aacaacATTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAat < 2:389913/76‑1 (MQ=255) aacaTTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat > 1:525585/1‑76 (MQ=255) caTTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta > 2:819944/1‑76 (MQ=255) caTTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta > 1:105699/1‑76 (MQ=255) aTTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa > 1:749077/1‑76 (MQ=255) aTTCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa < 2:1264447/76‑1 (MQ=255) tCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 2:385838/76‑1 (MQ=255) tCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 1:743372/76‑1 (MQ=255) tCATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 1:1030861/76‑1 (MQ=255) cATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg < 1:888999/76‑1 (MQ=255) cATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg > 1:510503/1‑76 (MQ=255) cATACTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg > 1:1051939/1‑76 (MQ=255) aCTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAACNCANCAATTTATATTGAANTTATTAAACGcat > 2:1188991/1‑76 (MQ=255) aCTGAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat > 1:287809/1‑76 (MQ=255) gAAATTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat > 2:976104/1‑76 (MQ=255) aaaTTGA‑ATTTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa > 2:833419/1‑76 (MQ=255) | CTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGA‑A‑TTTTTT‑T‑CACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATA > REL606/3911899‑3912032 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |