Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 2,131,466 | (A)9→8 | intergenic (‑125/+175) | gatY ← / ← fbaB | D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,131,458 | 0 | A | . | 91.3% | 87.8 / 2.9 | 23 | intergenic (‑117/+183) | gatY/fbaB | D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base . (15/6); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACA > REL606/2131383‑2131529 | aaatgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCa < 2:252085/76‑1 (MQ=255) gaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaa > 1:345373/1‑76 (MQ=255) aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa > 2:910146/1‑76 (MQ=255) aTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaa > 1:452014/1‑76 (MQ=255) tttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaa < 2:452014/76‑1 (MQ=255) ttGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaa > 2:283135/1‑76 (MQ=255) tGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaac > 2:110201/1‑75 (MQ=255) tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg > 1:1126454/1‑73 (MQ=255) ttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGcc > 1:85126/1‑76 (MQ=255) tAAATGTCGAGAGAACGATCACAAAAACGACAATATGAAAATTAATAGAGGNATGAAAGGC‑AAAAAAAACGGCCTc < 2:207061/76‑1 (MQ=255) aaTGTCGAGAAAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTccc > 1:162301/1‑76 (MQ=255) cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAgg > 1:735973/1‑76 (MQ=255) cGAGATAACGATCACAAAAACGAAAATATGAAAATTAATAGAGGAATGAAAGGC‑ATAAAAAACGGCCTCCCGGAgg > 1:865193/1‑76 (MQ=255) agaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGa > 1:822131/1‑76 (MQ=255) tAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGc > 1:485472/1‑76 (MQ=255) aaCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGC‑GAAAAAAACGGCCTCCCGGAGGGAAGcc > 1:682013/1‑76 (MQ=255) cGATCACAAAANCGACAATATGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGc > 1:175789/1‑76 (MQ=255) cacaAAAACAACAAAATAAAAATAAATAGAGAAATAAAAGGC‑AAAAAAAACGGCCCCCCGGAGGGAAGCCGcaaca < 2:827317/76‑1 (MQ=255) caAAAACGACAATATGAAAATTAATAGAGGAATGAAAGG‑AAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACaaa < 2:735973/76‑1 (MQ=255) cAATATGAAAATTAAGTGAGGAATGAAAGGC‑AAGTAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg < 2:772095/76‑1 (MQ=255) tGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt < 2:442361/76‑1 (MQ=255) tGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt > 1:245131/1‑76 (MQ=255) tGAAAATTAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGtttt > 2:690572/1‑76 (MQ=255) aaaaaTAAAAAGA‑GAAAGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt < 2:916121/75‑1 (MQ=255) aaaaTTAATAGAGGAATGAAGTGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTAt < 2:969732/76‑1 (MQ=255) ttAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGCTTTTATGCAg < 1:162773/76‑1 (MQ=255) ttAAAAGAGAAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg < 2:162301/76‑1 (MQ=255) tAATAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg > 2:231604/1‑76 (MQ=255) tAGAGGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGACGCAACAAAGTGCGCGTGGTTTTATGCAGGATg > 1:921341/1‑76 (MQ=255) gacgAATGAAAGGC‑AAAAAAAACGGCNTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcg > 2:998004/4‑76 (MQ=255) agGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAGGTGCGCGTGGTTTTATGCAGGATgcgg > 1:483893/1‑76 (MQ=255) agGAATGAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcgg < 2:715015/76‑1 (MQ=255) aaTGAAAGG‑AAAAACAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGt < 2:721583/76‑1 (MQ=39) gAAAGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTaaa < 2:17599/76‑1 (MQ=255) aaaGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAAc > 1:10692/1‑76 (MQ=255) aaGGC‑AAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACa > 2:164128/1‑76 (MQ=255) | AAATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGCGTAAACA > REL606/2131383‑2131529 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |