| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,029,134 | T→C | intergenic (‑673/+388) | ECB_02828 ← / ← yghD | Polysialic acid transport protein kpsM/predicted secretion pathway M‑type protein, membrane anchored |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,029,134 | 0 | T | C | 89.9% | 241.5 / 13.8 | 69 | intergenic (‑673/+388) | ECB_02828/yghD | Polysialic acid transport protein kpsM/predicted secretion pathway M‑type protein, membrane anchored |
| Reads supporting (aligned to +/- strand): ref base T (3/4); new base C (36/26); total (39/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.90e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01 | |||||||||||
AACCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTCCAA > REL606/3029062‑3029205 | aaCCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTc > 2:1562854/1‑76 (MQ=255)aaCCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTc < 1:964630/76‑1 (MQ=255)aaCCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTc < 1:417816/76‑1 (MQ=255) aCCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCa < 1:108292/76‑1 (MQ=255) aCCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCa < 1:64477/76‑1 (MQ=255) ccAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCAc > 1:693842/1‑76 (MQ=255) cAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACt > 2:999400/1‑76 (MQ=255) aGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTg < 1:618161/76‑1 (MQ=255) aGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAACAGTTTCACTg < 1:706473/76‑1 (MQ=255) gtTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGaa > 1:543580/1‑76 (MQ=255) ttACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGaaa > 2:49946/1‑76 (MQ=255) ttACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGaaa > 1:434224/1‑76 (MQ=255) tACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAat > 2:975351/1‑76 (MQ=255) tACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAat < 1:1260230/76‑1 (MQ=255) caaccaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAatta > 2:1885245/1‑76 (MQ=255) caaccaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAatta > 2:42480/1‑76 (MQ=255) aaccaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAattat > 2:1664465/1‑76 (MQ=255) aaccaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAattat > 2:199607/1‑76 (MQ=255) accaacGCCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTGGAAattatt < 2:1276845/76‑1 (MQ=255) accaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTGGAAattatt > 1:599145/1‑76 (MQ=255) accaacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAattatt < 2:324676/76‑1 (MQ=255) naacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATtattat > 1:1576907/2‑76 (MQ=255) caacACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATtattat > 2:439549/1‑76 (MQ=255) cacCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGca < 2:239683/76‑2 (MQ=255) tCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTangt < 2:1326220/76‑4 (MQ=255) tCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTAttt < 1:1699162/76‑1 (MQ=255) cccTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTGGAAATTATTATGCTTATTTg > 2:1702320/1‑76 (MQ=255) cccTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTg > 2:991397/1‑76 (MQ=255) aGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAAtt < 2:1422501/76‑1 (MQ=255) tAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTcc > 2:605460/1‑76 (MQ=255) tAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTcc > 2:1481122/1‑76 (MQ=255) tAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTcc > 2:1551384/1‑76 (MQ=255) atatACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCtt > 2:1307299/1‑76 (MQ=255) atatACCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCtt < 1:889373/76‑1 (MQ=255) ataCCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTa < 1:232810/60‑1 (MQ=255) ataCCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTa > 2:232810/1‑60 (MQ=255) ataCCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTc < 2:153755/76‑1 (MQ=255) aCCTGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTg < 1:592219/76‑1 (MQ=255) ccTGAGCAGACATGTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGa < 2:733654/76‑1 (MQ=255) tGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGAtt > 1:1197006/1‑76 (MQ=255) tGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGAtt < 2:1524274/76‑1 (MQ=255) tGAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGAtt > 1:928430/1‑76 (MQ=255) gAGCAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTa < 2:396382/76‑1 (MQ=255) cAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTAtct > 2:1552616/1‑76 (MQ=255) cAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTAtct > 1:711998/1‑76 (MQ=255) cAGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTAtct > 2:226769/1‑76 (MQ=255) aGACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTAtctc < 2:1845265/76‑1 (MQ=255) gACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTcc < 1:1232946/76‑1 (MQ=255) gACATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTcc < 2:377086/76‑1 (MQ=255) gACACTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTCTCTGATTATCTcc < 2:377097/76‑1 (MQ=255) aCATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCa < 2:741404/76‑1 (MQ=255) cATTTCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAg < 2:998410/76‑1 (MQ=255) tttCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGa > 2:816767/1‑76 (MQ=255) ttCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGaa < 1:69601/76‑1 (MQ=255) tCTTAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAt > 2:697501/1‑76 (MQ=255) tAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAtttt > 2:567756/1‑76 (MQ=255) tAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAtttt > 1:797168/1‑76 (MQ=255) tAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAtttt > 1:49684/1‑76 (MQ=255) tAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAtttt > 2:106660/1‑76 (MQ=255) tAAAATCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAAtttt > 1:1000240/1‑76 (MQ=255) aaaTCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTAc > 2:679997/1‑76 (MQ=255) aaTCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACt < 2:532160/76‑1 (MQ=255) tCAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTaa > 2:911947/1‑76 (MQ=255) cAACGAATAGTTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAAt > 2:910047/1‑76 (MQ=255) cAACGAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAAt < 1:1562854/76‑1 (MQ=255) gAATAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTcc > 2:668011/1‑76 (MQ=255) aaTAGTTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTCCa < 2:1592358/76‑1 (MQ=255) aaTAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTCCa > 1:1262527/1‑76 (MQ=255) aTAGCTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTCCaa > 2:860181/1‑76 (MQ=255) | AACCAGTGTTACAACCAACACCATCCCTCAGTAGATATACCTGAGCAGACATTTCTTAAAATCAACGAATAGTTTCACTGGAAATTATTATGCTTATTTGAAATTCCCTTTCTGATTATCTCCAGGAATTTTACTAATTTCCAA > REL606/3029062‑3029205 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |