| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,007,028 | 0 | C | A | 44.4% | 30.6 / 28.8 | 27 | intergenic (‑152/+503) | insC/ECB_01906 | insertion sequence 2 OrfA protein; b1997/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (12/3); new base A (0/12); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.66e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGGGAGCAAGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGCTACG > REL606/2006953‑2007103 | aGGGAGCAAGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATAc < 1:4522973/76‑1 (MQ=255) gggAGCAAGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATACa < 1:3053465/76‑1 (MQ=255) gAGCAAGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATACAAg < 1:6297556/76‑1 (MQ=255) gaaaGCGACTCGGAAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATACAAGcc < 1:2236813/74‑1 (MQ=255) aaGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGGCGGGAAAAACAAGCCaa > 1:6253959/1‑75 (MQ=255) gagaaaaaTTTT‑TCCGGAAAAAGGCG‑GGG‑TTCGGGCAGTACCAGGGGGGGACGAAAAAAAAAAGCCATAGCGTCTg < 2:4994349/73‑1 (MQ=255) acgaaaaaTTTT‑GAGGAAACAACGCC‑GGG‑TAACGGCAATACTAGGGGCGGCCGAAAAATAAAAGCCATAGCGTCTg < 2:655493/73‑1 (MQ=255) gCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGGCGGGAAAAACAAAACATAACGGCTGcc > 1:2554261/1‑76 (MQ=255) aaaaaTTTTAACGAAAAAAA‑GCG‑GGG‑TTACGGCAGTACCGGGGGGGGACGAAAAAAAAAAGCCATAGCGTCTGccc < 2:3636652/76‑1 (MQ=255) aaaaTTTT‑AAGGAAAAAAGGCG‑GGG‑TTACGGCAGTACCGGGGGGGGACGAAAAATAAAAGCCATAGCGTCTGCCCt < 2:3857871/76‑1 (MQ=255) ttttt‑ACGTGAAAAA‑GCG‑GGT‑TTACGGCAGTGCCAGGGGCGGACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGc < 2:3119388/75‑1 (MQ=255) aTTTT‑ACGGGAAAAAGGCG‑GGG‑TTACGGCAGTCCGGGGGGCGGACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAg < 1:4030762/76‑1 (MQ=255) tt‑ACGGAAAAAA‑GCCGGGG‑TTCCGGCGGTACCGGGGGGGGACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAg < 1:1222619/76‑1 (MQ=255) gaaaaCAG‑GCG‑GGG‑TGAGGGCAGTCCCGGGGGGGGAGGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAg < 2:2103082/73‑1 (MQ=255) gTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGGCGGGAAAAAAAAACCACAGAGCCTGTCCCCCGCAgggggg > 2:5429935/1‑73 (MQ=255) ggAAAAA‑GGC‑GGGTTTACGGCAGTACCGGGGGGGGACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAg < 2:3344162/75‑1 (MQ=255) aaccAA‑GCG‑GGG‑TGACGGCGGTACGAGTGGCGGACGAAAAAAACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGAt < 1:5424516/73‑1 (MQ=255) aaaa‑GGC‑GGGTTTACGGCATTACCGGGGGCGGACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATg < 1:4855777/74‑1 (MQ=255) aCAG‑GCG‑GGG‑TTACGGCAGTACCGGGGGCGGAGGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGc < 2:1814140/76‑1 (MQ=255) cAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGGCGGGAAATAAAAGCCATAGCGGCCGCCCCGAGCAGGGGAGATggt > 1:1448610/1‑74 (MQ=255) aa‑GCG‑GGG‑TTCCGGCAGCACCGGGGGAGGGCGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTg < 1:2217128/76‑1 (MQ=255) cc‑GGG‑TTGCGGCAGTACCGGGGGGGGACGAAAAATAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTAt < 1:6163196/76‑1 (MQ=255) c‑GGG‑TTACGGCGTTACCGGGGGGGGACAAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATc < 2:5267003/76‑1 (MQ=255) ggg‑TTCTGGCATTCCCAGGGGGGGAAAAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATcc < 1:4623928/76‑1 (MQ=255) ggg‑TTACGGCAGTACCAGTGGCGGGGGGGAAAAAAAAACCCAAACGGCTGCCCTGAGGCGGGGAGGGGCTGTAttc > 2:2750567/1‑74 (MQ=255) g‑TTACGGCAGTACCGGGGGGGGAGGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCgg < 1:741725/76‑1 (MQ=255) ttgcGGCAGCACAAGGGGGGAACGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGgc < 2:2521823/73‑1 (MQ=255) gCAGTACCAGTGGCGGGCGGGAAAAACAAACCATAGCGTCTTCCCCGAGCCGGGGGGAGGGTGTATTCCGCGcccg > 2:3426910/1‑74 (MQ=255) cAGTACCAGTGGCGGGCGGGAAAAACAAACCAAAGCGGCTGCCCCGCGCCGGGGAGATGGGGTATACGGCGcccgc > 2:2197519/1‑73 (MQ=255) ggTACCGGGGGGGGAAGAAAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCcgacg < 1:2145993/75‑1 (MQ=255) aggggCGGAAGAAAAAAAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACt < 2:3157100/73‑1 (MQ=255) gTGGCGGGCGGGAAAAACAAGCCAAAGCGGCTGCCCTGGGCCGGGGAGACGGTGTATCCGGCGCCGACGACAAcaa > 2:1579137/1‑74 (MQ=255) ggCGGACGGAAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATg > 1:3122031/1‑76 (MQ=255) ggacggaAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAAc > 1:512082/1‑76 (MQ=255) agaaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACaga < 1:230311/73‑1 (MQ=255) agaaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACaga < 1:3310792/73‑1 (MQ=255) ggaAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACagag > 2:4943722/1‑76 (MQ=255) ggaAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACagag > 2:4543232/1‑76 (MQ=255) ggaAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACagag > 2:1249940/1‑76 (MQ=255) gaaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACagag < 1:6301405/74‑1 (MQ=255) aaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGt < 1:3093715/75‑1 (MQ=255) aaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGt < 1:4401360/75‑1 (MQ=255) aaaaaaaaaAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGt < 1:6325151/75‑1 (MQ=255) aaaaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTc > 2:2212576/1‑76 (MQ=255) aaaaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTc > 2:5530649/1‑76 (MQ=255) aaaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCg > 2:2942029/1‑76 (MQ=255) aaaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCg > 1:559124/1‑76 (MQ=255) aaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGgcgc > 2:5529879/1‑76 (MQ=255) aaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGCCGc > 2:1862932/1‑76 (MQ=255) aaTACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGACGc > 1:4853491/1‑76 (MQ=255) aaTAAAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGc < 2:4711466/76‑1 (MQ=255) cAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGCTACg > 2:3636037/1‑76 (MQ=255) cAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGCTACg > 2:2047438/1‑76 (MQ=255) cAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGCTACg > 2:6108462/1‑76 (MQ=255) | AGGGAGCAAGCGACTCGGCAAAAATTTT‑ACGTGAACAA‑GCC‑GGG‑TTACGGCAGTACCAGTGGCGGACGGAAAATACAAGCCATAGCGTCTGCCCTGAGCAGGGTAGATGCTGTATCCGGCGCCGACGACAACTATGAACAGAGTCGCTACG > REL606/2006953‑2007103 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |