| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,028,373 | A→T | L30Q (CTG→CAG) | ECB_02828 ← | Polysialic acid transport protein kpsM |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,028,373 | 0 | A | T | 97.4% | 150.6 / ‑5.1 | 39 | L30Q (CTG→CAG) | ECB_02828 | Polysialic acid transport protein kpsM |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base T (20/18); total (20/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.87e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CATAATGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCAGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGACTTCAA > REL606/3028300‑3028445 | cATAATGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGa > 2:626762/1‑76 (MQ=255) aTAATGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGAc > 2:51882/1‑76 (MQ=255) aTGTAACCAAAAATGCCCCACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGa > 2:497852/1‑76 (MQ=255) aTGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGa < 1:95314/76‑1 (MQ=255) tGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGaa < 2:690584/76‑1 (MQ=255) gTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAAc < 2:691577/76‑1 (MQ=255) aaCCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACtt < 2:85221/76‑1 (MQ=255) aCCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTa < 1:689888/76‑1 (MQ=255) aaaaTGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCaa < 1:10507/76‑1 (MQ=255) aaaaTGCCCAACAGAATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCaa > 1:384656/1‑76 (MQ=255) aaTGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAg > 2:701226/1‑76 (MQ=255) aTGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGc > 2:695083/1‑76 (MQ=255) gCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCgt < 2:556858/76‑1 (MQ=255) cccAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCgtg < 2:59933/76‑1 (MQ=255) cAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGt > 2:143714/1‑76 (MQ=255) aTCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTAttt < 1:665122/76‑1 (MQ=255) cAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTtct > 2:526181/1‑76 (MQ=255) ncAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCg > 1:694124/2‑76 (MQ=255) aaaTGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTa > 2:221291/1‑76 (MQ=255) aaaTGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCAGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTa < 2:282550/76‑1 (MQ=255) aaTGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAg > 1:348440/1‑76 (MQ=255) tGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGaa < 1:25425/76‑1 (MQ=255) tGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGaa < 2:239356/76‑1 (MQ=255) cgcAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAaat > 1:180652/1‑76 (MQ=255) gcAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAaata > 2:556435/1‑76 (MQ=255) gcAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAaata < 1:121111/76‑1 (MQ=255) cAGAGGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAaataa > 1:271751/1‑76 (MQ=255) agagGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAtataat > 2:706387/1‑76 (MQ=255) agagGGTTCGAGAATCGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAaataat > 1:437779/1‑76 (MQ=255) cGCCCACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACt > 2:142207/1‑76 (MQ=255) cccACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACttt < 2:1966/76‑1 (MQ=255) cccACAAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACttt < 1:665795/76‑1 (MQ=255) cacaAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCt > 1:300849/1‑76 (MQ=255) caAATATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTgg < 2:103914/76‑1 (MQ=255) aaTATCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGAc < 1:593075/76‑1 (MQ=255) atCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTATCTGGACTTc > 2:502135/1‑76 (MQ=255) tCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGACTTCa < 1:661888/76‑1 (MQ=255) tCCCTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGACTTCa < 1:213205/76‑1 (MQ=255) cccTGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGACTTCaa > 1:644717/1‑76 (MQ=255) | CATAATGTAACCAAAAATGCCCAACAGTATCAGCAAATGCGCAGAGGGTTCGAGAATCGCCCACAAATATCCCAGACGGAACTTACCAAAGCGTGTTCGTATTTCTCGTAGAAATAATGCCTCTACGGTGACTTTCTGGACTTCAA > REL606/3028300‑3028445 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |