| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | 76 | 21,532 | 0 | T | C | 65.7% | 73.8 / 62.7 | 70 | intergenic (–/–) | –/– | –/– |
| Reads supporting (aligned to +/- strand): ref base T (14/10); new base C (19/27); total (33/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.12e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGATGAGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > 76/21434‑21559 | ctGATGAGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCttcttc < 2:1463152/101‑1 (MQ=255)ctGATGAGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCttcttc < 2:584512/101‑1 (MQ=255) gatgaGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCgg > 1:234101/1‑101 (MQ=255) gatgaGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCgg > 2:1209792/1‑101 (MQ=255) gatgaGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCgg < 2:869675/101‑1 (MQ=255) atgaCTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGCCGAACTTCATCATGACCCTTTCCTTATACACGGCCCCTGCTTCggc > 2:274767/1‑100 (MQ=255) tgaGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGa > 2:1598059/1‑101 (MQ=255) tgaGTTGTTAGGAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCcttctggttta < 2:1255219/101‑7 (MQ=255) gaGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAg < 2:1158084/101‑1 (MQ=255) aGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAgg < 2:1361879/101‑1 (MQ=255) aGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAgg < 2:1083849/101‑1 (MQ=255) tgttAGTAAGTCGGGAAACTTAACAGTAACAACACAATTTGTTTAATGACACGCCTCCCCCCAACCCTTTCCTTCTACACGGCCCCTTCGCCGGGacgatc < 2:1408655/101‑6 (MQ=255) ttAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCtt > 2:907959/1‑101 (MQ=255) tAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCttt < 2:412122/101‑1 (MQ=255) tAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCttt > 1:1136786/1‑101 (MQ=255) aGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTc > 2:1719859/1‑101 (MQ=255) aGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTc < 2:1502279/101‑1 (MQ=255) tAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTccc > 2:1030798/1‑101 (MQ=255) aaGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCg < 2:269377/101‑1 (MQ=255) aaGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCg > 2:676371/1‑101 (MQ=255) aGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGt > 2:70361/1‑101 (MQ=255) cGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTc < 1:1686587/101‑1 (MQ=255) cGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTc > 1:1364062/1‑101 (MQ=255) gggAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCa > 1:1374889/1‑101 (MQ=255) ggAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAg > 2:55625/1‑101 (MQ=255) ttAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGcccctt < 1:190119/101‑4 (MQ=255) aaCAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGccccttgc > 1:686808/1‑96 (MQ=255) gTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGccccttgctgcc > 1:10467/1‑92 (MQ=255) tAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGccccttgctgccc < 1:834048/101‑11 (MQ=255) ccacgatgatgatccaCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggctta > 1:292403‑M2/6‑78 (MQ=255) ccacgatgatgagtatCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttc < 1:1666818‑M2/96‑24 (MQ=255) ccacgatgatgagtatCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggctta > 1:1284430‑M2/6‑78 (MQ=255) cacgatgatgatccACTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttac < 1:24991‑M2/97‑25 (MQ=255) cacgatgatgagtaTCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttac > 2:684277‑M2/5‑77 (MQ=255) cacgatgatgagtaGCTTCATCATGACCCTTTCCTTATTTAAGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcttcggcttac > 1:946408‑M2/5‑77 (MQ=255) gatgatgaTCCACTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctc < 1:529865‑M2/100‑28 (MQ=255) atgatgAGTATCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctct < 1:1580890‑M2/101‑29 (MQ=255) atgatgAGTAGCTTCATCATGACCCTTTCCTTATTTAAGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcttcggcttacctcc < 2:172066‑M2/101‑29 (MQ=255) tgatgaTCCACTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctcct < 2:1605876‑M2/101‑30 (MQ=255) gatgaGTAGCTTCATCATGACCCTTTCCTTATTTAAGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcttcggcttacctcctt < 2:1405924‑M2/101‑31 (MQ=255) atgagtatcttcatcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaattgtcggcttacctctgta < 2:498949‑M2/93‑32 (MQ=255) tgagtatcttcatcaTGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtac > 2:652080‑M2/8‑69 (MQ=255) agtatcttcatCATGACCCTTGCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactg < 2:1713100‑M2/96‑35 (MQ=255) tatcttcATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgca < 1:847502‑M2/98‑37 (MQ=255) tatcttcATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgca < 1:346112‑M2/98‑37 (MQ=255) atcttCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcac < 1:139896‑M2/99‑38 (MQ=255) catcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccatt > 1:659568‑M2/1‑59 (MQ=255) catGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttg > 2:1248336‑M2/1‑56 (MQ=255) atGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgt < 2:1580309‑M2/101‑47 (MQ=255) tgaccCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgtt > 2:1652922‑M2/3‑54 (MQ=255) tgaccCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgtt < 2:1471531‑M2/99‑48 (MQ=255) tgaccCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgtt > 2:1086083‑M2/3‑54 (MQ=255) gacCCTTTCCTTATTTATGGCCCCTTCCTCGGGGGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttg < 2:1720266‑M2/100‑49 (MQ=255) aCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttgg < 2:850908‑M2/101‑50 (MQ=255) cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttgga < 1:329850‑M2/101‑51 (MQ=255) cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttgga < 2:784123‑M2/101‑51 (MQ=255) cccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttgga > 2:876815‑M2/1‑51 (MQ=255) ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggac > 1:1594825‑M2/1‑50 (MQ=255) ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggac < 1:1526846‑M2/101‑52 (MQ=255) ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggac < 1:970003‑M2/101‑52 (MQ=255) cTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttccctctgtactgcacccattttgttggacg > 1:1605882‑M2/1‑49 (MQ=255) cTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacg < 1:597339‑M2/101‑53 (MQ=255) tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacga < 2:862067‑M2/101‑54 (MQ=255) ttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgat > 2:929375‑M2/1‑47 (MQ=255) tCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatg > 1:356805‑M2/1‑46 (MQ=255) tCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggctcacctctgtactgcacccattttgttggacgatg > 1:777913‑M2/1‑46 (MQ=255) ttatttatggccccttccTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaa < 2:1078956‑M2/93‑59 (MQ=255) tatttatggccccttcCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaat > 2:1257467‑M2/8‑42 (MQ=255) atttatggccccttCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatg < 2:195018‑M2/95‑61 (MQ=255) atttatggccccttCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatg > 2:1642514‑M2/7‑41 (MQ=255) tttatggcccctTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatgg < 2:1717021‑M2/96‑62 (MQ=255) ntggccCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaat < 1:373776‑M2/99‑65 (MQ=255) atggccCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaat > 1:1673798‑M2/3‑37 (MQ=255) atggccCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaat > 1:1224312‑M2/3‑37 (MQ=255) ggCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatag < 2:1089572‑M2/101‑67 (MQ=255) ggCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatag < 1:876847‑M2/101‑67 (MQ=255) cccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatagccc > 1:1580315‑M2/1‑32 (MQ=255) ccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatagcccc > 1:670435‑M2/1‑31 (MQ=255) cttcctCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatagcccct > 1:1097645‑M2/1‑30 (MQ=255) cttcctCGGGAGGGGCTTTCCCGTTTCAGCgtcccgctgaaatcgtcggcttacctctgtactgcacccattttgttggacgatgaaatggaatagcccct < 1:61235‑M2/101‑72 (MQ=255) | CTGATGAGTTGTTAGTAAGTCGGGAAACTTAACAGTAACAACACAACCAGTATGATGACGAACTTCATCATAACCCTTTCCTTATACAAGGCCCCTTCTTCGGGAGGGGCTTTCCCGTTTCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ > 76/21434‑21559 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |