Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,733,719 | 0 | C | G | 32.3% | 50.7 / 17.5 | 31 | G63G (GGC→GGG) | gumD | ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base C (4/17); new base G (10/0); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTAT > REL606/3733645‑3733794 | tttGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCg < 2:901637/76‑1 (MQ=255) ngtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGTTCggg > 1:167991/2‑75 (MQ=255) tggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGa < 1:1404593/76‑1 (MQ=255) tggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGa < 2:1435354/76‑1 (MQ=255) ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt < 2:1022905/76‑1 (MQ=255) ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt < 2:1394372/76‑1 (MQ=255) gATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATg < 1:646546/76‑1 (MQ=255) tCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTc < 1:1027818/76‑1 (MQ=255) tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGGGGTGTCgg > 2:1193083/1‑75 (MQ=255) tGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGa < 2:1326472/76‑1 (MQ=255) gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGaa < 1:569857/76‑1 (MQ=255) gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGaa < 2:1171762/76‑1 (MQ=255) aaTATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGGCGGCGGTGGCGAAg > 2:580092/1‑76 (MQ=255) atatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGTTCGGCGATGTCGAAGc > 2:189031/1‑76 (MQ=255) tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTCGAAGCg > 1:19505/1‑76 (MQ=255) atCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGTTTGGCGATGGCGAAgggg > 1:75982/1‑73 (MQ=255) cGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGGCGAAGCgggt > 2:259338/1‑73 (MQ=255) cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGGGGCGGAGCGGtttt > 1:1176811/1‑76 (MQ=255) aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTGGGAGCGGTTTTTCa > 2:1487420/1‑76 (MQ=255) aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa < 1:1350910/76‑1 (MQ=255) cTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGGAGAGGTTGGCGGTGTCGGAGCGGGTTTCCaa > 2:402119/1‑75 (MQ=255) tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGTTCGGCGATGTCGAAGCGATTTTTCACg > 2:1019252/1‑76 (MQ=255) tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACg < 1:7159/76‑1 (MQ=255) atatGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGa < 1:615000/76‑1 (MQ=255) atGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGGAGAGGTCGGCGATGTCGGAGCGATTTTCCACGAAg > 1:216157/1‑76 (MQ=255) tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAgg < 2:1250011/76‑1 (MQ=255) ggATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTTGGGGGTGGCGAAGCGGATTTTCACGAAggg > 1:987879/1‑75 (MQ=255) aaGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTCGGCGATGTCGAAGCGGGTTTTCACGAAGGCGCAt > 1:1027511/1‑76 (MQ=255) aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc < 2:1397489/76‑1 (MQ=255) cTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGATTTTTCACGAAGGCGCATGCTCTTcc > 1:1321710/1‑76 (MQ=255) cTTCCTGATCCAGATTATGGCTGACGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcc > 1:780913/1‑76 (MQ=255) cTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccacca < 1:1502644/76‑1 (MQ=255) gATCCAGATTATGGCTGGCGGAGAGGTTGGCGGTGTCGAAGCGGTTTTTCACCAAAGCGCCTGCTCTTCCcccacc > 1:1475498/1‑76 (MQ=255) aTCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCg < 2:860014/76‑1 (MQ=255) aGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTg < 1:415989/76‑1 (MQ=255) aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg > 2:850120/1‑76 (MQ=255) ggctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGc < 1:420959/76‑1 (MQ=255) ggctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGc > 2:254752/1‑76 (MQ=255) ctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCaa < 2:227537/76‑1 (MQ=255) cGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGtat < 2:96310/76‑1 (MQ=255) | TTTGTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTAT > REL606/3733645‑3733794 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |