| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 4,180,586 | (A)7→6 | intergenic (+333/+77) | yjaZ → / ← thiH | heat‑shock protein C/2‑iminoacetate synthase ThiH |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,180,580 | 0 | A | . | 88.9% | 51.3 / 4.2 | 18 | intergenic (+327/+83) | yjaZ/thiH | heat‑shock protein C/2‑iminoacetate synthase ThiH |
| Reads supporting (aligned to +/- strand): ref base A (1/1); new base . (8/8); total (9/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.86e-01 | |||||||||||
TGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG > NZ_CP009273/4180493‑4180660 | tGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaa < 2:299606/90‑1 (MQ=255) cGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAcacaa < 1:328280/90‑2 (MQ=255) aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaaaacac < 1:264635/85‑4 (MQ=255) aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaaaacac > 2:264635/1‑82 (MQ=255) aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc < 2:276357/90‑1 (MQ=255) aaGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCtt < 1:14163/82‑1 (MQ=255) aaGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCtt > 2:14163/1‑82 (MQ=255) aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATc > 2:330859/1‑90 (MQ=255) tttCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTcaac > 2:387584/1‑60 (MQ=255) tttCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTcaac < 1:387584/60‑1 (MQ=255) tCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc < 2:302412/90‑1 (MQ=255) cGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCg > 1:380503/1‑90 (MQ=255) caTTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc < 2:239898/90‑1 (MQ=255) ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 2:85661/90‑1 (MQ=255) ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt < 2:335375/90‑1 (MQ=255) ccTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCTGAATAAATTTACCTTc > 2:381135/1‑90 (MQ=255) ccTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc > 1:7795/1‑90 (MQ=255) ccTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTc > 1:228063/1‑90 (MQ=255) tgACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa > 2:340652/1‑90 (MQ=255) cGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTTATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc > 2:18912/1‑90 (MQ=255) gCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTTATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg < 1:18912/90‑1 (MQ=255) | TGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG > NZ_CP009273/4180493‑4180660 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |