Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,371,314 | 0 | G | A | 58.3% | 2.9 / 27.7 | 24 | F475F (TTC→TTT) | menD | 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase |
Reads supporting (aligned to +/- strand): ref base G (10/0); new base A (0/14); total (10/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.10e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATCGTAAA > NZ_CP009273/2371225‑2371402 | tCGAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCgag > 2:4316/1‑90 (MQ=255) cgcgGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTg > 2:1157248/1‑90 (MQ=255) ggcggcgTGCTAAAAAGGGCCGTTTTGCGGCATAAGAAAAAAACGCCCCCGCTCGTTTTGCGGTGTGGGCAAAAGCGAAAAAATTTGccc < 1:1537233/90‑1 (MQ=255) gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt > 1:825012/1‑90 (MQ=255) aaaTGGCCTTTTTGGGGCATAAAAAAAAACCGCACCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTc < 1:1701493/90‑1 (MQ=255) aaaTGGACTTTTTGCGCCATAAAAAAAAACCCCCCCCCCTCTCTTTGCGGCGTTGCCAACAGCAAAAAATTTTCCCCGCCGTTGTTGTTc < 1:590237/90‑1 (MQ=255) aaaTGGACGTTTTGGGGAAAAAAAAAAAAACCCCCCCGCTCGTTTTGCGGCGTTGGCAACAGAAAAAAATTTTCCCCGCCGTTGTTGTTc < 1:1155678/90‑1 (MQ=255) ggACTTTTGGCGGCATAAGATAAAACCGCTCACGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacca < 2:825013/90‑1 (MQ=255) gccGTTTTGCGCCATCAGATAGAAACGCTCCCGCTCGCTTTGCGGCGTTGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 2:1145416/88‑1 (MQ=255) gatgTTTTGCGGCATCAGATAAAACCCCCCCCCCTCTTTTTTCGGCGTGGGCAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTcaccac < 1:172950/87‑1 (MQ=255) gACTTTTTGCGGCAAAAGAAAAAACCCCCCCCCCTCGTTTTGCGGCTTTGGAAACAGGAAAAAAATTTCCCCGCCGTTGTTGTTcaccac < 1:525478/90‑1 (MQ=255) gACTTTTGGCGGCATCAGATAAAAACGCCCCCGCTCGTTTTGGGGCTTGGGAAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTcaccac < 1:4316/90‑1 (MQ=255) gACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 1:935533/90‑1 (MQ=255) gACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac > 1:763610/1‑90 (MQ=255) gCGGCATCAGATAGAAACGCTCACGCTCGCTTTTCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATa > 2:172950/1‑90 (MQ=255) tCAGAAAGAACCGCTCCCCCCCGTTTTGCGGCGTGGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgcg < 2:776633/90‑1 (MQ=255) gAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGa > 2:1221163/1‑90 (MQ=255) gAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGa > 2:167709/1‑90 (MQ=255) gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt > 1:79648/1‑90 (MQ=255) cgttttGCGGCGTGGCAAACAGGAAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATa < 2:1683834/87‑1 (MQ=255) ccgtttGGGGGGTTGGCAACAGGAAAAAAATTTTCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATa < 2:902074/87‑1 (MQ=255) gTTGGAAACAGCGAAAAAATTCCCCCGCCGTTGTTGTTCACCACATTTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGcgcg < 1:1157247/90‑1 (MQ=255) ggCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTg > 2:188675/1‑90 (MQ=255) cAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTgag > 1:1187898/1‑90 (MQ=255) gCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATCGt > 2:905521/1‑90 (MQ=255) agaAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATCGTaaa < 1:905520/90‑1 (MQ=255) | TCGAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATCGTAAA > NZ_CP009273/2371225‑2371402 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |