| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,983,700 | 0 | T | G | 34.2% | 71.8 / 13.7 | 41 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (5/22); new base G (14/0); total (19/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.30e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAT > REL606/3983629‑3983772 | tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGcc < 1:651501/75‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATTGGCA‑GGGGTGCCt < 1:1442046/75‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGccc > 1:1139531/1‑74 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTc > 1:1498230/1‑75 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTccc < 2:1470884/75‑1 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTccc > 2:1212814/1‑75 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 1:284917/75‑1 (MQ=255) gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcaca < 1:727575/75‑1 (MQ=255) gcgTTACCGCCAAGGGTGAGTGCTTGCATCCTACGGCTGGTATTGGTGGCATCCGCA‑GGGGTGCCTCCCACAAAt < 2:124498/73‑1 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCaacct > 2:1504495/1‑72 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGGGGGATCCGCA‑GGGGTGCCGCCCACAAAgtt > 2:1443799/1‑75 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGAc < 2:656846/75‑1 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAATTTGAc > 1:1070641/1‑75 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAattttc > 2:926753/1‑72 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTTAc > 1:579092/1‑75 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAATTTAGGATc > 2:913160/1‑75 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCCAAATTT‑ACGAACGc > 2:1377984/1‑75 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCACAAAATTGACGATCGc > 1:674247/1‑75 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTTACGATCGc > 2:68931/1‑75 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTg < 2:1597354/75‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAATTTACGATCGCTgggg > 1:874590/1‑73 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:1274321/75‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAAATTTACGAACGCTGGGGATGa > 1:424064/1‑75 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGc < 1:930162/75‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGcc < 1:363123/75‑1 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAAATTGACGATCGCTGGGGATGACGcc > 1:1354197/1‑75 (MQ=255) cagccTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCAt < 2:412297/72‑1 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCAt < 1:1349754/75‑1 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 1:402363/75‑1 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCACAAAATTGAAAAACGCTCGGGATGAAGCCATcg > 1:372578/1‑75 (MQ=255) tgcGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa < 1:1030349/73‑1 (MQ=255) gcGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGaa < 2:1377952/74‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 2:949403/75‑1 (MQ=255) cTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca > 1:1157661/1‑75 (MQ=255) tGTTATTGGTGGCATCGGCA‑GGGGTGCCCCCCACAAAATTGACGATCCCTGGGGATGACGCCATCCCGAAAcccc > 1:1140351/1‑73 (MQ=255) gTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacc > 1:837433/1‑75 (MQ=255) ttATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc > 2:444899/1‑75 (MQ=255) ttATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCCGGGGATGACGCCATCGCGAAGccaccc > 2:317460/1‑75 (MQ=255) tATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccc > 2:962964/1‑74 (MQ=255) tggtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 1:296221/75‑1 (MQ=255) ggtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc < 1:444899/75‑1 (MQ=255) gtggCATCGGCA‑GGGGTGCCTCCCCCAAATTTGAAGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt > 1:160491/1‑75 (MQ=255) tggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg > 2:1483401/1‑75 (MQ=255) tggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg > 2:1095481/1‑75 (MQ=255) aTCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGt < 1:1232693/75‑1 (MQ=255) tggCA‑CGTGTGCATACAACAAATTTGACCATCGCTGGAGAGGACGCCAGCGCGAAGCCACCCACACCTGCGGttt < 2:181839/74‑1 (MQ=255) cGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 2:181850/75‑1 (MQ=255) ggCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 2:970971/75‑1 (MQ=255) gCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa > 1:571498/1‑75 (MQ=255) cA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg < 1:1035355/75‑1 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAt > 1:1254302/1‑75 (MQ=255) | TGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAT > REL606/3983629‑3983772 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |