Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,733,719 | 0 | C | G | 20.3% | 164.3 / 18.8 | 59 | G63G (GGC→GGG) | gumD | ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base C (1/46); new base G (12/0); total (13/46) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTA > REL606/3733648‑3733793 | gtgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAag < 1:909352/76‑1 (MQ=255) tgAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAaga < 1:1290597/76‑1 (MQ=255) gAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAagag < 2:207126/76‑1 (MQ=255) aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGagagg > 2:1719647/1‑75 (MQ=255) aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGagagg > 2:392747/1‑75 (MQ=255) aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAagagg > 2:1392110/1‑75 (MQ=255) aaCCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGt < 2:1989718/76‑1 (MQ=255) ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTc < 1:1200354/76‑1 (MQ=255) ccTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTc < 1:1257030/76‑1 (MQ=255) cTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCg < 1:847475/76‑1 (MQ=255) cTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCg < 1:107448/76‑1 (MQ=255) ggtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGc < 1:1192708/76‑1 (MQ=255) gtggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCg < 1:2039465/76‑1 (MQ=255) ggATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGAt < 2:1892763/76‑1 (MQ=255) aTCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGgtgt > 2:1041582/1‑76 (MQ=255) gAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTGGaa > 1:349965/1‑76 (MQ=255) aaTATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAg < 2:710336/76‑1 (MQ=255) tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCg < 2:236763/76‑1 (MQ=255) tatCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCg < 1:66025/76‑1 (MQ=255) cGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGAGGTCGGAGCGact > 1:1898166/1‑74 (MQ=255) cAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGAtttt < 2:593743/76‑1 (MQ=255) aGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGAAGAGGTCGGCGATGTCGACGCGAttttt > 2:969754/1‑75 (MQ=255) gACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGCGGAGCGGtttttc > 1:2013482/1‑74 (MQ=255) gACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTcc < 1:29644/76‑1 (MQ=255) aCTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCa < 1:1326506/76‑1 (MQ=255) aCTATATGGATAAGGAAGACTTCATGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGGGTCGGAGAGAgttttta > 1:518630/1‑73 (MQ=255) cTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTCGAAGGGGTTTTTCAc > 1:232603/1‑76 (MQ=255) tataTGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGGAGAGGTTGGCGGTGTCGAAGGGGTTTTCCACg > 2:115847/1‑76 (MQ=255) atGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTCGAAGCGCGTTTTTACGAAg > 2:668879/1‑76 (MQ=255) atGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAg < 1:740443/76‑1 (MQ=255) tGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAgg < 1:268086/76‑1 (MQ=255) aaGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCAt < 2:2029151/76‑1 (MQ=255) aGGAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGGCGGAGCGGTTTTTCACGAAAGCGCATg > 1:1167124/1‑76 (MQ=255) gAAGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGATGTCGAAGCGATTTTTCACGAAGGCGCATGct > 2:1993413/1‑76 (MQ=255) aaGACTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGGCGGAGGGAGTTTTCACGAAGGGGGATGctc > 1:1318490/1‑76 (MQ=255) aGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGctct < 2:149599/76‑1 (MQ=255) aGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGctct < 2:1280560/76‑1 (MQ=255) gACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCtt < 2:932002/76‑1 (MQ=255) gACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCtt < 2:2080614/76‑1 (MQ=255) aCTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGTCGAAGCGATTTTTCACGAAGGCGCACGCTCTTc > 1:2013206/1‑76 (MQ=255) aCTTCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGCGAAGCGATTTTTCACGAAGGCGCACGCTCTTc > 2:225210/1‑76 (MQ=255) aCTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTc < 2:1673289/76‑1 (MQ=255) cTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTTGAAGCGATTTTCCACGAAGGCGCATGCTCTTcc < 1:1896692/76‑1 (MQ=255) cTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcc < 2:1131736/76‑1 (MQ=255) cTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcc < 2:538867/76‑1 (MQ=255) ttCCTGATCCAGATTATGGCTGGGGGAGAGATCGGCGATGTCGAAGCGATTTTTCACGAAGGCGCATGCTCTTcca > 2:2113037/1‑76 (MQ=255) ttCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcca < 2:405343/76‑1 (MQ=255) ttCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTcca < 1:2008661/76‑1 (MQ=255) tCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGTGGAAGGGGTTTTTCACGAAGGGGGCTGGTCTTtcac > 2:546044/1‑76 (MQ=255) tCCTGATCCAGATTATGGCTGGGGGAGAGGTCGGGGGTGTCGGAGCGAGTTTTCACGCAGGGGCGCGCTCTTcccc > 2:927552/1‑74 (MQ=255) tCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccac < 2:1036205/76‑1 (MQ=255) cTGATCCAGATTATGGCTGGGGGAGAGGTTGGCGGTGTCGGAGCGGTTTTTCACGAAGGCGCATGCTCTTccccca > 2:1050079/1‑76 (MQ=255) cTGATCCAGATTATGGCTGGGGGAGAGGTCGGCGGTGGTGGAGCGGTTTTTCTCAAAAGCGCCTGCTCTTcccccc > 2:1601320/1‑75 (MQ=255) cTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccacca < 1:529790/76‑1 (MQ=255) cTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccacca < 2:1300716/76‑1 (MQ=255) cTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTccacca < 2:1965611/76‑1 (MQ=255) gATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCaccacc < 2:793124/76‑1 (MQ=255) gATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCaccacc < 1:95617/76‑1 (MQ=255) aTCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCg < 1:442475/76‑1 (MQ=255) aTCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCg < 2:1426240/76‑1 (MQ=255) ccAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAg < 1:1737202/76‑1 (MQ=255) cAGATTATGGCTGGGGAAGAGGTCGGCGGTGTTGAAGCGAATTTTCACCCAGGAGGCTGCTCTTTCCCCACCgcgg > 2:773804/1‑73 (MQ=255) gATTATGGCTGGGGAAGAGTTTGGCGATGTCGAAACGATTTTCTACGAAGGCGCAATCTTTTCCCCCACCGAGTgg > 2:1828821/1‑76 (MQ=255) aTTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTggg < 1:1242200/76‑1 (MQ=255) ttATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGa < 2:860314/76‑1 (MQ=255) tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAt < 1:1346156/76‑1 (MQ=255) tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAt < 1:1293774/76‑1 (MQ=255) tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAt < 1:801036/76‑1 (MQ=255) tATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGAt < 1:801048/76‑1 (MQ=255) aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg < 2:331701/76‑1 (MQ=255) aTGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATg < 2:1427895/76‑1 (MQ=255) gctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCa > 2:1333177/1‑76 (MQ=255) ctggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCaa > 1:1292900/1‑76 (MQ=255) tggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAg > 2:418782/1‑76 (MQ=255) ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt > 2:408903/1‑76 (MQ=255) ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt > 1:1612679/1‑76 (MQ=255) ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt < 1:741792/76‑1 (MQ=255) ggcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGt < 1:895669/76‑1 (MQ=255) gcGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGta < 2:1178697/76‑1 (MQ=255) | GTGAACCTGGTGGATCTGAATATCGCAGACTATATGGATAAGGAAGACTTCCTGATCCAGATTATGGCTGGCGAAGAGTTCGGCGATGTCGAAGCGATTTTCCACGAAGGCGCATGCTCTTCCACCACCGAGTGGGATGGCAAGTA > REL606/3733648‑3733793 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |