Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,717,969 | T→C | E237G (GAA→GGA) | glyQ ← | glycine‑‑tRNA ligase subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,717,969 | 0 | T | C | 86.8% | 86.3 / 8.2 | 38 | E237G (GAA→GGA) | glyQ | glycine‑‑tRNA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base T (3/2); new base C (17/16); total (20/18) | |||||||||||
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.22e-01 |
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑3718051 | gaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgcgc > 2:301853/1‑90 (MQ=255) tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 1:199058/1‑90 (MQ=255) tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca > 1:653015/1‑90 (MQ=255) ttACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCagcagc > 1:78719/1‑90 (MQ=255) gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCtgctgc > 2:781733/1‑90 (MQ=255) cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATGCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:718390/90‑1 (MQ=255) cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:810707/90‑1 (MQ=255) gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 1:551131/1‑84 (MQ=255) gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:551131/84‑1 (MQ=255) gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTtctc > 1:801587/1‑90 (MQ=255) aaGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAc < 2:653015/90‑1 (MQ=255) gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:631718/1‑90 (MQ=255) gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 2:67112/1‑90 (MQ=255) gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc > 1:689769/1‑90 (MQ=255) gggCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGa > 1:387456/1‑90 (MQ=255) cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc > 1:523339/1‑90 (MQ=255) cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc < 2:517303/90‑1 (MQ=255) gcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAg > 1:163107/1‑90 (MQ=255) ttCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTg > 2:327349/1‑90 (MQ=255) aaTACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAg < 2:501346/90‑1 (MQ=255) tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 2:288569/90‑1 (MQ=255) tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa < 2:10334/90‑1 (MQ=255) gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc > 1:631754/1‑71 (MQ=255) gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc < 2:631754/71‑1 (MQ=255) gTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCcac > 2:712828/1‑90 (MQ=255) aGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACAt < 2:689769/90‑1 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 2:340245/1‑44 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 2:445254/44‑1 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc < 1:340245/44‑1 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc > 1:445254/1‑44 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa > 2:825222/1‑73 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa < 1:825222/73‑1 (MQ=255) ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg > 1:121735/1‑90 (MQ=255) gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 1:298582/90‑1 (MQ=255) gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa < 2:265751/90‑1 (MQ=255) cggATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 2:523339/90‑1 (MQ=255) cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 1:508479/90‑1 (MQ=255) cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt < 2:45414/90‑1 (MQ=255) | GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT > NZ_CP009273/3717887‑3718051 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |