| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,638,720 | (A)8→7 | intergenic (‑53/+310) | cspI ← / ← rzpQ | Qin prophage; cold shock protein CspI/Qin prophage; DUF2514 domain‑containing protein RzpQ |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,638,713 | 0 | A | . | 90.9% | 33.1 / ‑1.7 | 11 | intergenic (‑46/+317) | cspI/rzpQ | Qin prophage; cold shock protein CspI/Qin prophage; DUF2514 domain‑containing protein RzpQ |
| Reads supporting (aligned to +/- strand): ref base A (0/1); new base . (6/4); total (6/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCT > NC_000913/1638640‑1638788 | tCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTCGGGCAGAATGGTCCGaaaaaaa < 2:298828/80‑1 (MQ=255) aCTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAat < 1:298829/80‑1 (MQ=255) aGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCagagag > 2:298845/1‑80 (MQ=255) ttttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAAc > 1:298838/1‑81 (MQ=255) ttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAAcc < 1:298830/80‑1 (MQ=255) gACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAgg < 1:298831/81‑1 (MQ=255) gACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAg < 1:298832/80‑1 (MQ=255) agaaATTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat > 2:298839/1‑80 (MQ=255) aaaTTAGCCCTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat < 2:298833/81‑1 (MQ=255) cTTGGGCAGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat > 2:298851/1‑80 (MQ=255) aGAATGGTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGa > 2:298841/1‑81 (MQ=255) ggTCCG‑AAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCt > 2:298850/1‑81 (MQ=255) | TCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCT > NC_000913/1638640‑1638788 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |