| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,041,166 | G→T | 13.9% | intergenic (+48/‑209) | yedZ → / → zinT | inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,041,166 | 0 | G | T | 13.9% | 112.7 / 9.1 | 36 | intergenic (+48/‑209) | yedZ/zinT | inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic |
| Reads supporting (aligned to +/- strand): ref base G (15/16); new base T (2/3); total (17/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01 | |||||||||||
TGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA > NC_000913/2041031‑2041301 | tgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCg > 3:243000/1‑149 (MQ=255)tgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCg < 4:23443/149‑1 (MQ=255) tGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATa > 2:228303/1‑149 (MQ=255) ttACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCa < 4:243000/149‑1 (MQ=255) cGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCa < 1:204773/149‑1 (MQ=255) gTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg > 3:197916/1‑149 (MQ=255) taAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa > 3:115917/1‑149 (MQ=255) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGt < 4:329939/149‑1 (MQ=255) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGt > 3:261026/1‑149 (MQ=255) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg > 4:442277/1‑125 (MQ=38) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg < 3:442277/127‑3 (MQ=38) ttCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGt > 2:127409/1‑149 (MQ=255) tttttAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatca > 1:354814/1‑149 (MQ=255) tAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCt < 3:470566/149‑1 (MQ=255) aaCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCaggg > 4:72944/1‑106 (MQ=25) aaCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCaggg < 3:72944/109‑4 (MQ=25) aCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCa > 3:151464/1‑149 (MQ=255) aaGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCaaaa < 4:197916/149‑1 (MQ=255) aGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataat < 2:354814/149‑1 (MQ=255) tCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc > 1:151648/1‑76 (MQ=255) tCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc < 2:151648/76‑1 (MQ=255) acaacaAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttattt < 1:228303/149‑1 (MQ=255) acaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttattttat < 4:115917/149‑1 (MQ=255) acaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttattttat < 4:151464/149‑1 (MQ=255) caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATttattttatt > 2:460455/1‑149 (MQ=255) tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATaa < 4:302124/149‑1 (MQ=255) tGTATAACTTCATGTTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga < 1:460455/149‑1 (MQ=255) atatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAg > 2:57018/1‑149 (MQ=255) atTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGaa > 1:326675/1‑149 (MQ=255) aGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGcc < 1:189987/149‑1 (MQ=255) gcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGct > 2:233542/1‑149 (MQ=255) ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCaaaa < 3:445578/89‑1 (MQ=255) ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCaaaa > 4:445578/1‑89 (MQ=255) ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAACATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTaa < 1:57018/149‑1 (MQ=255) gtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatagat < 4:107554/107‑1 (MQ=255) gtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatagat > 3:107554/1‑107 (MQ=255) | TGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA > NC_000913/2041031‑2041301 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |