Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,213,525 | 0 | C | T | 42.9% | 45.1 / 24.0 | 42 | R11K (AGA→AAA) | iclR | glyoxylate bypass operon transcriptional repressor IclR |
Reads supporting (aligned to +/- strand): ref base C (1/23); new base T (18/0); total (19/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.37e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG > NZ_CP009273/4213439‑4213613 | aTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGc > 1:1012886/1‑90 (MQ=255) tCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcg > 1:1221747/1‑90 (MQ=255) tCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGccc > 1:1332586/1‑89 (MQ=255) ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgcg > 2:807829/1‑90 (MQ=255) aTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGtttt > 2:995359/1‑89 (MQ=255) cAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt < 2:1491644/90‑1 (MQ=255) aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATg < 2:1349463/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 2:218800/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 1:529276/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 1:562564/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 1:874291/90‑1 (MQ=255) ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCg < 1:491012/90‑1 (MQ=255) ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGTTTCGCCGGGAgggggggg > 2:1426203/1‑86 (MQ=255) cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACa < 2:466738/90‑1 (MQ=255) cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 1:172149/90‑1 (MQ=255) ttAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct < 2:926870/90‑1 (MQ=255) tAAAGACTGAAACTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc < 1:775032/90‑1 (MQ=255) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt < 2:952483/90‑1 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 2:1012886/90‑1 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 2:1529365/90‑1 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGGTTTGCGGGGATGGGGGGGGCCATGACAGGCTCttttttt > 2:1256057/1‑89 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAGGGGGGGGGCCCTGACAGGCTCttttttt > 2:1466011/1‑89 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGCTTCGCGGGAATGGGGGGGCCCCTGGAAGGCTCttttttt > 1:167325/1‑89 (MQ=255) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGGGGCCCCGACAGTCTCtttttttt > 2:1533599/1‑88 (MQ=255) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGGGGCCCCGACAGGCTCtttttttt > 2:918649/1‑88 (MQ=255) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGACCCCGACAGTCTCtttttttt > 1:911242/1‑88 (MQ=255) aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGTTTTGGGGGGAAGGGGGGGGGCCTGAGAGGCTCttttttttt > 2:1142537/1‑87 (MQ=255) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTAt < 1:1533599/90‑1 (MQ=255) tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATc < 2:1326411/90‑1 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTCCCCCGCTTCGCGGGGAAGGGGGGGGCCATGACAGCCCCTTTTTTTTGttttt > 1:662760/1‑86 (MQ=255) cGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGTTTTGCGGGGATGGGGGGGGCCCCGACAGGCTCTTTTTTTTTTTTCGTGGgaaa > 2:214125/1‑89 (MQ=255) tGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGATGGGGGGGGCCATGATAGTGTCTTTTTTTTGTTTTGTGGGAAtcat > 2:383184/1‑90 (MQ=255) gTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGATGGGGGGGGCCCTGACAGGCTCTTTTTTTTTTTTTGTGGGAAtcattt > 2:826875/1‑90 (MQ=255) gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCTACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcattttc < 1:568224/90‑1 (MQ=255) cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGGGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAAtcattttta > 1:412328/1‑88 (MQ=255) cAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGATGGGGGGGGCCATGACAGGCTCTTTTTTTTTTTTCGTGGGAAATATttttttttttt > 2:170217/1‑89 (MQ=255) ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt < 2:1332586/90‑1 (MQ=255) ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt < 1:269014/90‑1 (MQ=255) gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt < 2:672470/90‑1 (MQ=255) gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt < 2:644944/90‑1 (MQ=255) gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt < 1:1112016/90‑1 (MQ=255) tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg < 2:302021/90‑1 (MQ=255) | ATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATG > NZ_CP009273/4213439‑4213613 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |