Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,954,706 | T→G | K107T (AAA→ACA) | motB ← | flagellar motor protein MotB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,954,706 | 0 | T | G | 97.1% | 136.1 / ‑3.7 | 35 | K107T (AAA→ACA) | motB | flagellar motor protein MotB |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (19/15); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.57e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01 |
TTTCGGATCGGACTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTTTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCACCGCCGG > REL606/1954633‑1954775 | tttCGGATCGGTCTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTc < 2:394978/76‑1 (MQ=255) tttCGGATCGGACTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTc > 2:69267/1‑76 (MQ=255) tttCGGATCGGACTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTc > 2:622981/1‑76 (MQ=255) tctcggatcggaCTCTATAAACTGGTTGACCTCACCCCACTATTNCCGCAATCGANTTTGCTCCATGCGTTTTGTc < 1:122035/74‑1 (MQ=255) ttCGGATCGGACTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCa < 2:514742/76‑1 (MQ=255) cggatcggaCTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAga > 2:122035/1‑75 (MQ=255) atcggaCTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCNGctct > 2:963683/1‑76 (MQ=255) tcggaCTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCtt > 2:621691/1‑76 (MQ=255) aCTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGAt < 2:830379/76‑1 (MQ=255) tCAACTGGTCGAGATCACCCCGCANTTTCCGAAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCgg > 2:199539/1‑76 (MQ=255) tCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCgg > 1:1154321/1‑76 (MQ=255) cTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGgctgc > 2:882880/1‑76 (MQ=255) tGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGgctgct > 2:273719/1‑76 (MQ=255) tCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTAt > 1:437185/1‑76 (MQ=255) cGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTAtt > 1:69297/1‑76 (MQ=255) cGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTAtt > 2:279211/1‑76 (MQ=255) gagaTCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTc < 2:569146/76‑1 (MQ=255) tCACCCAGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTTTTCACtt > 2:823445/1‑76 (MQ=255) cACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTTACTTc > 2:409762/1‑76 (MQ=255) cGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCt < 1:788551/76‑1 (MQ=255) cGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCt < 1:530577/76‑1 (MQ=255) ncAAGTTCCGCAATCGACTTTNCTCCATNCGTTTGGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTg < 1:1004467/75‑1 (MQ=255) ttCCGCAATCNACTCTGCTCCATGCGTTTTGTCNGCTCTTCGNTGTTCGGCTGNTTATTCACTTCCCCCAGGCtgt > 2:1041355/1‑74 (MQ=255) tCCGCAATCGACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTg < 2:101915/76‑1 (MQ=255) cAATCGACTTTGATCCATGCGTTTTGTCAGCTCTTCGATGTTCGACTGCTTATTCACTTCCCCCTGGCCTTGGgtg > 2:180434/1‑76 (MQ=255) gACTTTGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAatc < 1:511503/76‑1 (MQ=255) ttGCTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAAtcatca < 1:23160/76‑1 (MQ=255) cTCCATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATTATcacca < 1:320408/76‑1 (MQ=255) tCCATGCGTTTTTTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATcaccac < 1:814527/76‑1 (MQ=255) ccATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCaccacc > 1:819868/1‑76 (MQ=255) cATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCAccg < 2:501431/76‑1 (MQ=255) cATGCGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCAccg < 1:199539/76‑1 (MQ=255) cGTTTTGTCAGCTCTTCGATGTTCTGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATAACCACCGCCgg > 2:831076/1‑76 (MQ=255) cGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCACCGCCgg < 2:438576/76‑1 (MQ=255) cGTTTTGTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCACCGCCgg > 2:152979/1‑76 (MQ=255) | TTTCGGATCGGACTCTATCAACTGGTCGAGATCACCCCGCAATTTCCGCAATCGACTTTGCTCCATGCGTTTTTTCAGCTCTTCGATGTTCGGCTGCTTATTCACTTCCCCCTGGCTTTGGGTGTAATCATCACCACCGCCGG > REL606/1954633‑1954775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |