| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,370,044 | (T)9→8 | intergenic (+41/‑172) | pspE → / → ycjM | thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,370,036 | 0 | T | . | 92.0% | 86.3 / 2.4 | 25 | intergenic (+33/‑180) | pspE/ycjM | thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (10/13); total (11/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
AAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCG > NC_000913/1369912‑1370129 | aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac < 2:10043/139‑1 (MQ=255) caaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc < 1:478747/139‑1 (MQ=255) tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc < 1:398832/139‑1 (MQ=255) gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc > 1:207520/1‑128 (MQ=255) gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATccc < 2:207520/128‑1 (MQ=255) ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa > 1:304044/1‑139 (MQ=255) aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTa > 1:173987/1‑113 (MQ=255) aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTa < 2:173987/113‑1 (MQ=255) gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc < 1:294369/112‑1 (MQ=255) gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc > 2:294369/1‑112 (MQ=255) aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac > 1:749669/1‑74 (MQ=255) aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac < 2:749669/74‑1 (MQ=255) cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa > 1:261621/1‑139 (MQ=255) ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATc < 1:314605/79‑1 (MQ=255) ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATc > 2:314605/1‑79 (MQ=255) gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg < 2:805230/139‑1 (MQ=255) gggTCAAAGGTTTGAGTAATTCAGCATGAGCCGCGCTGATGCGCCG‑GTTTTTTTCTGTCACCCCCTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTCCTGCg < 2:534802/138‑1 (MQ=255) aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca < 1:494247/54‑1 (MQ=255) aaGGTTAAAAGAATTCACCAGGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca > 2:494247/1‑54 (MQ=255) aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa > 1:373597/1‑60 (MQ=255) aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa < 2:373597/60‑1 (MQ=255) ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac > 1:80003/1‑43 (MQ=255) ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac < 2:80003/43‑1 (MQ=255) gcTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt < 2:94525/93‑1 (MQ=255) gcTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt > 1:94525/1‑93 (MQ=255) | AAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCG > NC_000913/1369912‑1370129 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |