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 | . | 89.4% | 179.0 / 19.7 | 47 | intergenic (‑167/‑156) | yifB/ilvL | predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide |
Reads supporting (aligned to +/- strand): ref base T (0/5); new base . (13/29); total (13/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT > REL606/3911891‑3912037 | gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt < 1:179163/76‑1 (MQ=255) aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt > 1:1319516/1‑76 (MQ=255) tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt > 2:712850/1‑76 (MQ=255) tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt > 2:712832/1‑76 (MQ=255) tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt < 2:1432783/76‑1 (MQ=255) atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt > 1:1072004/1‑76 (MQ=255) atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt < 1:80882/76‑1 (MQ=255) tgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttc < 2:74251/76‑2 (MQ=255) cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac > 1:261881/1‑73 (MQ=255) cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTtcactcac < 1:1063074/76‑1 (MQ=255) aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTcactcact > 2:931373/1‑76 (MQ=255) aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTcactcact < 1:1762029/76‑1 (MQ=255) aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTcactcact < 1:1278279/76‑1 (MQ=255) aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTcactcact < 1:1261223/76‑1 (MQ=255) aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtat < 2:1028851/76‑1 (MQ=255) aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtat < 1:759758/76‑1 (MQ=255) aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtat < 2:1734148/76‑1 (MQ=255) tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatttt < 2:273998/76‑1 (MQ=255) tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatttt > 1:226822/1‑76 (MQ=255) tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtatttt < 1:1329031/76‑1 (MQ=255) cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACtattttat > 2:68856/1‑76 (MQ=255) ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTAttttt < 2:1259332/76‑1 (MQ=255) ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTAttttt > 2:21827/1‑76 (MQ=255) ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTAttttt < 2:1056376/76‑1 (MQ=255) cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTaaa > 1:273731/1‑76 (MQ=255) ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAaac < 1:1126178/76‑1 (MQ=255) gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAacaaca < 1:1299646/76‑1 (MQ=255) gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAacaaca < 1:271242/76‑1 (MQ=255) aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAacaaca < 1:1754885/76‑1 (MQ=255) aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa < 1:563339/76‑1 (MQ=255) aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt < 2:1429860/76‑1 (MQ=255) caaCGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt < 2:1394935/74‑1 (MQ=255) cTACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt < 1:1667924/76‑1 (MQ=255) tACGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt < 1:1743421/76‑1 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat < 2:522216/76‑1 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat > 2:1129723/1‑76 (MQ=255) cGAGGAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat > 1:489959/1‑76 (MQ=255) aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata < 2:632568/76‑1 (MQ=255) gAAGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg > 2:572485/1‑76 (MQ=255) aGGGAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa > 1:492310/1‑76 (MQ=255) gAACAACATTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt > 1:383863/1‑76 (MQ=255) caacaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt > 2:1017221/1‑76 (MQ=255) aacaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta > 2:1059869/1‑76 (MQ=255) acaTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa < 1:402657/76‑1 (MQ=255) aTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 2:434042/76‑1 (MQ=255) aTTCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc < 1:705407/76‑1 (MQ=255) ttCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg > 2:546037/1‑63 (MQ=255) ttCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg < 1:546037/63‑1 (MQ=255) tCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc < 1:243227/76‑1 (MQ=255) tCATACTGAAATTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc < 2:381721/76‑1 (MQ=255) aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc < 2:122092/76‑1 (MQ=255) aaaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa < 1:1123167/76‑1 (MQ=255) aaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa < 1:495771/76‑1 (MQ=255) aaTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa < 1:174509/76‑1 (MQ=255) aTTGAA‑TTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa < 2:399428/76‑1 (MQ=255) tGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt < 1:198328/76‑1 (MQ=255) | GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT > REL606/3911891‑3912037 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |