Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,622,796 | T→G | L1202R (CTG→CGG) | rhsB → | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,622,796 | 0 | T | G | 97.6% | 105.9 / ‑3.5 | 42 | L1202R (CTG→CGG) | rhsB | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (25/16); total (25/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01 |
TGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGA > NC_000913/3622701‑3622890 | tgatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt < 2:190780/100‑1 (MQ=25) gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt > 1:129541/1‑100 (MQ=18) atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa > 1:709812/1‑100 (MQ=18) tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc < 2:215216/100‑1 (MQ=18) gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa < 2:537474/100‑1 (MQ=18) gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa > 2:1095124/1‑100 (MQ=17) ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg > 1:554967/1‑100 (MQ=18) ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg > 1:244165/1‑100 (MQ=18) ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc > 1:790779/1‑100 (MQ=17) ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc > 1:748701/1‑100 (MQ=18) cAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc > 1:3202/1‑100 (MQ=18) ccTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGc < 1:734407/100‑1 (MQ=18) gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat > 1:657762/1‑100 (MQ=18) gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGGGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat < 2:701004/100‑1 (MQ=255) gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat < 1:72020/100‑1 (MQ=18) atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa > 2:728607/1‑100 (MQ=18) aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc > 1:178682/1‑100 (MQ=18) ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg > 1:1022203/1‑100 (MQ=18) cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa < 2:940885/97‑1 (MQ=17) cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc > 1:265168/1‑100 (MQ=18) aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg > 1:73135/1‑100 (MQ=18) tCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa < 1:583578/100‑1 (MQ=17) aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata < 1:285740/100‑1 (MQ=18) gCTGCAGCAGCTTATACGGCTGCTGGGGCAGCAGTATGATGAGGAGTCCGGACGGAATTACAACCGCCACCGCTCTTATGACCCGCGGCAGGAGCGata > 2:317525/1‑99 (MQ=11) gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc > 1:1017011/1‑100 (MQ=18) agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa < 2:270078/99‑1 (MQ=17) agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc < 2:265165/100‑1 (MQ=18) tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc > 2:43423/1‑100 (MQ=17) ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt > 1:1032740/1‑100 (MQ=14) ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt < 1:461007/100‑1 (MQ=14) ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc < 2:1096649/100‑1 (MQ=14) agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg > 1:593125/1‑100 (MQ=25) atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg > 1:135970/1‑100 (MQ=255) atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg > 1:758828/1‑100 (MQ=255) tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg < 2:193865/100‑1 (MQ=255) tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg > 1:230699/1‑100 (MQ=255) gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGGTGGAAcc < 1:1201256/100‑1 (MQ=255) gagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGt > 2:103129/1‑100 (MQ=255) agTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa < 2:677177/100‑1 (MQ=255) gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt > 1:438490/1‑100 (MQ=255) cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa > 2:294974/1‑100 (MQ=255) cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa < 2:709807/100‑1 (MQ=255) | TGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGA > NC_000913/3622701‑3622890 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |