Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 349,993 | 0 | T | G | 52.3% | 47.6 / 56.0 | 44 | V526G (GTT→GGT) | prpE | propionate‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base T (7/14); new base G (23/0); total (30/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATG > NZ_CP009273/349905‑350076 | gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGtt < 1:384416/90‑1 (MQ=255) aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgg > 2:116135/1‑90 (MQ=255) aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg > 2:88914/1‑88 (MQ=255) aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg > 1:417207/1‑88 (MQ=255) cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTgggggg > 1:150485/1‑88 (MQ=255) cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTggggggga > 1:170986/1‑86 (MQ=255) ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGaaaa > 1:234344/1‑89 (MQ=255) ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAg > 2:79262/1‑90 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCGGAAGTGGCGGTGGGTGGGGGGGAAGaa > 1:423385/1‑89 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt < 1:358085/90‑1 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAt > 1:65983/1‑90 (MQ=255) gCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAt > 1:85580/1‑90 (MQ=255) tCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAATGGGCTgga > 1:300808/1‑88 (MQ=255) gcgcgAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGGAAGGGGAgggggg > 2:18816/1‑89 (MQ=255) gcgcgAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGGAgggggg > 2:126343/1‑89 (MQ=255) gagaTTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtg < 1:116135/90‑1 (MQ=255) tGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGtt < 1:128876/90‑1 (MQ=255) gA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGTGGCGGGGGCGgtt > 1:391153/1‑90 (MQ=255) a‑‑aGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt > 1:174659/1‑76 (MQ=255) a‑‑aGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggt < 2:174659/76‑1 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCAt > 1:394461/1‑90 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCAt > 2:125873/1‑90 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTTGGGGGGAAGGTGGGCCGGAAGGGGAGGGGGGGGGGGGGGTTgggct > 2:201283/1‑86 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCAt > 2:96808/1‑90 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCAt > 2:427104/1‑90 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAAGCGCTGAAAGGGCGGGGGGGGGGGGGGGTTGtgct > 2:446470/1‑87 (MQ=255) gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGGAAGATGGGCTGAAAAGGGAGGGGGGGGGGGGGGTTTTCAt > 2:316228/1‑90 (MQ=255) gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCGGGGGGGGGGGGCGGTTTTCATTc > 1:350831/1‑90 (MQ=255) gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAAATGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCATTc > 2:28561/1‑90 (MQ=255) ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa < 1:18331/90‑1 (MQ=255) ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa < 1:99641/90‑1 (MQ=255) ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa < 1:308672/90‑1 (MQ=255) aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAagag > 1:166983/1‑90 (MQ=255) aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGGGCCGGAAGGGGAGGGGGGGGGGGCGGTTGTTGTTTCGAAagag > 2:285598/1‑90 (MQ=255) aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTGCAAGGGGAGGGGGCGGGGGGGGTTGTCATTTCGAAAGAGAGCg > 1:462129/1‑90 (MQ=255) tCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGAAAGATGCGCTGAAAAGGGAGGGGGGGGGGGGGGTTGTCCTTTCGCAAAAAAGCGa > 1:342777/1‑90 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 1:208745/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 1:245441/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 2:242035/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 1:271700/90‑1 (MQ=255) ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc < 2:394461/90‑1 (MQ=255) gggCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg < 2:279210/90‑1 (MQ=255) ccGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg < 1:316228/90‑1 (MQ=255) cGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCgtggtg > 1:264671/1‑90 (MQ=255) cGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtg > 2:248511/1‑90 (MQ=255) | GTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGA‑‑AGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATG > NZ_CP009273/349905‑350076 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |