Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,180,249 | 0 | T | G | 46.6% | 140.7 / 52.5 | 117 | G24G (GGT→GGG) | yjaH | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (14/48); new base G (53/1); total (68/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.06e-18 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 230 total aligned reads displayed.
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG > REL606/4180174‑4180324 | ggCAGTGAACAGTTGTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGt < 2:6959132/76‑1 (MQ=255) gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTg < 2:4756638/76‑1 (MQ=255) gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTg < 1:7058795/76‑1 (MQ=255) gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTg < 2:1490905/76‑1 (MQ=255) cAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGgc > 1:7217555/1‑76 (MQ=255) aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTgcg < 1:8937255/76‑1 (MQ=255) tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGgcggc > 1:538266/1‑74 (MQ=255) gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTg < 2:8996293/76‑1 (MQ=255) gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGgcggcg > 1:4604377/1‑73 (MQ=255) gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGgcggcg > 1:4331873/1‑73 (MQ=255) aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGt < 2:2346110/76‑1 (MQ=255) aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGt < 2:8070424/76‑1 (MQ=255) aCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGtt < 2:4380608/76‑1 (MQ=255) gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTtct < 2:312569/76‑1 (MQ=255) gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTtct < 2:3509420/76‑1 (MQ=255) ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTtctc < 2:1762413/76‑1 (MQ=255) ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCCGGTtctt > 2:3928957/1‑75 (MQ=255) tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCa < 1:6667410/76‑1 (MQ=255) ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCAc < 2:520241/76‑1 (MQ=255) ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCGGCTTTTCAc > 2:2749814/1‑76 (MQ=255) ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTtctccc > 1:8688121/1‑74 (MQ=255) ttAACGAAGGGGGGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGGTtctctc > 1:4978036/1‑74 (MQ=255) tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGGGCCGGTTTTCACt > 1:3617663/1‑76 (MQ=255) tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTtctccat > 1:8515104/1‑73 (MQ=255) aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTg < 2:4181538/76‑1 (MQ=255) aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCGCGTTCTCACTg > 2:8754753/1‑76 (MQ=255) aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGGCGGTTCTCACttc > 2:7040191/1‑74 (MQ=255) aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGc < 2:6831347/76‑1 (MQ=255) aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACtcc > 2:2251986/1‑74 (MQ=255) aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGc > 2:8885716/1‑76 (MQ=255) gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTTGCTGCCt > 2:5148368/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGGTGTTTCTCTCTTCCTg > 2:5430027/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGCCTg > 1:2932123/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCGGGTTCTCACTGCCTg > 2:822390/1‑76 (MQ=255) aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGt < 2:7325140/76‑1 (MQ=255) aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGt < 2:6331310/76‑1 (MQ=255) aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGt < 2:6987009/76‑1 (MQ=255) aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTTTCTGCCCGt > 1:399715/1‑76 (MQ=255) ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTa < 1:1075940/76‑1 (MQ=255) ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTa < 1:8804406/76‑1 (MQ=255) gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg < 2:2809627/76‑1 (MQ=255) ggtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGGTTCTCCCTGCCTGTAGc > 2:6256428/1‑76 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGcc > 1:7610794/1‑76 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTTCCCGTAggc > 2:3895676/1‑74 (MQ=255) ttttCCCCCGTTTTGTCTTTCTGGCGTCGATCTTTGATG‑CTGGCAGGGGCGCTGTTTCCCACTCCCTTTAGCcaca < 2:6653860/75‑1 (MQ=255) gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTTTTCCCGCCTGTCGCcaaa > 1:3177443/1‑74 (MQ=255) tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCCGTAGCCACaa > 2:6553617/1‑76 (MQ=255) tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACaa > 1:5359836/1‑76 (MQ=255) tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCCGTTTTTTACTGCCTGTAGGCCCaa > 2:3544911/1‑76 (MQ=255) ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTTGCCACAAc > 1:1901702/1‑76 (MQ=255) ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCCGGTTCTCACTGCCCGTAGCCACAAc > 2:1991458/1‑76 (MQ=255) ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGAGCCGTTTCTCACTTCCCGTGGAcaccac > 1:3856317/1‑76 (MQ=255) tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAAct < 1:6599830/76‑1 (MQ=255) cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGGTGTTTATCACAACCCGTAGCCACAAccc > 2:500096/1‑74 (MQ=255) cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACAActc > 2:2912700/1‑76 (MQ=255) cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGCTTCTCACTGGCTGTTGCCACAActc > 1:6464217/1‑76 (MQ=255) cccctTTTTGTCTTTCTGGCGTCGATCATTGAGG‑CTGGCGGGGGCGTTTTTTCTCCCTCCTTTTAGCCACAActct < 2:774908/75‑1 (MQ=255) cctttttttttttCGGGCGTCGATCTTTGATG‑CTGGCAGGTGCTTTGTTTCTCCCTGCTTTTAGCCACAACTCTTc < 2:7516481/73‑1 (MQ=255) cctttttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTc < 1:8371370/73‑1 (MQ=255) cccTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAATCTTc > 1:1924200/1‑76 (MQ=255) tcTTTTGTTTTTCTGGCTTCGTTCATTAATG‑CTGGCGGTTGCGCTTTTTCTCCCTGTTTTTAGCCACAACTCTTCa < 2:8601460/75‑1 (MQ=255) cTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACACATCTTCAc > 2:6022462/1‑76 (MQ=255) ttttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACACATCTTCAC‑t > 1:247878/1‑76 (MQ=255) ttttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGGTTTTCACCGCCTGTCGCCACAAACCCTCAC‑t > 1:2347645/1‑76 (MQ=255) ttttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCCCCGCCTGGAGCCCCAACTCTTCAC‑t > 1:1566820/1‑76 (MQ=255) tttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑tt > 1:4206198/1‑76 (MQ=255) ttGTTTTTCGGGCGTCGATCATTGATG‑CTGGCGGGTGCGTTGTTTCCCCCTGCCTTTAGCCACAACTCTTCAC‑Ttc < 2:3496317/76‑1 (MQ=255) ttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑Ttc > 1:2019671/1‑76 (MQ=255) ttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACCACTCTTCAC‑Ttc > 1:227042/1‑76 (MQ=255) ttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCCCTGCCTGTAGCCCCAACTCTTTAC‑Ttc > 2:5826432/1‑76 (MQ=255) ttGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTTAC‑Ttc > 1:8896405/1‑76 (MQ=255) gTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACAACTATTCAC‑Ttcct > 1:3348022/1‑76 (MQ=255) gTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCCCAAATCTtccc‑tcccc > 2:7182146/1‑75 (MQ=255) cTTTCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCGGTTTCTTCATGCCTGTAGACACAACTCTTCAC‑Ttcctcc > 2:2085381/1‑76 (MQ=255) cTTTCGGGCGTGGATCATTGATG‑CTGGCGGGTGCGCTGTTTCTCCCTGCCTTTAGCCACAACTCTTCAC‑Ttcctcc < 2:1063291/76‑1 (MQ=255) tttCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCCGTTTCTCACTGCCTGTAGCCACCAATCTTCAC‑TTCCTccc > 1:6588928/1‑76 (MQ=255) ttCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCCCTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCt > 1:8327566/1‑76 (MQ=255) ttCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCt > 1:574060/1‑76 (MQ=255) ttCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACCACTCTTCAC‑TTCCTCCCt > 1:314129/1‑76 (MQ=255) ttCTGGCGTCGATCATTGATG‑CTGGCAGGCGAGCGGTTTCTCACCTGCTGTTGGCACAACTCTTCAC‑TTTCTCCCt > 2:6123191/1‑76 (MQ=255) ttCTGGCGGCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCt < 2:6730160/76‑1 (MQ=255) tCTGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCtt < 2:4470970/76‑1 (MQ=255) cTGGCGTCGATCATTGATT‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCttt < 1:1461027/76‑1 (MQ=255) tGGCGTCGATCATTGATG‑CTGGCAGGTGCGGTGTTTTTCCACGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCCTTa > 1:598064/1‑76 (MQ=255) tGGCGTCGATCATTGATG‑CTGGCAGGTGCCCGGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTa < 1:6164763/76‑1 (MQ=255) tGGCGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACCGCCTGTAGCCAAAACCCTTCAC‑TTTCTCCCCTTa > 2:7891442/1‑76 (MQ=255) gggCGTCTATAA‑TGATGCCTGGCGGGGGCCCTGTTTCTCACTGCTTTTAGCCACAACTCTTCAC‑TTCCTCCCTTTa < 2:1718225/75‑1 (MQ=255) ggCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTCTCACTGCCTGTGGCCCCAACTCTTCAC‑TTCCTCCCttttc > 2:2137234/1‑74 (MQ=255) gCGTCGATCATTGATG‑CTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACAAATCTTCAC‑TTCCCCCCCTTatc > 2:1224910/1‑74 (MQ=255) cGTCGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAATCTTCAC‑TTTCTCCCTTTatcc > 2:2842065/1‑73 (MQ=255) tagATCTTTGAGT‑CGGGCGGGGGCTCTTTTTCTCCCGGCTTTTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGcc < 2:7593139/74‑1 (MQ=255) cGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCa > 1:2734702/1‑76 (MQ=255) cGATCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGACACACATCTTCAC‑TTCCTCCCTTTATCGCCa > 2:437558/1‑76 (MQ=255) gATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAg < 1:2770106/76‑1 (MQ=255) aTCATTGATG‑CTGGCAGGGGCGCTGTTTTTCACCGCCTGTTGCCACAACCCTTCAC‑TTTCTCCCTTTACCGCCAGt > 2:3560760/1‑76 (MQ=255) aTCATTGATG‑CTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCCGt > 1:3032759/1‑76 (MQ=255) aTG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTTTAGCCCACACACTTCACTTTCCTCCC‑TTACCGCCAGTGGAgttt > 2:3275170/1‑75 (MQ=255) gcgggTGCGCTGTTTCTCCCTGCTTTTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAGTGGATTTGCTGAAg < 1:1245275/73‑1 (MQ=255) ggTGCGCTGTTTCTCACTGCCTTTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc < 1:2178508/76‑1 (MQ=255) gTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa > 2:8510976/1‑76 (MQ=255) | GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATG‑CTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCAC‑TTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG > REL606/4180174‑4180324 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |