Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,029,711 | 0 | G | C | 24.1% | 73.5 / 10.2 | 29 | R207P (CGC→CCC) | pqiA | paraquat‑inducible membrane protein A |
Reads supporting (aligned to +/- strand): ref base G (21/1); new base C (0/7); total (21/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.13e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTC > REL606/1029636‑1029785 | tCCTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTgcg > 2:915075/1‑76 (MQ=255) tCCTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTgcg > 1:54196/1‑76 (MQ=255) cTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGcc > 1:976791/1‑76 (MQ=255) cTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGcc < 1:2061193/76‑1 (MQ=255) tGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCa > 1:22807/1‑76 (MQ=255) gCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAg > 1:1761975/1‑76 (MQ=255) tGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCt > 2:1998653/1‑76 (MQ=255) ttttCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAAccg > 1:1508581/1‑75 (MQ=255) gTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTc > 2:2061916/1‑76 (MQ=255) gCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCAc > 2:283019/1‑76 (MQ=255) ttGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACggg > 2:152057/1‑76 (MQ=255) ttGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGCCACggg > 2:731361/1‑76 (MQ=255) tgttcgtggCTGTTTAGGGGCCGCCCCCCCCCCGAGGCCAGACCGCCCCCAGCCGCTAAAACCAGGCGTCACgggg < 1:331801/73‑1 (MQ=255) tGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACgggg > 1:900874/1‑76 (MQ=255) tGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACgggg > 1:1559439/1‑76 (MQ=255) aTCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGAt > 2:1646365/1‑76 (MQ=255) cgtggCTGTTTAGGGGACGACACCCCCCCGATGCCAGACCCGCCCCAGCCGCTAAAACCAGGCGTCACGGGGATAc < 2:588923/76‑1 (MQ=255) cgtcgtTGTTTAGGGGCCGCCCCCCCCCCGACGCCAGACCGGCCCCAGCCGCTAAAACCAGGCGTCACGGGGATAc < 2:828239/76‑1 (MQ=255) cgtcgCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATAc > 1:2031640/1‑76 (MQ=255) tGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAgg > 1:165118/1‑76 (MQ=255) ttAGGGGCCGACACCGCCCCGATGCCAGACCCCCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGc < 1:1596995/76‑1 (MQ=255) tATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAggggag > 1:338989/1‑74 (MQ=255) aTGGGCCGCCACCGCCCCAAGGCCAGACCGGCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGNTg < 2:1462718/76‑1 (MQ=255) aTGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTg > 1:1809162/1‑76 (MQ=255) acgacCCCGCCCCGCCGCCACACCCGCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTc < 1:1800887/76‑1 (MQ=255) acgacATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTc > 2:1624068/1‑76 (MQ=255) acgacATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGCCACGGGGGTACGTCAGGGGCTGCGTTc > 1:863681/1‑76 (MQ=255) gacATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGttagg > 2:851092/1‑73 (MQ=255) cccccGACGCCAGACCCCCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCt < 1:1301276/75‑1 (MQ=255) cccccGACCCCAGACCGCCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCt < 2:1122199/75‑1 (MQ=255) ccccGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCtg > 2:1639954/1‑76 (MQ=255) ccccGACCCCAAACCCGCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCtg < 2:593551/76‑1 (MQ=255) ccGATGCCAGACCCCCCCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCtgtt < 1:1967974/76‑1 (MQ=255) tGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTAc > 2:1222411/1‑76 (MQ=255) gAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCaa > 2:1070237/1‑76 (MQ=255) gAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCaa > 2:871613/1‑76 (MQ=255) cgcCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAAtccttc > 1:1430020/1‑76 (MQ=255) | TCCTGCAACTGCGCGCTTTTCAGTGCGTTGATCGTCGCTGGTTATGGGACGACATCGCCCCGATGCCAGAACTGCGCCAGCCGCTAAAACCAGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTC > REL606/1029636‑1029785 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |