Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,367,400 | 0 | A | C | 54.2% | 27.9 / 50.7 | 48 | G361G (GGT→GGG) | menE | o‑succinylbenzoate‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base A (14/8); new base C (0/26); total (14/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.63e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGC > NZ_CP009273/2367316‑2367486 | cTCCTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCctc < 2:1108916/90‑1 (MQ=255) cTTGTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTcc < 1:443528/90‑1 (MQ=255) cTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc > 1:1426935/1‑90 (MQ=255) gTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGcc < 1:1312030/90‑1 (MQ=255) cGGCAACGGGGACGATAACCCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 2:395315/90‑1 (MQ=255) gCAACGGGGACGATAACCACCTGCAAACCCGCAGGATGTGCACCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 2:389749/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 2:1456378/90‑1 (MQ=255) gCAACGGGAACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 1:345894/90‑1 (MQ=255) ggACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATa < 1:1016073/90‑1 (MQ=255) cGATAAACCCCTGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 2:132044/90‑1 (MQ=255) cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 2:1174059/90‑1 (MQ=255) cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 1:770744/90‑1 (MQ=255) ccctcAAAACCGCCCAATGTGCCGCAATTACGGGCGCGCCTTCTTCCGCTTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc < 2:1023260/86‑1 (MQ=255) tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt > 1:322750/1‑90 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAc < 2:487244/54‑1 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAc > 1:487244/1‑54 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc > 1:1602715/1‑90 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc > 2:1425934/1‑90 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc > 1:1350304/1‑90 (MQ=255) cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 2:302989/90‑1 (MQ=255) cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 1:1547060/90‑1 (MQ=255) aCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACa > 1:43026/1‑90 (MQ=255) cAGGATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATgg < 2:17559/90‑1 (MQ=255) ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 1:403693/90‑1 (MQ=255) ggATGTGCAGAAATTCCGCGCTCGTCTTCTTCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 2:1470887/90‑1 (MQ=255) ggATGTGCACCAATTACGCGCTCGCCTTCTTCCGGCTAAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 2:1350304/90‑1 (MQ=255) tgCAGCATTTCCGCGCTGGCCTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt < 2:745816/90‑1 (MQ=255) tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt < 1:1549870/90‑1 (MQ=255) cagcaATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGc < 1:1085966/90‑1 (MQ=255) caATTACGCGCTCGACTTTTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt < 1:1399273/90‑1 (MQ=255) aTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa > 2:623654/1‑90 (MQ=255) tACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt < 1:147734/81‑1 (MQ=255) tACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt > 2:147734/1‑81 (MQ=255) tACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg > 2:722932/1‑90 (MQ=255) gcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 2:1415940/36‑1 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc > 1:1415940/1‑36 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg > 2:1418280/1‑38 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 1:1418280/38‑1 (MQ=255) cTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtctc < 1:722932/90‑1 (MQ=255) cTCGACTTCTTCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtctc < 2:532351/90‑1 (MQ=255) tttcttcCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc < 2:34857/89‑1 (MQ=255) tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGa < 1:659054/90‑1 (MQ=255) tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGa < 1:397780/90‑1 (MQ=255) tcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGa < 1:1425934/90‑1 (MQ=255) gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 2:177095/90‑1 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 2:1602715/90‑1 (MQ=255) aaTACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATTCATCTCACCGCGATCGCgcgtagc > 1:1352831/1‑90 (MQ=255) aaTACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtagc > 1:448962/1‑90 (MQ=255) | CTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAGC > NZ_CP009273/2367316‑2367486 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |