Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,822,106 | 0 | G | C | 38.8% | 240.8 / 91.1 | 129 | pseudogene (62/198 nt) | fucA | L‑fuculose‑1‑phosphate aldolase; b2800_1 |
Reads supporting (aligned to +/- strand): ref base G (78/1); new base C (0/50); total (78/51) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.79e-35 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 187 total aligned reads displayed.
TCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGA > REL606/2822031‑2822180 | gCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgcc < 2:4869797/76‑1 (MQ=255) gCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgcc > 2:360632/1‑76 (MQ=255) gCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgcc < 2:6823391/76‑1 (MQ=255) gCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgcc < 2:1246851/76‑1 (MQ=255) cTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgcca > 1:7418630/1‑76 (MQ=255) cTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGATGCGCCAGCAC‑TTCAACTTCATGCgcca > 1:2205185/1‑76 (MQ=255) tCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccag < 1:5772478/76‑1 (MQ=255) tCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccag < 2:6883142/76‑1 (MQ=255) cAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccagc < 2:4982341/76‑1 (MQ=255) cAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccagc < 2:3600405/76‑1 (MQ=255) aGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccagcc < 2:2147429/76‑1 (MQ=255) aGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccagcc > 2:4485819/1‑76 (MQ=255) aGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCgccagcc < 1:3340847/76‑1 (MQ=255) cACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCAt > 1:3628692/1‑76 (MQ=255) aCTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATa > 1:4263050/1‑76 (MQ=255) cgggCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCACCAC‑TTCAACTCCACCCCCCAGCCATaa < 2:4666312/74‑1 (MQ=255) cTGGCACCGGGTCCGTAATCGCCAGGGCCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATaa > 1:3270666/1‑76 (MQ=255) tGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAAc > 1:6628496/1‑76 (MQ=255) cACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCt > 2:3930439/1‑76 (MQ=255) cACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCt > 2:637051/1‑76 (MQ=255) cACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCAGGCCCCAGCCATAACGCt < 2:1067059/76‑1 (MQ=255) cGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCtttt > 2:2923976/1‑76 (MQ=255) gggCCCGTAAGCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCttttt < 2:555651/76‑1 (MQ=255) gTCCGTAATGGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTcc > 1:1514409/1‑76 (MQ=255) gTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTcc > 1:24068/1‑76 (MQ=255) tCCGTAATCGCCAGGGTCGTCAGGTAAAGTCGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCa < 1:5723906/76‑1 (MQ=255) ccGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAg > 1:1659720/1‑76 (MQ=255) ccGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAg > 1:4658533/1‑76 (MQ=255) ccGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTCTTCCAg > 2:4114360/1‑76 (MQ=255) cGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGa > 2:2381563/1‑76 (MQ=255) tAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTTAACTTCATGCGCCAGCCATAACGCTTTTTCCGGAtt > 1:5190863/1‑76 (MQ=255) aaTCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTc > 2:5284662/1‑76 (MQ=255) aaTCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTc > 1:3394410/1‑76 (MQ=255) accGCCGGGGTCGTCAGGTAAGGTTGCGCCAGCCC‑TTCACCTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCa < 2:4378965/74‑1 (MQ=255) tCGCCAGGGTGGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCAc < 1:6123271/76‑1 (MQ=255) tCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCAc < 1:1332577/76‑1 (MQ=255) cGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCAcc > 1:7031366/1‑76 (MQ=255) cgaGGGTCGTCATGTAAAGTTGCCCCAGCTC‑TCCAACTTCATCCCCCAGCCATAACGCTTTTTCCAGATTCACCTc < 2:6535622/74‑1 (MQ=255) ccggggTCGTCGGGTAAAGTTGCGCCAGCCC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTc < 1:4392630/73‑1 (MQ=255) ccAGGGTCTTCAGGTAAAGTGCCGCCACCCC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTc < 1:680565/76‑1 (MQ=255) cAGGGTCGTGAGGTAAAGTGGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTca < 2:7141389/76‑1 (MQ=255) cAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca > 2:2320076/1‑76 (MQ=255) cAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca > 2:3238032/1‑76 (MQ=255) aGGGTCGTCAGGTAAAGTTGCGCCAGCCC‑TCCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTcac < 1:2722406/76‑1 (MQ=255) gggTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTcaca < 1:5411178/76‑1 (MQ=255) ggTCGTCGGTTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACaa < 1:3919913/76‑1 (MQ=255) ggTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCAGGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACaa < 1:6222077/76‑1 (MQ=255) ggTCGTCAGGAAAAGTTGCCCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACaa < 1:1201449/76‑1 (MQ=255) gtcgtcAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAg < 2:4489446/76‑1 (MQ=255) tcgtcAGTTAAAGTGGCGCCAGCAC‑TTCAACTTCACGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGc < 1:4784473/76‑1 (MQ=255) cgtcAGGTAAAGTTGCCCCACCCC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCg < 1:5943878/76‑1 (MQ=255) tcAGGTAAAGTTGCGCCAGCAC‑TTCACCTTCAGGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGAt < 1:3383724/76‑1 (MQ=255) tcAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGAt > 2:4257768/1‑76 (MQ=255) taaGGAAAAGTTGCCCCACCAC‑TCAACCTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGAt < 1:6762676/74‑1 (MQ=255) cAGGTAAAGTTGCGCCAGCCC‑TTCAACTTCATCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATa < 2:598562/76‑1 (MQ=255) cAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATa > 2:6736776/1‑76 (MQ=255) cAGGTAAAGTTGCCCCACCCC‑TTCACCTTCATCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATa < 1:6283399/76‑1 (MQ=255) aGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATaa > 2:5952395/1‑76 (MQ=255) aGGTAAAGTGCCCCCAGCAC‑TTCACCTCCACCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATaa < 1:1154201/76‑1 (MQ=255) aGGTAAAGGTGCGCCAGCCC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATaa < 1:1790075/76‑1 (MQ=255) ggTAAAGTCGCCCCAGCAC‑TTCACCTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAg < 1:5245602/76‑1 (MQ=255) gTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGc > 2:343734/1‑76 (MQ=255) aaaaGGTGCCCCACCAC‑TCCACCTTCAGCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGcc < 2:617614/75‑1 (MQ=255) aaaGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc > 1:1052233/1‑76 (MQ=255) aaaGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc > 2:4997191/1‑76 (MQ=255) aaaGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc > 2:6666900/1‑76 (MQ=255) aaaGTTGCGCCAGCAC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc < 1:1879040/76‑1 (MQ=255) aaaGTTGCCCCACCAC‑TTCACCTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc < 2:2233586/76‑1 (MQ=255) aaaGTGGCCCCAGCCC‑TCAAACTCCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc < 2:5082210/76‑1 (MQ=255) aaGTTGCGCCAGCAC‑TTCGACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCa > 1:1267522/1‑76 (MQ=255) aaGTTGCCCCCGCACTTTCACCTTCAGCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGccc < 2:3821727/76‑1 (MQ=255) aGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCat > 1:4461305/1‑76 (MQ=255) aGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCat > 1:5424516/1‑76 (MQ=255) aGTTGCCCCAGCCC‑TTCAACTTCATGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCat < 1:5636934/76‑1 (MQ=255) gTTCCGCCAGCAT‑TCCACCTCCACCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatg < 1:6837057/76‑1 (MQ=255) ttGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatga > 1:2774081/1‑76 (MQ=255) ttGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatga > 2:85955/1‑76 (MQ=255) gtGCCCCACCAC‑TTCACCTCCAGGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatga < 2:1434660/75‑1 (MQ=255) tGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatgat > 1:5872131/1‑76 (MQ=255) gcgcCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCatgatg > 2:771673/1‑76 (MQ=255) cAGCTC‑TTCATCTTCACCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGAtgttgt < 2:6053981/76‑1 (MQ=255) cAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGAtgttgt > 2:1583348/1‑76 (MQ=255) cAGCAC‑TCCACCTTCACCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGAtgttgt < 1:4153027/76‑1 (MQ=255) gCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTaa > 1:2262672/1‑76 (MQ=255) cAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAAc > 2:1366594/1‑76 (MQ=255) tc‑ttcaccatcaGGCGCCAGCCATAACGCTTTTTCCAGATCCACCTCACAAGCGATAAGCCCATGATGTTGTAACa < 2:113146/75‑1 (MQ=255) ctctcaacttcaGCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACa < 1:2983503/73‑1 (MQ=255) ttcaactccaGCCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACaaa < 2:3544986/76‑1 (MQ=255) aacttcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGtt > 2:1900042/1‑76 (MQ=255) ccttcaTGCCCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTg < 2:3307019/75‑1 (MQ=255) cttcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGc > 1:5538143/1‑76 (MQ=255) cttcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGc > 1:3838397/1‑76 (MQ=255) ttcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGcc > 2:6699750/1‑76 (MQ=255) ttcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGcc > 1:6413844/1‑76 (MQ=255) tcaTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCt > 2:3128719/1‑76 (MQ=255) aTGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTa > 2:1927068/1‑76 (MQ=255) gcgcCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACg > 1:3986421/1‑76 (MQ=255) ccgcCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACg < 2:7292353/75‑1 (MQ=255) cgcCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGa > 1:7112584/1‑76 (MQ=255) cgcCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGa > 1:458161/1‑76 (MQ=255) | TCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCAC‑TTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGA > REL606/2822031‑2822180 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |