Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,536,106 | T→C | E261G (GAA→GGA) | cysP ← | thiosulfate/sulfate ABC transporter substrate‑binding protein CysP |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,536,106 | 0 | T | C | 96.4% | 85.1 / ‑4.1 | 28 | E261G (GAA→GGA) | cysP | thiosulfate/sulfate ABC transporter substrate‑binding protein CysP |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (20/7); total (20/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01 |
GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA > NZ_CP009273/2536021‑2536190 | gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt > 2:877766/1‑90 (MQ=255) tATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAc < 1:665081/90‑1 (MQ=255) aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGAAcc > 1:872742/1‑90 (MQ=255) cgcgGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTgg > 2:289515/1‑90 (MQ=255) gtaATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCAc < 2:373214/90‑1 (MQ=255) tAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGttt < 2:832845/90‑1 (MQ=255) aGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtttt > 1:271488/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:825259/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:100215/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:996079/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:914081/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:412107/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 2:415129/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:577452/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:420050/1‑90 (MQ=255) tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa > 1:328112/1‑90 (MQ=255) ggTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGc < 1:547790/90‑1 (MQ=255) tGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAt > 2:291061/1‑90 (MQ=255) ataGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAg < 1:877766/90‑1 (MQ=255) cAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTg < 1:291061/90‑1 (MQ=255) aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 2:924025/1‑90 (MQ=255) aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 1:91394/1‑90 (MQ=255) aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg > 1:846286/1‑90 (MQ=255) gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc > 2:755517/1‑90 (MQ=255) gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACCTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCTGAATc > 1:588890/1‑90 (MQ=255) gCTTTGACGGCTTTTCCCGTACCGCTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc > 2:429090/1‑90 (MQ=255) ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc < 1:635920/90‑1 (MQ=255) cTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCaa < 1:580090/90‑1 (MQ=255) | GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA > NZ_CP009273/2536021‑2536190 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |