| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,426,592 | T→A | 35.3% | intergenic (+35/+36) | rplI → / ← yjfZ | 50S ribosomal subunit protein L9/uncharacterized protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,426,592 | 0 | T | A | 35.3% | 13.4 / 15.7 | 17 | intergenic (+35/+36) | rplI/yjfZ | 50S ribosomal subunit protein L9/uncharacterized protein |
| Reads supporting (aligned to +/- strand): ref base T (6/5); new base A (3/3); total (9/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01 | |||||||||||
CGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCT > NC_000913/4426453‑4426740 | cGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCg < 3:1152519/149‑1 (MQ=255) gTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGc < 3:1329246/126‑1 (MQ=255) gTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGc > 4:1329246/1‑126 (MQ=255) cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGc > 3:1860418/1‑149 (MQ=255) aaGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAtt < 4:70185/101‑1 (MQ=255) aaGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAtt > 3:70185/1‑101 (MQ=255) gAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGAt < 4:1860418/149‑1 (MQ=255) ttCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt < 1:827230/103‑1 (MQ=255) ttCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt > 2:827230/1‑103 (MQ=255) aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCt < 4:466576/110‑1 (MQ=255) aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCt > 3:466576/1‑110 (MQ=255) aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa < 4:925915/99‑1 (MQ=255) aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa > 3:925915/1‑99 (MQ=255) gACCATTGGTTGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAc > 1:941492/1‑149 (MQ=255) aCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAACATTCTCtt > 3:1168220/1‑75 (MQ=255) aCCATTGGTCGGCGTTCTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt < 4:1168220/75‑1 (MQ=255) cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGata > 4:492412/1‑110 (MQ=255) cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcg > 3:1420932/1‑149 (MQ=255) cATTGGTCGACGTTTTGCTTTCTATTTTTCGTCTGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATACGTCAGGATGGTTTGCTGGCGTTAGata < 3:492412/110‑1 (MQ=255) ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgct > 2:1165832/1‑149 (MQ=255) | CGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCT > NC_000913/4426453‑4426740 |
| 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: ‑ |