Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 4,413,553 | 0 | A | C | 34.0% | 59.4 / 10.5 | 47 | I7L (ATC→CTC) | CJ019_04183 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (21/10); new base C (0/16); total (21/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.49e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAACCCCCCTAAGGCTA > ACB122/4413424‑4413688 | cccTGACGGGCCGCGCTACAAGGCGAGCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCGGGAGCGCCCGGCGCATCCAACAACAACTTTAACGTCTTGCTTGAGGTTACCCATGCCCACATAAAACAAAATCGACt < 3:309824/136‑1 (MQ=255) cgcTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGca > 5:304899/1‑136 (MQ=255) gcTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCCCAGAAAACAAAATCGACTGCCCGGCATTGc > 2:315192/1‑134 (MQ=255) cTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGTTCGGCCGGCGCATCCAACAACAACTTGACCGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGcaca < 7:243677/136‑1 (MQ=255) aaGGCGATCGGAAACAGCAAGGCCGTGTTCGGCGTCCGCTGGATCGGCCGGCGCAGCCAACAACAACTTGACCGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAACTCGACTGCCCGGCATTGCACAAGCg < 8:253321/136‑1 (MQ=255) tgggaaaCAGCACGGCCGGGTGGGTCGGCCGCGGGATCGGCGGGCGCATCCAACAACAACTTCAAGGTCTTGCTTGTGGTTACCCAGGCCCACATAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTc < 3:327468/131‑1 (MQ=255) gCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGc > 4:209051/1‑136 (MQ=255) aGCATGGCCGTGTTGGTCGTCCGCGGTTTCGGCGGGCGCATCCAACAACAACTTGAAGTTCTTGCTTGAGGTAACCAATGCCCACATAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTAccc < 6:121574/136‑1 (MQ=255) ccGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGc > 5:287255/1‑136 (MQ=255) ttCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcg > 5:169958/1‑136 (MQ=255) ttCGTCGTCCGCTGGAGCGGCCGGCGCATCCAACAAAAACTTGAACTTTTTGCTTGAGGATACCAATCCCCACACAAAAAAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcg < 6:177214/136‑1 (MQ=255) tcgtcgTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgc > 6:159356/1‑136 (MQ=255) tcgtcgTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgc > 5:7165/1‑136 (MQ=255) tcgtcgTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgc > 6:318468/1‑136 (MQ=255) tcgtcgTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcg > 2:59831/1‑135 (MQ=255) ggggggccGCTGGAGGGGCGGGCGCAAACAACAACAAGTTCACCGGATTGGTTGCGGTAACCCATCCCCCCACAAAAAAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgcg < 5:328840/130‑1 (MQ=255) cgccgcccGCGGGATCGGCCGGCGCATCCAACACAAACTTCAACGGCTTGCTTGCGGTAACCCATGCCCACACAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgc < 8:44819/129‑1 (MQ=255) gtcggcCGGTGGAGCGGCCGGCGCAGCCAACAACAACTTGAACGGCTTGCTTGCGGTAACCCATGCCCAAATAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcgcg < 2:228147/135‑1 (MQ=255) cGCTGGATCGGCCGGCGCATCCAACAACAAGTTGAACGTCTTGCTTGCGGTTCCCCCTGCCCACCCAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCcgcg < 6:142479/136‑1 (MQ=255) gcggccgggGCATCCAACAACAACTTGACCGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTc < 6:150767/135‑1 (MQ=255) gccggcGCTTCAAACAACAACTTGACGGTCTTGCTTCTGGTTACCCATGCCCACAGAACACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGc < 4:26093/136‑1 (MQ=255) caacTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCa > 2:286624/1‑135 (MQ=255) aacTTGAACGTCTTGCTTGAGGTTACCCTGCCCCACACAAAAAAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCAcc < 5:134628/136‑1 (MQ=255) aaattcaaCGTCTTGCTTGCGGTACCCCAGGCCCACACAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCAcc < 5:175697/130‑1 (MQ=255) gAACGTCTTGCTTGAGGTTCCCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCAccgcc < 3:71991/134‑1 (MQ=255) gAACGTCTTGCTTGAGGTTACACATGCACACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCAccgccg > 7:272601/1‑135 (MQ=255) aCGTCTAGCTTGAGGTTACCCATGCCCACACAAAACAAAATCGACTGCCCGGCATTGCACAACCGAATCTTCGGCTACCCGTACGGCGATCGCGTACCGCGTACCGTCCGCATGCTGGTTACCGTCACCGCCGAAc > 7:45823/1‑136 (MQ=255) tgcttgATTTTCCCCTTCCACACATAAAACAAAATCGACTTCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCTAGTTACGCTTACTTTCATCATGCTTTTTAACGTCACCGCCGAACGCAt < 2:276057/134‑1 (MQ=255) tgcttgAGGTTACCCATGCCCACAGAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCa < 3:294499/106‑1 (MQ=255) tgcttgAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCa > 4:294499/1‑106 (MQ=255) taaCCCATCCCCACATAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCg < 6:169958/134‑1 (MQ=255) aCCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTc > 8:270097/1‑133 (MQ=255) cATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGAGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATg > 8:58401/1‑136 (MQ=255) aTGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATgg > 4:129701/1‑136 (MQ=255) cccACAGAAAAAAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCa < 3:89814/136‑1 (MQ=255) cccACAAAAAAAAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGcc < 6:7165/135‑1 (MQ=255) cacaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcg > 1:226433/1‑66 (MQ=255) cacaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATcg < 2:226433/66‑1 (MQ=255) cacaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCacc > 7:76227/1‑134 (MQ=255) cacaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAc > 5:267627/1‑136 (MQ=255) cacaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAc > 4:43001/1‑136 (MQ=255) acaGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcc > 7:243649/1‑136 (MQ=255) aGAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcccc < 8:277066/136‑1 (MQ=255) aGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTTACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcccc < 6:104475/136‑1 (MQ=255) gAAAACAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcccc < 2:87759/135‑1 (MQ=255) caaaaaaaaATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcccc < 2:206197/134‑1 (MQ=255) caaaaCAAACTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAAcccc < 1:58768/134‑1 (MQ=255) aTCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAACCCCCctaaggcta < 4:178854/136‑1 (MQ=255) | CCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTTGAACGTCTTGCTTGAGGTTACCCATGCCCACAGAAAACAAAATCGACTGCCCGGCATTGCACAAGCGAATGTTCGGCTACCCGTACGGCGATCGCGTTCCGCGTACGGTCCGCATGCTGGTTACCGTCACCGCCGAACGCATCCCTGGTGTCGAGATGGCCAACCCCCCTAAGGCTA > ACB122/4413424‑4413688 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |