| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | REL606 | 3,684,366 | C→A | 100% | R216S (CGT→AGT) | ECB_03438 → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,684,366 | 0 | C | A | 97.4% | 122.9 / ‑2.8 | 38 | R216S (CGT→AGT) | ECB_03438 | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (20/17); total (21/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CATATCGACCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACATTTGTCCCT > REL606/3684222‑3684490 | cATATCGACCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggt > 2:15894/1‑150 (MQ=255) tCGACCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgcgc < 1:432016/150‑1 (MQ=255) tCGACCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgcgc < 2:584268/150‑1 (MQ=255) aCCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAAc < 2:225314/150‑1 (MQ=255) ttgTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGGATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTa < 1:449662/148‑1 (MQ=255) gTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg > 2:653786/1‑150 (MQ=255) gTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg > 2:924/1‑150 (MQ=255) tataagagacagaagtTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAc < 1:641486/138‑1 (MQ=255) cAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAcc < 1:80530/150‑1 (MQ=255) aGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCa > 2:607652/1‑150 (MQ=255) gATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCaa > 2:262550/1‑150 (MQ=255) aagtTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAcctgtctcttata > 2:643370/1‑139 (MQ=255) aagtTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaa > 2:20049/1‑150 (MQ=255) gtTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaaaa > 2:395679/1‑150 (MQ=255) acacGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAgcg < 1:312589/150‑1 (MQ=255) cGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCgg > 2:400231/1‑150 (MQ=255) cGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCgg > 2:184535/1‑150 (MQ=255) cGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCgg > 2:502431/1‑150 (MQ=255) cGGTTACCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCgg > 1:553458/1‑150 (MQ=255) aTCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGa > 1:597017/1‑150 (MQ=255) cTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACAt < 1:606545/150‑1 (MQ=255) aaTTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCacac < 2:598775/150‑1 (MQ=255) ttGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCt < 1:647408/150‑1 (MQ=255) tGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGCGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGcc < 2:555062/150‑1 (MQ=255) gcgcCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGAt < 1:261804/150‑1 (MQ=255) gcgcCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGAAGGGAGTGGAt < 1:651865/150‑1 (MQ=255) tGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGa > 2:565686/1‑150 (MQ=255) tGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACGCTTACGACCAAGAATATGACAAGCGCGGGCAGACCTCGGACTGGCACACCTACGGCCGGGCGTGGATGGTGa > 1:30865/1‑150 (MQ=255) gTATGAGGTGACTGAAGAGCCGCGCTACCTGGCGCTGCCGAACTACTTTGTCGAACAGAGTGGTGCGCACCCGCACACCTACGACCAAGAACAAGAAAAGCGCGGGCAGCAGACGCCCTGGCACACACACGGCCCGCCGGGGTAGGTGaa > 1:152562/1‑150 (MQ=255) gTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGaa > 2:206755/1‑150 (MQ=255) gAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTAc > 2:365278/1‑150 (MQ=255) taagagacagaCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGc < 1:195235/140‑1 (MQ=255) gagCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGc < 2:230678/150‑1 (MQ=255) aCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGcctgtctctta > 2:195798/1‑141 (MQ=255) aCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACAtt > 2:312583/1‑150 (MQ=255) aCCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAACAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACAtt > 2:94710/1‑150 (MQ=255) gcgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGAACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACATTTGTcc < 1:15843/150‑1 (MQ=255) gcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACATTTGTCCCt < 1:311701/150‑1 (MQ=255) | CATATCGACCGCGTATTTGGTCCAGATGAAAGTAAGTTACACGGTTATCCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAAAGCCTACAGCCAGGCACATTTGTCCCT > REL606/3684222‑3684490 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |