Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,841,672 | 0 | A | G | 57.1% | 27.5 / 83.1 | 63 | H198H (CAT→CAC) | selD | selenide, water dikinase SelD |
Reads supporting (aligned to +/- strand): ref base A (14/13); new base G (16/20); total (30/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.16e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCTCGGTCGGTACGATCCCCGTT > NZ_CP009273/1841535‑1841807 | gCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGAt > 2:297903/1‑139 (MQ=255) gCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGAt > 1:341625/1‑139 (MQ=255) cccTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGtt < 1:628448/139‑1 (MQ=255) cTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCt < 2:708431/139‑1 (MQ=255) tGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTg > 1:1222587/1‑139 (MQ=255) cGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAg > 1:74768/1‑139 (MQ=255) cTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACa > 1:774209/1‑136 (MQ=255) cTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACa < 2:774209/136‑1 (MQ=255) gCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTtc > 1:62507/1‑139 (MQ=255) cccAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCTCTTCCTTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTtct < 1:1229435/139‑1 (MQ=255) cccAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTtct < 2:405760/139‑1 (MQ=255) cccAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTtct < 2:1222587/139‑1 (MQ=255) cccAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTtct < 1:1074710/139‑1 (MQ=255) cagcagACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCa < 2:181098/139‑1 (MQ=255) aCCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg > 2:481471/1‑139 (MQ=255) ccAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTgg > 1:426994/1‑139 (MQ=255) cAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGCGCTTTTTTCTCAGCCGTgg > 2:711586/1‑138 (MQ=255) cccGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGaa > 2:495454/1‑139 (MQ=255) gTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGc > 2:1009782/1‑139 (MQ=255) gTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGc > 2:734962/1‑139 (MQ=255) tAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc > 2:1198090/1‑139 (MQ=255) ggTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATccc < 1:132783/138‑1 (MQ=255) gTCATCGCTTTTACGCCTTCGATGTTGGCCAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATccc > 2:1105832/1‑137 (MQ=255) cTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGttt > 1:28426/1‑139 (MQ=255) ttttACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTc > 2:856193/1‑139 (MQ=255) tACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGt > 2:332555/1‑138 (MQ=255) gCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAgg < 1:297903/139‑1 (MQ=255) ccTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAgg > 1:600762/1‑138 (MQ=255) tCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGAt < 1:735957/87‑1 (MQ=255) tCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGAt > 2:735957/1‑87 (MQ=255) tCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGaa > 2:78474/1‑137 (MQ=255) gATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGt < 2:1024417/139‑1 (MQ=255) aTGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGtt < 1:540472/139‑1 (MQ=255) aTGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGtt < 1:1014662/139‑1 (MQ=255) aTGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGtt < 1:332555/139‑1 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTtct < 1:412766/86‑1 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTtct > 2:338418/1‑86 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTtct < 1:338418/86‑1 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTtct > 2:412766/1‑86 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAg < 1:243999/89‑1 (MQ=255) ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAg > 2:243999/1‑89 (MQ=255) gCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATc > 2:348013/1‑139 (MQ=255) aGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGc < 2:62507/139‑1 (MQ=255) gcgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTtgtg < 1:1009782/139‑1 (MQ=255) gcgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTtgtg > 2:1179601/1‑139 (MQ=255) cgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGc > 2:88493/1‑139 (MQ=255) cgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGc < 1:778842/133‑1 (MQ=255) cgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGc > 2:778842/1‑133 (MQ=255) cgcCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGc > 1:344946/1‑139 (MQ=255) ccTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCgg < 2:57995/139‑1 (MQ=255) cGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCGGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTa < 2:341625/137‑1 (MQ=255) ttCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGtttt > 1:52680/1‑139 (MQ=255) cGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTttcttc < 1:811168/138‑1 (MQ=255) cGGCACATCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCa < 1:734962/139‑1 (MQ=255) aTCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCt < 1:795893/139‑1 (MQ=255) tCACTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCtc < 1:856193/139‑1 (MQ=255) cTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCtcggt < 1:88493/139‑1 (MQ=255) cTTCCGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCGCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCTGTACTGTTTTTCTTCACCCGCtcggt < 1:1007163/139‑1 (MQ=255) cTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCtcggt < 1:1198090/139‑1 (MQ=255) ccGTCGCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCtcggtc < 2:28426/137‑1 (MQ=255) gCCAGTCCCTGGTGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCTCGGTCGGTACGa > 2:791546/1‑139 (MQ=255) ccAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCTCGGTCGGTACGAt < 1:495454/139‑1 (MQ=255) cTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCTCGGTCGGTACGATCCCCGtt < 2:52680/139‑1 (MQ=255) | GCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGTCAGGAACAGTTTGCATCCGGCTTGTGCGGTACTGTTTTTCTTCACCCGCTCGGTCGGTACGATCCCCGTT > NZ_CP009273/1841535‑1841807 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |