| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,348,325 | 0 | A | C | 27.8% | 118.4 / 21.0 | 54 | G113G (GGT→GGG) | nuoA | NADH dehydrogenase subunit A |
| Reads supporting (aligned to +/- strand): ref base A (27/12); new base C (0/15); total (27/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.02e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGCCTACCC > REL606/2348275‑2348374 | gAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGa > 1:5142188/1‑51 (MQ=255)gAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGa < 1:512463/51‑1 (MQ=255)gAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGa < 1:6798819/51‑1 (MQ=255)gAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGa > 1:6232503/1‑51 (MQ=255) aaCGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGAc > 1:1098634/1‑51 (MQ=255) aCGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGAcc > 1:2727818/1‑51 (MQ=255) cgcgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCt < 1:472835/51‑1 (MQ=255) gcgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTg < 1:5398982/51‑1 (MQ=255) gcgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTg < 1:5411935/51‑1 (MQ=255) gcgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTg < 1:2746509/51‑1 (MQ=255) gcgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTg > 1:4896949/1‑51 (MQ=255) cgcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGc < 1:7295907/51‑1 (MQ=255) gcgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGcc > 1:6441283/1‑51 (MQ=255) cgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGa < 1:3878828/44‑1 (MQ=255) cgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCa < 1:5088394/51‑1 (MQ=255) cgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCa < 1:3573412/51‑1 (MQ=255) cgGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCa > 1:2469675/1‑51 (MQ=255) gggCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCa < 1:270895/50‑1 (MQ=255) gggCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCAg < 1:475020/51‑1 (MQ=255) gggCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCAg < 1:1872816/51‑1 (MQ=255) cGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTaa > 1:4789961/1‑51 (MQ=255) tccagtccagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAaca > 1:3751032/1‑51 (MQ=255) ccagtccagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAacac > 1:7144085/1‑51 (MQ=255) gtccagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaa > 1:5646675/1‑51 (MQ=255) tccagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaa > 1:5276932/1‑51 (MQ=255) cagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaaat > 1:777380/1‑51 (MQ=255) cagCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaaat > 1:3136911/1‑51 (MQ=255) gcgcCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaaataaaa > 1:3271941/1‑51 (MQ=255) gcgcCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaaataaaa > 1:4303833/1‑51 (MQ=255) gcgcCAATACGCACCAGATAAACCAGACCTGCCAGTAACACaaaaataaaa > 1:2296665/1‑51 (MQ=255) ccAATACGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAAtt > 1:6046670/1‑51 (MQ=255) cAATACGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTg > 1:6208699/1‑51 (MQ=255) aaTACGCACCAGATAAACCAGCCCTGCCAGTAACACAAAAATAAAAATTg < 1:6192942/50‑1 (MQ=255) aaTACGCACCAGAAAAACCAAACCTGCCAGTAACACAAAAATAAAAAtt < 1:3184885/49‑1 (MQ=39) tACGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAg > 1:2687671/1‑51 (MQ=255) cGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCt > 1:6870995/1‑51 (MQ=255) gCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCtt > 1:1273943/1‑51 (MQ=255) gCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCtt > 1:4706586/1‑51 (MQ=255) cACCAGATAAACCAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCtt < 1:3742814/50‑1 (MQ=255) cACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTc > 1:2087640/1‑51 (MQ=255) aGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCcaca > 1:683277/1‑51 (MQ=255) aaaaaCCAGCCCTGCCAGTAACACaaaaataaaaat < 1:661824/34‑1 (MQ=37) aaaaaCCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTcc < 1:744524/44‑1 (MQ=255) tAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAg > 1:2729410/1‑51 (MQ=255) aacccAGCCCTGCCAGTAACACaaaaataaa < 1:6254771/28‑1 (MQ=25) aacccAGCCCTGCCAGTAACACAAAAATAAAAATTGc < 1:872601/34‑1 (MQ=37) aacccAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCtt < 1:4515242/39‑1 (MQ=38) aacccAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCtt < 1:4142288/39‑1 (MQ=38) aacccAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACaaa < 1:5617382/46‑1 (MQ=39) aacccAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGc < 1:7280012/48‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCa < 1:2802604/38‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCa < 1:5421122/38‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCAg < 1:2696224/39‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGc < 1:887963/40‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCtt < 1:6515699/42‑1 (MQ=255) aaaCCAGCCCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCca < 1:3141487/45‑1 (MQ=255) ccATACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGCCTa > 1:3800924/1‑51 (MQ=255) ccAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGCCTa > 1:1298988/1‑51 (MQ=255) gACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGCCTAccc > 1:2962145/1‑51 (MQ=255) | GAACGCGCGGGCGTCCAGTCCAGCGCGCCAATACGCACCAGATAAACCAGACCTGCCAGTAACACAAAAATAAAAATTGCAGCTTCCACAAAGCCTACCC > REL606/2348275‑2348374 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |