| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 684,444 | G→C | 46.1% | T9R (ACA→AGA) | insH1 ← | IS5 transposase and trans‑activator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 684,444 | 0 | G | C | 46.1% | ‑4.4 / 13.2 | 13 | T9R (ACA→AGA) | insH1 | IS5 transposase and trans‑activator |
| Reads supporting (aligned to +/- strand): ref base G (5/2); new base C (4/2); total (9/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01 | |||||||||||
AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAA > CP009273/684360‑684588 | aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTa > 1:384079/1‑149 (MQ=11) taTCTCTTTTCTGGTCTGACGGCGCTTACTGATGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTa > 2:259697/2‑144 (MQ=2) ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGGAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCCGGGACTTGTTCGCACCTTCCTTAGGTAACAtt > 2:315110/1‑144 (MQ=16) aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaa > 1:25263/1‑149 (MQ=255) ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTAATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGAACCTTCCTTAGGTAACATGTAGTGTAGGTACATGTCAAg > 1:42529/1‑149 (MQ=11) cTTACTGCTGAATCCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaaa > 2:189460/1‑144 (MQ=255) tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATGTTg > 1:119677/1‑149 (MQ=255) cccTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTCTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtttt < 2:301765/142‑1 (MQ=255) ttGATGACCCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgcgc < 1:49228/149‑1 (MQ=255) gatgaCTCATGATGAACCCTCTCCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGcgccg < 1:134521/149‑1 (MQ=255) tgaGCCCTGTTCTATGGCTCCAAATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTTCTTAGGTAACCTTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGaa < 1:246398/149‑1 (MQ=255) cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGTGAGATGCGCCGCACCATTCCGAAGCaa > 1:36416/1‑149 (MQ=255) cccTGTTCTATGGCTCCAGATGACAAACAAGATCTCATATCAGGGACTTGTTCGCACCTTCCTGAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGCCGCGCCCCACCATTCCGa > 2:323078/1‑144 (MQ=25) | AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCCGCACCATTCCGAAGCAA > CP009273/684360‑684588 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |