Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,008,240 | 0 | T | G | 30.4% | 81.9 / 19.9 | 46 | T41P (ACC→CCC) | ECB_01907 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (7/25); new base G (14/0); total (21/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.85e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGTGTTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAG > REL606/2008166‑2008313 | aGTGTGTTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGgtg < 2:174389/76‑1 (MQ=255) tgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGt < 2:745309/76‑1 (MQ=255) tgtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGt > 2:1445832/1‑76 (MQ=255) gtTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGtt < 2:943125/76‑1 (MQ=255) ttCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTa < 2:411606/76‑1 (MQ=255) tCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAg < 2:770512/76‑1 (MQ=255) tCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGtttg > 2:1454349/1‑74 (MQ=255) cAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGc < 1:1406739/76‑1 (MQ=255) cAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTGGTTTGc > 2:443727/1‑76 (MQ=255) aaCCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGGTAGGAt > 1:1228847/1‑76 (MQ=255) cAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAg < 2:342225/76‑1 (MQ=255) gggAAGTTCAGTGTCGTATTGCGGATGTTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGt < 2:1520973/76‑1 (MQ=255) gggAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGt < 2:1016778/76‑1 (MQ=255) ggAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTg < 2:1624452/76‑1 (MQ=255) ggAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTg < 2:511759/76‑1 (MQ=255) aaGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCg < 2:1023634/76‑1 (MQ=255) aaGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGGTTGCATAAAGGGCg > 2:1146981/1‑76 (MQ=255) aGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTTGGATAAGGGGCgg > 2:1696395/1‑76 (MQ=255) tCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAg < 2:1711883/76‑1 (MQ=255) cGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGt < 1:1584839/76‑1 (MQ=255) cGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGt < 2:1712085/76‑1 (MQ=255) gTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTa < 2:845206/76‑1 (MQ=255) tATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTAc < 1:906731/76‑1 (MQ=255) tATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGGGTTTGGATAAGGTGGGGCCCGTGGAGGAc > 2:903228/1‑76 (MQ=255) aTTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACg < 1:239963/76‑1 (MQ=255) tGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCATAAGGTGCGGCCTGTGCAGTACggg > 2:1720081/1‑76 (MQ=255) tGCGGATGGTATCTGAGTGCACTCCCACACAAGTGGCAGGGGGGGTTAGGATAAGGGGGGGGCCGTGCCGTACggg > 1:1182925/1‑76 (MQ=255) cGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCATAAGGGGGGGCCCGTGGAGGAAggggg > 2:523126/1‑74 (MQ=255) gATGGTATCTGAGGGCACTGCAACACAAGGGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGcc < 1:810576/76‑1 (MQ=255) aTGGTATCTGAGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGGCCTTGCCGTCCAGGCGccc > 2:337512/1‑76 (MQ=255) ggTATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTTGCATAAGGTGCGGCCTTGGCAGGACGGGGGccccc > 2:401293/1‑74 (MQ=255) tATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCa < 2:473290/76‑1 (MQ=255) tATCTGAGTGCACTGCCACACAAGTGGCAGGGGGAGTTAGCATAAGGGGGGGCCTGTGCAGGACGGGCGCCCGCCa > 1:371913/1‑76 (MQ=255) tGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGTAGTACGGGCGCCCGTCAGACg < 1:199130/76‑1 (MQ=255) aGTGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGCCAGACggc > 2:1490940/1‑74 (MQ=255) gTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTcc < 2:1516688/76‑1 (MQ=255) tGCACTGCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGGCAGACGTCCg > 2:1560397/1‑76 (MQ=255) tGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGCTGcc > 2:156321/1‑76 (MQ=255) tGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACCGGCGCCCGTCAGACGGCCGGTggc > 2:353903/1‑74 (MQ=255) gCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGCCAGGCGTCCGTTGccc > 2:355212/1‑75 (MQ=255) gCCACACAAGTGGCAGGGGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGGGCCCGTCAGACGTCCGTTGccc > 2:10623/1‑75 (MQ=255) gCCACACAAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCTGTGCAGTACGGGCGCCCGCCAGACGGCCGTTGCCa > 2:1522226/1‑76 (MQ=255) cgcactAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATc < 2:591615/74‑1 (MQ=255) cacacaAGTGGCAGTTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATc < 1:592485/76‑1 (MQ=255) acacaAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCa < 1:1278509/76‑1 (MQ=255) acaAGTGGCAGGGGTAGTTAGCATAAGGGGCGGCCCGTGCAGTAAGGGCGCCCCCCCGACGTCCGTTGCCACCaaa > 2:1296793/1‑75 (MQ=255) aaGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAagag < 2:228925/76‑1 (MQ=255) gTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCg < 1:1499332/76‑1 (MQ=255) tGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGa < 2:1477875/76‑1 (MQ=255) ggCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGAt > 1:651497/1‑76 (MQ=255) cAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATaa > 2:402466/1‑76 (MQ=255) aGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATaaa > 1:377841/1‑76 (MQ=255) ggTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAg > 1:1313169/1‑76 (MQ=255) | AGTGTGTTCAAACCAGGGAAGTTCAGTGTCGTATTGCGGATGGTATCTGAGTGCACTGCCACACAAGTGGCAGGTGTAGTTAGCATAAGGTGCGGCCTGTGCAGTACGGGCGCCCGTCAGACGTCCGTTGCCATCAAGAGCGATAAAG > REL606/2008166‑2008313 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |