Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 214,743 | A→C | Y174D (TAC→GAC) | proS ← | proline‑‑tRNA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 214,743 | 0 | A | C | 95.0% | 70.4 / ‑3.1 | 20 | Y174D (TAC→GAC) | proS | proline‑‑tRNA ligase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (10/9); total (11/9) | |||||||||||
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.95e-01 |
ACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGC > NZ_CP009273/214679‑214828 | aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 2:502037/90‑1 (MQ=255) ggCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGt < 1:589744/90‑1 (MQ=255) gcgcGG‑AATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt > 1:386223/1‑77 (MQ=255) gcgcGG‑AATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGAt < 2:386223/77‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 2:465032/90‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 1:245916/90‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa > 1:421106/1‑90 (MQ=255) aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 1:255417/75‑1 (MQ=255) aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa > 2:255417/1‑75 (MQ=255) aGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTaa > 1:571831/1‑90 (MQ=255) ttGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcg > 2:324727/1‑90 (MQ=255) gCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTcgcg < 1:579657/90‑1 (MQ=255) gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg > 1:4759/1‑90 (MQ=255) gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTAGCGGGAACGCATg > 1:623461/1‑90 (MQ=255) ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGaa > 1:1731/1‑77 (MQ=255) ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGaa < 2:1731/77‑1 (MQ=255) ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa > 1:464507/1‑90 (MQ=255) aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 2:604844/39‑1 (MQ=39) aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa > 1:604844/1‑39 (MQ=39) aTCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcg > 2:497515/1‑90 (MQ=255) tcgtcgGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAAcgc < 2:594708/90‑1 (MQ=255) | ACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAAACGC > NZ_CP009273/214679‑214828 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |