Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 1,721,869 | 0 | G | C | 26.0% | 87.0 / 11.8 | 50 | V98V (GTC→GTG) | CJ019_01617 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (30/7); new base C (0/13); total (30/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.18e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCT > ACB122/1721733‑1721999 | gCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCg < 5:156023/137‑1 (MQ=255) tATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCa < 2:222039/137‑1 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGcc > 5:120458/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGcc > 5:149211/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGcc > 5:170455/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGcc > 5:363104/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCACCGTCTTACTCGACCCAGcc > 8:139715/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGGCAGCAT‑CTTATTCCCCGTCTGACTCGACCCCGcc > 7:127463/1‑137 (MQ=255) tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGAAGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGcc > 5:305284/1‑137 (MQ=255) ggTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCAcccc > 7:94314/1‑136 (MQ=255) gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATTATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCAGCAT‑CTTCTTACCCCTCTTCCTCGACCCCGCCAccccg > 3:16240/1‑135 (MQ=255) gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCgg > 3:364127/1‑137 (MQ=255) aaaaCCAATTCATACAGAAACTAGAAAGCCTGCGTTTTAATGTTGTACGCCAGGTGTTTTCCTGTTCACGGTCGGCGCTGATGCGGACCCTGCGGTCAGCTT‑CTTCTTCCCCGTCTTCCTCGACCCCGCCACCCggg < 4:62787/136‑2 (MQ=255) aCCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGgcgcg > 5:326921/1‑137 (MQ=255) ccTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGCCCCCGCCCCCCGGcgcgccgcgcc > 2:380310/1‑136 (MQ=255) tACAGAAACGAGCAAGCCAGCCGTTTCTTGATGTAAGCCAGGTGAATTCCTCTTCACCGTCGGCGCAGATTCGGACCCAGCCGTGAGGAA‑CTTCTTCCCCCTCTTACTCGACCCATCCACCCGGCCCGCCGCCCAcg > 5:254092/1‑136 (MQ=255) aaaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCcac > 1:52263/1‑137 (MQ=255) aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCcccaccca > 8:40435/1‑137 (MQ=255) aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCATCGACACATCCACCCGCCGCGCCGCGCACTTAAACAccaa > 1:43162/1‑135 (MQ=255) aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAA‑CTTCTTACCCGTCTTCCTCGACCCCGCCACCCCGCGCGCAGCGCACTTACACAccag > 1:330009/1‑135 (MQ=255) aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGCCCCAGCCACCCGGCGCGCAGCGCCCTTCCACACCCAgcg > 2:150091/1‑137 (MQ=255) aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAgcc > 5:286343/1‑136 (MQ=255) tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTACCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGgcgc > 2:89294/1‑137 (MQ=255) gCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAg > 1:44487/1‑137 (MQ=255) cccggTCATTCCCTCTCCCCGGGCGGGGATGATGCGGACCCAGGGGTCAGCAT‑CTTCTGCCCCGTCTACCTCGCCCCACCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAgg < 1:250437/134‑1 (MQ=255) taTTCCCTCTTCACGGGGGGCGGTGTGGCGCCCCCAGGGGTCA‑CATGGTTCTCCCCCGCCTCCCTCCCCCCAGCCCCCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGt < 4:88838/136‑1 (MQ=255) ccTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCg > 5:394115/1‑137 (MQ=255) cACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCAc > 6:372819/1‑137 (MQ=255) cGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCCGCGCACTTCCACCCCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACggg > 8:59284/1‑136 (MQ=255) cGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCCCGGCAAATGCGGCGCGGGCGCCGGCAACCTCATCGGCCCACCGc > 1:171612/1‑137 (MQ=255) gcTGATGCGGACCCCGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCCCCCCCCCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt < 5:386319/137‑1 (MQ=255) ttggggACCCAGCGGTCTGCTT‑CTTATTCCCCGTCTTCCTCCACCCAGCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt < 6:326921/133‑1 (MQ=255) tggcgcACCCCGCGGTCAGATT‑GTTCTACCCGGTCTTCCTCCACCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt < 4:189800/135‑1 (MQ=255) tgccGGCCCCCGCGGTCAGCAT‑GTTCTCCCCCGTCTCCCCCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt < 3:359209/134‑1 (MQ=255) ctGCGGACCCCGCGGTCATCAT‑CTTCTGCCCCGTCTTCCTCGCCCCCGCCCCCCGGCGCGCCCCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt < 6:316923/136‑1 (MQ=255) aggCGGCCCCCGGGGCCACCAT‑TTTCTCCCCCGGCTCCCCCCCCCCCCCCCCCCGGCGCGCACCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt < 3:353543/135‑1 (MQ=255) gCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCttt > 7:69017/1‑137 (MQ=255) ggTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTCTGGTCGCTGTCCAGATACCgg > 6:337482/1‑137 (MQ=255) ggTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCATATACAgg > 3:133294/1‑137 (MQ=255) ggTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAgg > 3:139250/1‑137 (MQ=255) ggTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGAGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCGTTCCAGATACAgg > 6:332736/1‑137 (MQ=255) gaaT‑CTTCTTCCCCGTCTTCCTCCACCCACCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGtt < 8:94314/135‑1 (MQ=255) att‑cttctTCCCGTTCTCCCTCGCCCCCCCCCCCGGGGGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTc < 5:302884/135‑1 (MQ=255) aT‑CTGCTCCCCCGTCTACCTCCCCCCCCCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTcc < 8:42200/137‑1 (MQ=255) t‑cttctCCCCCGTCTCCCCCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCa < 4:259714/137‑1 (MQ=255) cttcCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTc < 6:363104/137‑1 (MQ=255) ttcCCCTTCTTCCTCCACCCCCCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc < 4:170503/137‑1 (MQ=255) ttcCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc < 2:44487/137‑1 (MQ=255) taccccGCCTCCCACCACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc < 6:120458/135‑1 (MQ=255) ggTCTTCCTCCCCCCACCCCCCGGGCGGCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACa < 4:16240/136‑1 (MQ=255) ttCCTCCACCCCCCCACCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCt < 1:286626/137‑1 (MQ=255) | GCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAT‑CTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCT > ACB122/1721733‑1721999 |
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: ‑ |