Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,367,400 | 0 | A | C | 53.6% | 17.8 / 31.3 | 28 | G361G (GGT→GGG) | menE | o‑succinylbenzoate‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base A (6/7); new base C (0/15); total (6/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGT > NZ_CP009273/2367313‑2367483 | aaaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAca > 1:585194/1‑89 (MQ=255) aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc < 1:229801/90‑1 (MQ=255) aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc < 1:571294/90‑1 (MQ=255) cTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTcc < 2:53893/90‑1 (MQ=255) cTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc > 1:687542/1‑90 (MQ=255) gTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGcc < 2:776104/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 2:792028/90‑1 (MQ=255) cAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTTTTCCGGCTGAATCCCCTCTCCGCCACTga < 1:735774/90‑1 (MQ=255) cGGGGACGATAACCCCCTGCAAAACCCCAGCATGTCCACCAATTACGCGCTCGTCTTCTTCCGGCTGAATACCCTCTCCGCCACTgaaga < 1:699422/90‑1 (MQ=255) cGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaaga < 1:301536/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 1:587913/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 1:388720/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 2:553572/90‑1 (MQ=255) ccTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc < 2:628678/90‑1 (MQ=255) tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt < 2:685197/90‑1 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg > 2:489507/1‑68 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 1:489507/68‑1 (MQ=255) cAAAACCCCAGGATGTGCACCAATTACGCGCTCGATTTCTTCCGGCCGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 1:599631/90‑1 (MQ=255) ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 2:165749/90‑1 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt > 2:188875/1‑90 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt > 2:481142/1‑90 (MQ=255) tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt < 2:28745/90‑1 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt < 1:684909/90‑1 (MQ=255) ccgccCTCGACTTCTCCCGGCTAAATCCCCCCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGc < 1:248662/89‑1 (MQ=255) ttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc < 1:188875/90‑1 (MQ=255) cGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcg > 1:306949/1‑90 (MQ=255) ggCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcgc < 2:687542/90‑1 (MQ=255) cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt < 2:240139/90‑1 (MQ=255) | AAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGT > NZ_CP009273/2367313‑2367483 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |