| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 3,981,822 | A→G | I202T (ATC→ACC) | hemX ← | uroporphyrinogen‑III C‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,981,822 | 0 | A | G | 88.9% | 94.6 / 5.6 | 36 | I202T (ATC→ACC) | hemX | uroporphyrinogen‑III C‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (12/20); total (15/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.87e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01 | |||||||||||
TGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGTCATTCATATCCGCCAGGCTGGCGTCTGCACTTTTCAGCAACGCTGCAG > NZ_CP009273/3981687‑3981954 | tGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatggtg < 2:301551/138‑1 (MQ=255) aggttGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGt < 2:224744/136‑1 (MQ=255) gATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGt > 2:401811/1‑138 (MQ=255) ttCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATa < 1:299934/129‑1 (MQ=255) ttCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATa > 2:299934/1‑129 (MQ=255) cGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCg < 1:187785/138‑1 (MQ=255) cTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATa < 2:383052/139‑1 (MQ=255) cTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATa < 1:560424/139‑1 (MQ=255) aTGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATAc < 2:521152/137‑1 (MQ=255) aTGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATAc < 1:452947/137‑1 (MQ=255) cTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGaaa < 2:66980/139‑1 (MQ=255) gAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc < 2:375557/139‑1 (MQ=255) ctctTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGc > 1:156825/1‑138 (MQ=255) ttCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACG‑AGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGatat > 1:542673/1‑139 (MQ=255) cgccGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGat < 2:459297/134‑1 (MQ=255) ccGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATatcat < 1:278784/139‑1 (MQ=255) cGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATcc > 2:545800/1‑104 (MQ=255) cGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATcc < 1:545800/104‑1 (MQ=255) tCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCgg > 2:509034/1‑139 (MQ=255) aaTCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGt < 1:185905/136‑1 (MQ=255) tCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGat < 2:545075/121‑1 (MQ=255) tCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGat > 1:545075/1‑121 (MQ=255) tCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGAc > 1:301666/1‑139 (MQ=255) aTTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGc < 1:448547/62‑1 (MQ=255) aTTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGc > 2:448547/1‑62 (MQ=255) tCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCggg < 1:252338/138‑1 (MQ=255) tCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCggg < 2:544561/138‑1 (MQ=255) tCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCggg < 1:432150/138‑1 (MQ=255) ccAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGtcatt > 1:316788/1‑139 (MQ=255) gACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGtcattcat < 2:301666/139‑1 (MQ=255) gATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCGCCTg < 2:536923/49‑1 (MQ=255) gATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTg > 1:536923/1‑49 (MQ=255) gATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGTCATTCATATCCGCCAGGCTGGCGTc > 1:521862/1‑139 (MQ=255) ttGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGaa > 2:37174/1‑93 (MQ=255) ttGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGaa < 1:37174/93‑1 (MQ=255) aTTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGTCATTCATATCCGCCAGGCTGGCGTCTGCACTTTTCAg > 2:42490/1‑138 (MQ=255) aGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGTCATTCATATCCGCCAGGCTGGCGTCTGCACTTTTCAGCAACGCTGCAg > 2:255192/1‑139 (MQ=255) | TGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGTCATTCATATCCGCCAGGCTGGCGTCTGCACTTTTCAGCAACGCTGCAG > NZ_CP009273/3981687‑3981954 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |