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 | 97.2% | 132.3 / ‑3.6 | 36 | Y174D (TAC→GAC) | proS | proline‑‑tRNA ligase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (16/19); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.72e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01 |
TAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAA > NZ_CP009273/214659‑214824 | tAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGttt < 1:174049/90‑1 (MQ=255) cGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCa < 1:504910/90‑1 (MQ=255) gTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgag > 1:502635/1‑90 (MQ=255) tACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga < 2:367447/90‑1 (MQ=255) aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 1:677794/90‑1 (MQ=255) aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 1:646816/90‑1 (MQ=255) aCGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 1:44262/90‑1 (MQ=255) gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa > 1:754092/1‑90 (MQ=255) gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa > 1:244565/1‑90 (MQ=255) gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa > 1:549619/1‑90 (MQ=255) gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa > 1:540020/1‑90 (MQ=255) gcgcGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTa > 1:507115/1‑90 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 2:781641/90‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 1:125544/90‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 2:549619/90‑1 (MQ=255) gAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGaa < 1:556935/90‑1 (MQ=255) cAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTa < 1:705270/90‑1 (MQ=255) ccccATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGc < 1:406390/90‑1 (MQ=255) aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt > 1:63388/1‑90 (MQ=255) aTGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCt > 2:778967/1‑90 (MQ=255) gAAGATTTTGCTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAg > 1:188200/1‑90 (MQ=255) cTGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAATCATCTTTCATCAGGAATTCGCgg < 1:1270/90‑1 (MQ=255) tGTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCggg < 2:754092/90‑1 (MQ=255) gTAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGa > 1:476668/1‑90 (MQ=255) tAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga < 1:125665/47‑1 (MQ=255) tAGGCCGCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTgaga > 2:125665/1‑47 (MQ=255) ccGCATACATTGCATCGTCGGTTTCCTGCAGGGAt < 2:341857/35‑1 (MQ=38) ccGCATACATTGCATCGTCGGTTTCCTGCAGGGAt > 1:341857/1‑35 (MQ=38) gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg > 1:207477/1‑90 (MQ=255) gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg > 1:307317/1‑90 (MQ=255) gCATACATTGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATg > 1:339745/1‑90 (MQ=255) ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa > 1:753554/1‑35 (MQ=38) ttGCATCGTCGGTTTCCTGCAGGGATTCCTGAGaa < 2:753554/35‑1 (MQ=255) ttGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGa > 1:451883/1‑90 (MQ=255) tGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaa < 2:502635/90‑1 (MQ=255) tGCATCGTCGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGaa < 1:778967/90‑1 (MQ=255) | TAGAACCGGTGTCGGCTTGTACGGCGCGGAAATCCAGCCCCATGCGGCTGAAGATTTTGCTGTAGGCCGCATACATTGCATCGTAGGTTTCCTGCAGGGATTCCTGAGAAGTATGGAAAGAGTAAGCATCTTTCATCAGGAATTCGCGGGAACGCATGACGCCGAA > NZ_CP009273/214659‑214824 |
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: ‑ |