Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 2,835,167 | 0 | C | G | 30.4% | 51.7 / 19.3 | 47 | R171T (AGA→ACA) | CJ019_02688 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (15/17); new base G (13/1); total (28/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.44e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACTTCGCTTTGTAGCAGAAGGA > ACB122/2835035‑2835299 | tCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACt > 4:149408/1‑137 (MQ=255) ccTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACtt < 2:287574/137‑1 (MQ=255) cccGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTTGGtt > 3:68685/1‑137 (MQ=255) ccGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTtg < 5:5038/137‑1 (MQ=255) acGAGACAAAGTCTTTGCCGCCAAGCAGAGTAATCTCCGCTCGATGTCGAGGGTTGTTGATAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTATACTTTCGGTtg < 6:79701/136‑1 (MQ=255) gAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTtgtg < 1:5651/137‑1 (MQ=255) cAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTtg > 1:108449/1‑131 (MQ=255) cAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTtg < 2:108449/131‑1 (MQ=255) tGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACa < 1:203232/137‑1 (MQ=255) gcagcGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGGCAGTTCGGTTGTGGGTAATCTATACAACGCTGttt > 6:7147/1‑137 (MQ=255) gcagcGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGttt < 2:299296/137‑1 (MQ=255) agcGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACGTTAGGTTGTGGTTTATCTATACAACGCTGTTTCg > 8:204713/1‑137 (MQ=255) cGTAAGCTCCGCGCGATGTCGAGTGTTGTTGCTAGCCACTGATGTTAGTGCATCGGGAGTCATAGGTTTACCGTATGAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATAGAACGCTGTTTCGAt < 2:328064/137‑1 (MQ=255) aaGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCTTGTTCATTGGGTTGTGGTGTATCTTTACAACGCGGTTTAGATGGt > 5:2015/1‑137 (MQ=255) cgcgATGTCGAGTGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATATGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATg < 7:322486/137‑1 (MQ=255) cgcgATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATg < 2:11266/137‑1 (MQ=255) cgATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGaa > 5:151791/1‑137 (MQ=255) cgATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTAGATGGTTGCTAATGaa > 3:115565/1‑137 (MQ=255) cgATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTAATTTAGGTTGTTGAATATATATAAAAAGCTGTTTAGAAGGTTTCTAATGaa > 5:167684/1‑137 (MQ=255) gATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCGACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGGTAATGaat > 3:244983/1‑136 (MQ=255) tgCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGGTAATCTATAAAACGCTGTTTAGATGGTTGCTAATGAATCAGCAGAAGCTTAAt > 8:49791/1‑137 (MQ=255) cTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTAATTTCGGTTGTGGTTAATCTATAAAAAGCTGTTTAGATGGTTTATAATAAATCAGCAGAATCTTAAttt > 1:259194/1‑136 (MQ=255) cTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTAGGGTGTGGGTAATATATACAACCCTGTTTCGATGGTTGCTAATTAAGCAGCAGAAGCTTCAttt > 3:303428/1‑136 (MQ=255) cTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGGCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTg < 8:95750/137‑1 (MQ=255) cACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGTTGTGGATTATATATAAAAAGCTGTTTCGATGGTTGCTAATTAAGCAGCAGAAGCTTCATTGGCacc > 5:37352/1‑135 (MQ=255) aCTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGTTGTGGTTAATCTATACAACGCTGTTTAGATGGTTGCTAATTAAGCAGCAGAAGCTTCATTGGCAgcg > 7:108840/1‑137 (MQ=255) tgatAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGGACTGTAGGTTGTGGTTATTATATACAACACTGTTTGGTTGGTTGCTAATACATCAGCCGAAGATTCATTTGAAgcgaaagg > 2:162910/1‑132 (MQ=255) atAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGGGCTTTAGGTTGGGGTTAATATATTCAAAGCTGTTTCGATGGTTGGTAATGAAGCAGCAGCAGCGTAATTGGCCGCGCATCGAc > 4:100173/1‑137 (MQ=255) gCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACGTTGGGTTGTTGGTTATATATAAAAAGCTGTTTGGATTGTTTCTAATTAAGCAGCAGAAGCTTAATTGGCAGCGCATCGACTGATg > 1:49915/1‑137 (MQ=255) aTCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCa < 6:151791/137‑1 (MQ=255) gggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACtt > 4:37122/1‑48 (MQ=255) gggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGGTGTGGATTATCTATAAAACGCTGTTTAGATGGTTGGTAATTAAGCAGCAGAAGCTTCATTTGAAGCGCATAGACTGATGCAGcc > 5:270429/1‑137 (MQ=255) gggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACtt < 3:37122/48‑1 (MQ=255) gggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATTAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGcc > 5:332007/1‑137 (MQ=255) gggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGcc > 2:35707/1‑137 (MQ=255) ggAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGTTGTGGTTTATCTATAAAACGCGGTTTAGATGGTTGGTAATGTAGCAGCAGAAGCTTCATTTGCAGCGCATCGACTGATGGAGCCa > 2:298412/1‑137 (MQ=255) gAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGGTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCaa > 7:213465/1‑137 (MQ=255) ttcttAGGTTTACCGTATCTATAAATTAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGGAGCGCATCGACTGATGCAGCCAAGc < 5:188440/133‑1 (MQ=255) tAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACGTTGGGTTGTGGATAATCTATACAACCCTGTTTGGATGCTTGCTAAATAAGCAGCAGAATCTTCATTTGCAGCGCATCGACTGAATCAGCCAAGCGttt > 2:6784/1‑136 (MQ=255) ggTTTACCGTATCAAAAAAATAGACCGGTCGTGTACTTTCGGGTGTGGATAATATATACAACGCAGTTTCGATGGTTGCTAATTAAACAGCAGAAGCTTCATTGGCAGCGCATCGACGTCTGCACCCAAGCGTTAcc > 8:228880/1‑137 (MQ=255) ggTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGTGTTAcc < 2:49915/137‑1 (MQ=255) cGTATCAAAAAAATAGACCGGTCGTCGAGTTTAGGTTGTGGATTATTTATAAAACCCGGTTTCGTTGGTTTGTAATTTATCATCAGAAGCGTCAATGTCAGCGCAGAGCCTGAGGCAGCCAAGCCTTAACACGCggg > 4:3867/1‑136 (MQ=255) tATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACt < 6:37352/137‑1 (MQ=255) tATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACt < 4:115565/137‑1 (MQ=255) aaaaaaTAGACCGGTCGTGTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTAGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCCGCCAAGCGTTACCACGCGGACTTCGCt > 6:156848/1‑137 (MQ=255) tAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGGTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCGAAGCGTTACCACGCGGACTTCACTTTGTAg < 6:167684/137‑1 (MQ=255) gtgagtcTACTTTCGGTTTTGGATAATCTATACAACTCTGTTTCGATGGTTGCTAATGATGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTAGCACGCGGACTTCGCTTTGTAGCAGAAg < 8:94283/133‑1 (MQ=255) ggTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACTTCGCTTTGTAGCAGAAg < 8:133271/137‑1 (MQ=255) tcgtcTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACTTCGCTTTGTAGCAGAAGGa < 2:283433/137‑1 (MQ=255) | TCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTACCGTATCAAAAAAATAGACCGGTCGTCTACTTTCGGTTGTGGATAATCTATACAACGCTGTTTCGATGGTTGCTAATGAAGCAGCAGAAGCTTCATTGGCAGCGCATCGACTGATGCAGCCAAGCGTTACCACGCGGACTTCGCTTTGTAGCAGAAGGA > ACB122/2835035‑2835299 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |