Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,617,873 | 0 | A | T | 37.1% | 48.4 / 19.6 | 35 | I491F (ATT→TTT) | ppk1 | polyphosphate kinase 1 |
Reads supporting (aligned to +/- strand): ref base A (5/17); new base T (13/0); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.42e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAG > NZ_CP009273/2617738‑2618010 | aaCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTg < 1:444967/139‑1 (MQ=255) cGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGaa < 2:796215/139‑1 (MQ=255) tggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATa < 2:921604/139‑1 (MQ=255) gtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAAATTTATTGAAAACCCATAcc > 1:843860/1‑139 (MQ=255) ggtcGTTTTGCACACACCGGGCCCGGGAACTTTAACAAAAAACCCGCGCTTCTTTATACTTCCTATTCGTTGCTTACCGCCGATGCCCTCATCACCAACGAAGTACGTCGGGTA‑TTTAACTTTATTGAAAACCCATAcc < 2:576184/136‑1 (MQ=255) tGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGA‑TTTTACTTTTTTTAAAAACCATACCGTCCCGTg > 1:463606/1‑139 (MQ=255) acacaTCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGA‑TTTTAATTTTTTTAAAAACCCTACCGTCCGGGGACa > 1:260543/1‑139 (MQ=255) aTCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTTACTTTTTTTAAAACCCATACCGCCCGGTGGCATTTg > 2:888267/1‑139 (MQ=255) aTCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTTTAAAACCC‑TACCCTCCGGTGaaacttt > 2:616950/1‑133 (MQ=255) tCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGa < 1:427048/139‑1 (MQ=255) gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTG‑TTTAACTTTTTTTAAAACCCACACCCTCCGGGGACATTTTatt > 2:940548/1‑139 (MQ=255) ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAAATTTTTTGAAAAACCCTACCCTCCGGGGACATTTTtttattttat > 1:558858/1‑139 (MQ=255) aaaaCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGT‑TTTTACTTTATTTTAAACCCATACCGTCCGGGGACATTTGATTATTTTATGGTATcgccgcca > 1:614064/1‑136 (MQ=255) cgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGGGGCATTTGATTATTTTATGGTATCGCCGCAAAACAcccc > 1:228407/1‑138 (MQ=255) gcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTG‑TTTAACTTTTTTTAAAAACCACAACGCCCGGTGGCCTTTGATTTTTTAATTGTATTGCCGCAAAACTCccgc > 2:888552/1‑139 (MQ=255) tCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCcgccgc < 2:843860/139‑1 (MQ=255) tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTa < 1:838567/138‑1 (MQ=255) tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTa < 1:888267/138‑1 (MQ=255) tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTa < 2:187279/138‑1 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTTACTTTTTTGTAAACCCATACCCTGCGGGGGCCTTTGTTTATTTTATTGTGTTGCCCCCAAAATCCCGCCGCCtttttgatg > 2:120649/1‑133 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTTACTTTTTTGAAAACCCACACCCTCCGGGGACATTTTATTTTTTTATGTGATCGCCGCCAAAATCCCCCCGCCTATTGTTTg > 2:169245/1‑138 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATg < 2:560092/138‑1 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGA‑TTTTACTTTTTTTAAAACCCATACCGTCCGGGGACATTTGATTTTTTTATGGTATTGCCGCAAAAATCCCGCCGCCCATTGTTTg > 1:52905/1‑138 (MQ=255) gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTG‑TTTTACTTTTTTGAAAACCCCTACCCCCCGGGGGGCTTTTTTTTTTTTATTGGATTCCCGCAAAAATCCCCCCGCCCCTTTTTTGAAATTGGGGAc > 2:865092/1‑139 (MQ=255) gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTTACTTTTTTTTAAAACCCTACCCTCCGGGGGCCTTTTTTTATTTAATGGGATAGCCGCAAAAATCCCCCCCCCCCTTTTTTGAAAaggggggc > 2:554073/1‑137 (MQ=255) cgccgATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCgag < 1:169245/139‑1 (MQ=255) aCCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCgcagc < 2:52905/139‑1 (MQ=255) cGAAGTACGGCGGGGA‑TTTTACTTTATTTAAAAACCATACCGCCCGGGGAAATTTTATTATTTAATGGGATCGCCGCAAAACTCCCCCCCCCTCTTGTATGAAAAGGTGGACCGCGAGATCGCCAACCCCCAGCAAggg > 1:431267/1‑139 (MQ=255) cGAAGTACGGCGGGGA‑TTTAAATTTTTTGGAAACCCCTACCCTCCGGGGGCATTTTGTTATTTTATGGGATCGCCGCCAAAATCCCCCCCCCCATTGTATGAAATGGTGGACCGCGGGATAGCCAACCCCCAGCAAggg > 1:150493/1‑139 (MQ=255) gAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGc < 2:463606/139‑1 (MQ=255) gTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGc < 2:558858/139‑1 (MQ=255) ggcggGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAg < 2:150493/139‑1 (MQ=255) tA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTAt < 2:516295/139‑1 (MQ=255) tA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTAt < 2:797627/139‑1 (MQ=255) tATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAg > 1:611238/1‑139 (MQ=255) | AACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTA‑TTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAG > NZ_CP009273/2617738‑2618010 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |