Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,061,648 | 0 | A | C | 37.2% | 282.3 / 39.2 | 145 | G40G (GGT→GGG) | yccY | phosphotyrosine‑protein phosphatase |
Reads supporting (aligned to +/- strand): ref base A (79/12); new base C (0/54); total (79/66) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.91e-28 | |||||||||||
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 198 total aligned reads displayed.
TTCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAATGTT > REL606/1061573‑1061722 | ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAaca < 1:2347292/76‑1 (MQ=255) ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAaca < 1:7767184/76‑1 (MQ=255) ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAaca < 1:6796192/76‑1 (MQ=255) ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAaca < 1:5686982/76‑1 (MQ=255) ttCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAaca > 1:3117906/1‑76 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCCT‑GAaaac > 1:6894577/1‑76 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAacac > 1:332790/1‑76 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAacac < 2:3178290/76‑1 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAacac < 1:7942514/76‑1 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAacac < 1:3263212/76‑1 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAacac < 1:5395100/76‑1 (MQ=255) tCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGGTTTACC‑AGACCAT‑GAacac > 2:847529/1‑76 (MQ=255) cAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGGATATGCAGGGTGGTTTAAC‑AGACCAT‑GAACAcc > 2:1064783/1‑76 (MQ=255) cAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACAcc < 1:7004958/76‑1 (MQ=255) cAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACAcc > 1:4566016/1‑76 (MQ=255) aaGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACAccc > 2:8006246/1‑75 (MQ=255) aaGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCt < 1:4417987/76‑1 (MQ=255) gagaTACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGTTTTACC‑AGACCAT‑GAACACCTGc > 2:8056964/1‑76 (MQ=255) agaTACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGcc < 1:8200540/76‑1 (MQ=255) agaTACGCCATGATTAGCGGCAACATCTGCTGCGGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACAACTGcc > 1:5901850/1‑76 (MQ=255) aCGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGTTTTACC‑AGACCAT‑GAACACCCGCCGAtt > 2:2854444/1‑76 (MQ=255) gCCATGGTTAGCGGCAACCTTTGCTGCGGTCGCATCTGCAGGGTGTTTTACC‑AGCCCAT‑GAACCCCTGCCGATTTg < 2:4112472/76‑1 (MQ=255) ccATGTTTGGCGCCAACTTCTGCTGCTGTCGCATCTGCAGGGTGTTTTCCC‑AGCCCAT‑AACCCCCTGCCGATTTGa < 2:1587643/76‑1 (MQ=255) ccATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGGTTTACC‑CGACCAT‑GAACACCTGCCGGTTTGa > 2:2141949/1‑76 (MQ=255) aTGATTGGGGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAG‑GAACCCCTGCCGATTTGACt < 2:5178328/76‑1 (MQ=255) aTGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCCGGGGGTTTTACC‑AGACCAT‑GAACACCTGCCGCTTTGACt > 1:8201917/1‑76 (MQ=255) aTGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACt < 2:2155939/76‑1 (MQ=255) tGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGGTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACtt > 2:6848213/1‑76 (MQ=255) aTTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGTTTTACC‑AGACCCT‑GAACACCTGCCGCTTTTACtttt > 2:6578423/1‑76 (MQ=255) ttAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCCT‑GAACACCTGCCGATTTGACTTTTa > 2:4414735/1‑76 (MQ=255) ttAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTa < 1:4235529/76‑1 (MQ=255) ttAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTa < 2:176697/76‑1 (MQ=255) ttAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTa < 1:2787690/76‑1 (MQ=255) ttAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTa < 2:5067355/76‑1 (MQ=255) ttAGCGGCAACATCGGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTa < 1:2432573/76‑1 (MQ=255) tAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTAc > 1:1036302/1‑76 (MQ=255) aGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGGTTTGACTTTTACg > 2:2405641/1‑76 (MQ=255) aGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGGGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACg > 2:204796/1‑76 (MQ=255) aGCGCCACCATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACg < 2:4705478/76‑1 (MQ=255) gCGGCAACATCTGCTTCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGc < 1:4113206/76‑1 (MQ=255) gCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGc < 2:3356278/76‑1 (MQ=255) cGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGcc > 1:1027956/1‑76 (MQ=255) gccAACATCTGCTGCTGCCGCATCTGCGGGGTGTTTTCCC‑AGCCCAT‑GAACCCCTGCCGATTTGACTTTTACGCCt < 1:4996554/74‑1 (MQ=255) gCAACATCGGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTg < 1:5270399/76‑1 (MQ=255) cAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTgg < 1:3549013/76‑1 (MQ=255) aaCATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGc < 2:5680808/76‑1 (MQ=255) aaCATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGc > 1:317968/1‑76 (MQ=255) aTCTGCTGCTGTCGCATCTGCGGGGTTTTTTACC‑AGCCCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGa < 2:2508186/76‑1 (MQ=255) aTCTGCTGCTGTCGCATCTGCAGGGGGTTTTACC‑AGACCAT‑GAACAACTGCCGATTTGACTTTTACGCCTGGCAGa > 2:5868501/1‑76 (MQ=255) ttggctgctgTCGCCTTTGCAGTGTGTTTTACC‑AGACTAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGAc < 2:6620810/73‑1 (MQ=255) tCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGAc > 1:6744909/1‑76 (MQ=255) tCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGAc > 1:6217826/1‑76 (MQ=255) ctgctgctgTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACg > 1:7311618/1‑76 (MQ=255) gctgctgTCGCATCTGCAGGGTTTTTTACC‑AGCCCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGtt < 2:304669/76‑1 (MQ=255) gctgctgTCGCATCTGCAGGGGGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGtt > 1:5139352/1‑76 (MQ=255) ctgctgTCGCATCTGCGGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGttt < 1:5581546/76‑1 (MQ=255) cgtctgTCGCATCGGCAGGGTGTTTTACC‑AGCCCAG‑AAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGttt < 1:5741893/73‑1 (MQ=255) ctgTCGCATCTGCAGGGGGTTTTAAC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTAGGCCTGGCAGACGTCTACg > 2:3675301/1‑76 (MQ=255) tgTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACgc < 1:6125451/76‑1 (MQ=255) gccGCATCTGCGGGTTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACgca < 1:6532301/74‑1 (MQ=255) tCGCATCTGCAGGTTGTTTTACC‑AGCCCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACgcag < 2:3546462/76‑1 (MQ=255) tCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACgcag > 2:5905932/1‑76 (MQ=255) cGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACgcagc > 1:3058548/1‑76 (MQ=255) cATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCaa < 2:6159408/76‑1 (MQ=255) cATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCaa > 1:2650987/1‑76 (MQ=255) cATCTGCA‑GGTGTTTTTCC‑AAGCAATCGACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCaa < 2:6499283/76‑1 (MQ=255) aTCTGCAGGGTGTTTTACC‑AGACCAT‑GACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAg < 1:216877/76‑1 (MQ=255) tCTGCAGGGTTTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAgc < 2:7373676/76‑1 (MQ=255) tCTGCAGGGTGTTTTACC‑AGACCAT‑GACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAgc < 2:5434039/76‑1 (MQ=255) tCTGCAGGGTGTTTTACC‑AGACCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAgg < 1:2742606/76‑2 (MQ=255) tGCAGGGTGTTTTACC‑AGCCCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAgcgc < 2:5262518/76‑1 (MQ=255) gCAGGTTGTTTTACC‑AGACCAT‑GACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCt < 2:7092249/76‑1 (MQ=255) gCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCt > 1:4726130/1‑76 (MQ=255) cAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTc > 1:186258/1‑76 (MQ=255) cAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTc > 1:2827584/1‑76 (MQ=255) aGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCg > 1:4864490/1‑76 (MQ=255) gtgtTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCa > 1:1492038/1‑76 (MQ=255) ttttttCCC‑AATCCAT‑AACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCaa < 1:5442618/74‑1 (MQ=255) tgtTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGTGATCGCCaa > 1:8077541/1‑76 (MQ=255) tgtTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCaa > 2:5619956/1‑76 (MQ=255) tttttACC‑CGACCAT‑ACCCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAAt < 2:2631151/75‑1 (MQ=255) ttttACC‑AGACCAT‑AACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATg < 1:3711942/76‑1 (MQ=255) tttACCAAGACAAT‑GACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATg < 1:4655178/76‑1 (MQ=255) tttACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATgg > 1:1584887/1‑76 (MQ=255) ttACC‑AGACAAT‑AAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATggg < 2:4264785/76‑1 (MQ=255) tACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGc > 1:7250496/1‑76 (MQ=255) ccc‑AGACCAT‑GACCCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCg < 2:1475771/75‑1 (MQ=255) aCC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCg > 1:597225/1‑76 (MQ=255) c‑AGCCCAT‑GAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGaa < 2:7271101/76‑1 (MQ=255) c‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGaa > 2:4569187/1‑76 (MQ=255) c‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGaa > 2:3842325/1‑76 (MQ=255) ccAT‑AAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGc < 2:1819047/76‑1 (MQ=255) cAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCa > 2:1004262/1‑76 (MQ=255) cAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCa > 2:7417985/1‑76 (MQ=255) t‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCaaa > 1:8268621/1‑76 (MQ=255) gAACCCCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAAt < 1:386434/76‑1 (MQ=255) gAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAAt > 2:6662085/1‑76 (MQ=255) aaCACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAATg > 1:7432052/1‑76 (MQ=255) aaCACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAATg > 2:1127228/1‑76 (MQ=255) acacCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAATGt > 2:1402144/1‑76 (MQ=255) | TTCAAGAGATACGCCATGATTAGCGGCAACATCTGCTGCTGTCGCATCTGCAGGGTGTTTTACC‑AGACCAT‑GAACACCTGCCGATTTGACTTTTACGCCTGGCAGACGTTTACGCAGCAAGCGCTCGCCAATGGGCGAACGGCAAATGTT > REL606/1061573‑1061722 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |