Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,217,434 | 0 | G | C | 37.5% | 47.0 / 38.2 | 72 | A76P (GCC→CCC) | agaS | tagatose‑6‑phosphate ketose/aldose isomerase |
Reads supporting (aligned to +/- strand): ref base G (42/3); new base C (18/9); total (60/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 126 total aligned reads displayed.
AAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAAC > REL606/3217359‑3217508 | aaaGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCCTGGCTCg > 1:837489/1‑76 (MQ=255) aaGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGGGCCGTGGCTCgc > 1:2256280/1‑76 (MQ=255) aaGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgc < 1:6310145/76‑1 (MQ=255) agagAATCTGCGGATCATCCTGACCGGTGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgcc > 1:1541158/1‑76 (MQ=255) agagAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgcc > 2:3765699/1‑76 (MQ=255) gagaGTCTGCGGATCATCCTGAGCGGAGCCGGAACGGCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgcca < 2:4952336/76‑1 (MQ=255) gagaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgcca < 1:3033596/76‑1 (MQ=255) gagaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgcca < 2:3480607/76‑1 (MQ=255) agaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTcccccg > 2:6288846/1‑74 (MQ=255) agaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTcccccc > 1:3900017/1‑74 (MQ=255) agaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgccag < 2:1534764/76‑1 (MQ=255) agaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgccag < 1:6246871/76‑1 (MQ=255) gaATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGGGGCTccccgcc > 2:2404136/1‑73 (MQ=255) aaTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTcgcccgcc > 1:1681121/1‑76 (MQ=255) aaTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgccagcc < 1:1972492/76‑1 (MQ=255) aaTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCgccagcc < 2:1770543/76‑1 (MQ=255) aTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCccccgccc > 2:4937815/1‑75 (MQ=255) aTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCccccgcca > 1:6860335/1‑76 (MQ=255) tCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTACCGGTGACATCACCGCGCCGTGGCTCGCCAGCCAt < 1:520131/76‑1 (MQ=255) tGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATAc > 2:5543894/1‑76 (MQ=255) tGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCCTAc > 2:2851423/1‑76 (MQ=255) gCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCGCCATAcc > 2:5110335/1‑76 (MQ=255) atcatcCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCCCCAAACC‑gcca > 1:7383665/1‑73 (MQ=255) tcatcCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCACCCCCACC‑GGCaa > 1:2505918/1‑76 (MQ=255) catcCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCGCCATACC‑CGGaaa > 2:6083116/1‑76 (MQ=255) atcCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTaaaa > 2:7469965/1‑76 (MQ=255) tcCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCACCATACC‑CGGaaaaa > 2:5326043/1‑76 (MQ=255) ccTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGGAAAAAc > 1:5460977/1‑76 (MQ=255) ccTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑CGTaaaaaa > 2:4663608/1‑75 (MQ=255) ccTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCCTACC‑GGGAAAAAc > 2:4482012/1‑76 (MQ=255) ccTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCACCCCTACCGGGTaaaaa > 2:6886092/1‑76 (MQ=255) tGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCATACC‑CGGAAAAAAtt > 2:4724876/1‑76 (MQ=255) gACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCCACAC‑GGCAAAAACtac > 1:695321/1‑74 (MQ=255) gACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCATACC‑CGGAAAAAAttt > 2:6078224/1‑75 (MQ=255) gACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCCCAAC‑GGTAAAAAAtta > 2:306129/1‑75 (MQ=255) gACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCACAAC‑CGGAAAAAATTc > 2:2931731/1‑76 (MQ=255) gccGGGGCCGGGACGGCCGGATTTTATGGGGAACATCACCGCGCCGTGCCCCCCCAGCCATACC‑GGTAAAAACTTc < 2:4574894/75‑1 (MQ=255) aCCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCCTACC‑GGCAAAAACATCa > 1:5026794/1‑76 (MQ=255) ggAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCATACC‑CGTAAAAACCTCAgcg > 2:119223/1‑76 (MQ=255) accGAAACGCGGGAA‑TTTATGGGAGACATCGCCGCGCCGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCg < 1:5630110/75‑1 (MQ=255) cGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCACCCATACC‑GGGAAAAACTTCAGCGCCGTa > 1:4243601/1‑76 (MQ=255) ggAACGTCGGCT‑TTTGTCGGTGACATCATCGCGCCGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTAc < 2:1522486/76‑1 (MQ=255) ggAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCACCCATACC‑GGTAAAAACTTCACCGCCGCAc > 1:741082/1‑76 (MQ=255) accGTCGGCT‑TTTAGGGGTGACACCCCCGCGCCGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCg < 1:382795/74‑1 (MQ=255) aaCGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCGCCATACC‑CGGAAAAAACTCAGCGCCGGAccc > 1:6787739/1‑75 (MQ=255) aaCGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCAAAAC‑GGCAAAAAACTCAGCGCCGTAccc > 2:1512866/1‑75 (MQ=255) aCGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAAATTCAGCGCCGTACCGa > 2:6172616/1‑76 (MQ=255) cGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCCGCCATACC‑GGTAAAAAATTCAGCGCCGGACCGac > 2:3456848/1‑76 (MQ=255) gTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCCCCACACC‑GGGAAAAACTTCACCGCCGTACCGccc > 1:497994/1‑76 (MQ=255) tCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCCACCC‑GGTAAAAACTTCACCGCCGTAACGacca > 1:1603459/1‑76 (MQ=255) tCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGGGCCCTACCGaccc > 2:3744871/1‑75 (MQ=255) tCGGCA‑TTTATCGGTGACATCATCGCGCCATGGCTCGCCCACCACACC‑GGGAAAAACTTCCGCGCCGGACCGacca > 2:6114959/1‑76 (MQ=255) cGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCCTACC‑GGTAAAAAATTCCGCGCCGTACCGcccac > 1:3390195/1‑76 (MQ=255) cGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGaccac > 2:5581408/1‑76 (MQ=255) cGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGaccac > 1:4925337/1‑76 (MQ=255) cGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGaccac > 2:5727333/1‑76 (MQ=255) cGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCACACC‑GGTAAAAACTTCAGCGCCCTACCGacccc > 2:6605069/1‑74 (MQ=255) gCA‑TTTCGCGGCGCCACCACCGCCCTGTGGCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCg < 1:7569113/76‑1 (MQ=255) cttttAGCGGTGACACCCCCGCGCCGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGa < 1:4849893/74‑1 (MQ=255) a‑TTTATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCATAAC‑GGTAAAAACTTCAGCGCCGTACCGaccacccgt > 1:2771108/1‑73 (MQ=255) ttATCGGTGACATCATCGCGCCGTGGCTCGCCCGCCCTACC‑GGGAAAAACTTCAGCGCCGTCCCGACCCCCCATCt > 2:3232052/1‑76 (MQ=255) ttATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTAACGACCACCCATCt > 1:6725135/1‑76 (MQ=255) ttATCGGTGACATCATCGCGCCGTGGCTCCCCCGCCCTACC‑GGGAAAAACTTCAGCGCCGTACCGACCCCCCATCt > 1:2841902/1‑76 (MQ=255) tATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCCACCACCGATCTg > 1:7542459/1‑76 (MQ=255) tATCGGTGACATCATCGCGCCGTGGCTCCCCAGCCATACC‑CGTAAAAACTTCAGCGCCGTACCGACCACCCATCTg > 2:4596157/1‑76 (MQ=255) tATCGGTGACATCATCGCGCCGTGGCTCCCCACCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCCATCTg > 1:6288507/1‑76 (MQ=255) tATCGGCGCCATCGTCGCGCCGTGCCCCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTg < 2:4192587/76‑1 (MQ=255) aTCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCCCCCATCTgg > 2:1284708/1‑76 (MQ=255) aTCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGGAAAAACTTCAGCGCCGTAACGACCCCCGCTCTgg > 2:179925/1‑76 (MQ=255) cGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGGAAAAACTTCAGCGCCGTACCGACCACCCATCTGGTc > 1:6323905/1‑76 (MQ=255) cGGTGACATCATCGCGCCGTGGCTCGCCAGCCATAAC‑GGTAAAAAATTCAGCGCCGTACCGACCACCGATCTGGTc > 1:3390426/1‑76 (MQ=255) ggTGCCATCCGCGCGCGGGGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCa < 2:279503/76‑1 (MQ=25) tGACATCATCGCGCCGTGGCTCCCCCGGCAAACC‑GGGAAAAACTTCCGCGCCGTACCGACCACCGATCTGGTcccc > 2:2628407/1‑76 (MQ=37) gACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGGAAAAACTTCAGCGCCGGACCCACCACCCATCTGGTCACCa > 2:283544/1‑76 (MQ=255) aCATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCaa > 1:3199966/1‑76 (MQ=255) aCATCATCGCGCCGTGGCTCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCaa > 1:2233209/1‑76 (MQ=255) atcatcGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATc > 2:795980/1‑76 (MQ=255) atcatcGCGCCGTGGCTCCCCAGCCATAAC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATc > 2:6442662/1‑76 (MQ=255) tcatcGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATcc > 2:833674/1‑76 (MQ=255) catcGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCCAACCg > 2:3248121/1‑76 (MQ=255) catcGCGCCGTGGCTCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCg > 2:2628529/1‑76 (MQ=255) catcGCGCCGTGGCTCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCg > 2:1337485/1‑76 (MQ=255) atcGCGCCGTGGCTCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGCTACCACTCCGa > 2:6533741/1‑76 (MQ=255) cgcgCCGTGGCTCGCCAGCCCTACC‑GGTAAAAACTTCAGCGCCGTAACGACCACCGATCTGGTCACCAATCCGATg > 2:7136340/1‑76 (MQ=255) cgcgCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATg > 1:4580702/1‑76 (MQ=255) cgcgCCGTGGCTCCCCCGCCATACC‑GGTAAAAACTTCAGCGCCGTACCCACCACCGATCTGGTCACCAATCCGCTg > 2:1190174/1‑76 (MQ=255) gcgcCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATgg > 1:2006478/1‑76 (MQ=255) gccccGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATgg < 1:1568517/73‑1 (MQ=38) cccGTGCCCCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGAc < 2:1453330/75‑1 (MQ=38) ccGTGGCTCCCCAGCCATACC‑GGGAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGCTGGACt > 2:1430843/1‑76 (MQ=39) cGTGGCTCGCCCGCCATACC‑GGGAAAAACTTCAGCGCCGTACCCACCCCCGATCTGGTCACCAATCCCATGGACTa > 2:1924111/1‑76 (MQ=255) cGTGGCTCCCCCGCCATACC‑GGTAAAAAATTCAGCGCCGTACCGACCACCGATCTGGTCACCCATCCGATGGACTa > 2:306437/1‑76 (MQ=38) cGTGGCTCCCCCGCCAAACC‑CGGAAAAAACTCACCGCCGCACCCCCCACCCATCTGGTCACCAATCCGATGGACTa > 2:1383874/1‑76 (MQ=255) tGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTAcc > 2:4648100/1‑76 (MQ=255) gccTCCCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCt < 2:5845113/74‑1 (MQ=255) cTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGa > 2:6279725/1‑76 (MQ=255) cTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGa > 2:1928999/1‑76 (MQ=255) tagccagccaTACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGaa > 2:5763852/3‑76 (MQ=255) tCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGaa > 1:2534476/1‑76 (MQ=255) cGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAAc > 2:2310104/1‑76 (MQ=255) | AAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCA‑TTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACC‑GGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAAC > REL606/3217359‑3217508 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |