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 | . | 98.0% | 236.9 / ‑1.4 | 50 | intergenic (‑167/‑156) | yifB/ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (22/27); total (23/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
ATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA > REL606/3911893‑3912034 | atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt > 2:385774/1‑76 (MQ=255) tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt < 2:480968/76‑1 (MQ=255) agatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt > 1:1172186/2‑76 (MQ=255) atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACTAGGAAGGGAACCACAGNCATACTGAAATTGAAtttttt > 2:279661/1‑76 (MQ=255) atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACTAGGAAGGGAACCACAGNCATAATGAAATTGACtttttt > 2:279680/1‑76 (MQ=255) cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac < 1:214347/76‑4 (MQ=255) gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTtcactca > 2:18095/1‑76 (MQ=255) aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACta < 1:1033998/76‑1 (MQ=255) aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtat < 1:992248/76‑1 (MQ=255) aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatt > 1:1621702/1‑76 (MQ=255) aTGCAGCGGACAAAGGATGAANCACGCGGAAGGGAACATCATTCATACTGAAATTGAA‑TTTTTTCACTCACtattt < 2:269460/76‑1 (MQ=255) aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtattt < 1:1657498/76‑1 (MQ=255) aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtattt < 1:28627/76‑1 (MQ=255) gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatttta < 1:275895/76‑1 (MQ=255) gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatttta < 1:1066573/76‑1 (MQ=255) cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtattttat < 2:352898/76‑1 (MQ=255) ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTAttttt > 2:873784/1‑76 (MQ=255) cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaa < 2:953324/76‑1 (MQ=255) cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaa < 2:586227/76‑1 (MQ=255) aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaaa > 2:633077/1‑76 (MQ=255) aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaaaa > 1:679414/1‑76 (MQ=255) ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAaac < 1:455879/76‑1 (MQ=255) gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAaaca < 1:170615/76‑1 (MQ=255) gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAacaaca < 1:1252735/76‑1 (MQ=255) gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAacaaca < 2:1173713/76‑1 (MQ=255) aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa > 2:1432086/1‑76 (MQ=255) aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa > 2:1282727/1‑76 (MQ=255) tACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt > 1:949676/1‑76 (MQ=255) aCGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta < 2:1324192/76‑1 (MQ=255) gAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata > 2:1369186/1‑76 (MQ=255) gAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata > 2:146041/1‑76 (MQ=255) gAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata < 1:1371135/76‑1 (MQ=255) gAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata < 2:221075/76‑1 (MQ=255) aGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat < 1:1267173/76‑1 (MQ=255) aGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat > 2:203903/1‑76 (MQ=255) gAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg < 1:1051201/76‑1 (MQ=255) aaGGGAACAACATTCATACTGAAATTGGA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa > 2:191794/1‑76 (MQ=255) aaGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa > 1:75225/1‑76 (MQ=255) aGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa < 1:1526946/76‑1 (MQ=255) aacaacATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta < 1:1035936/76‑1 (MQ=255) aacaacATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta > 2:461073/1‑76 (MQ=255) acaacaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat > 1:838724/1‑76 (MQ=255) acaacaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat > 2:1573867/1‑76 (MQ=255) acaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa < 1:1407300/76‑1 (MQ=255) caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa > 2:541341/1‑76 (MQ=255) aTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 2:1471653/76‑1 (MQ=255) aTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 2:595932/76‑1 (MQ=255) aTACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat < 1:385774/76‑1 (MQ=255) aTACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat > 2:607382/1‑76 (MQ=255) aTACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat > 1:993487/1‑76 (MQ=255) aCTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca > 1:227902/1‑76 (MQ=255) cTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat < 2:804847/76‑1 (MQ=255) tGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa > 2:898504/1‑76 (MQ=255) aaaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa > 2:1438307/1‑76 (MQ=255) aaaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa > 1:1173052/1‑76 (MQ=255) aaaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa < 2:1334479/76‑1 (MQ=255) | ATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA > REL606/3911893‑3912034 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |