Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 52.0% | 8.3 / 35.5 | 25 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (1/11); new base A (13/0); total (14/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-06 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380205‑2380478 | taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCg < 1:351452/139‑1 (MQ=255) aataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACg > 1:510350/1‑138 (MQ=255) tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACAACAAGcacaa > 1:309274/1‑138 (MQ=255) ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTcc < 1:196356/139‑1 (MQ=255) cagaCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCAt < 1:598479/136‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCAt < 2:596220/138‑1 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACAAAATGCACAT‑TTTATTTTTCCTCATCTTGTTTTagactg > 2:598479/1‑134 (MQ=255) tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAT‑TTTTTTTTTCCTCATC‑TGTTTTTAGATGCGAAACATGAAAATCAACAGAGGaaaaaa > 1:393439/1‑137 (MQ=255) gacaggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataa < 1:181674/135‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATggg < 2:223169/139‑1 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGGACAT‑TTTTTTTTTCCTCATC‑TGTTTTTAGATGCGGAAAAAGAAAATCAAAGGAGGAAAAATAAATGGGAATGATTGGc > 1:404949/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACAA‑TTTATTTTTCCTCATC‑TGTTTTAAGATGGGAAAAATGAAAATAAAAGGAGGGAAAAAAAAAGGGAATGATTggg > 1:124714/1‑137 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGACCAGGAAAAAAGACATG‑ACATCAATATTTTTTCCCATT‑TCGTTTTTGATACGAAAAAGAAAATTCCAAGGAGGATAAATAAATAGGAATGATTg > 2:542145/1‑136 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAATACACAA‑TTTTTTTTTCCTCCTC‑TGTTTTAAGAAGGGGAAAAAGAAAATTAAAGGGGGGAAAAAAAATGGGAATGATTGGCTATTTTg > 1:595303/1‑137 (MQ=255) aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAT‑TTTTTTTTTCCTCCTC‑TGTTTTTAGATGCGGAAAATGAAATTCAACGGAGGGAAAATAAATGGGGATGATTGGCTATTTTg > 1:210636/1‑136 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGCGCCCTAAAAAACGAGAAGCACAT‑TTTATTTTCCCTCCTC‑CGTTTTAAGACGGGGAACATGAAAATTAAACGAGGGAAAATAAAAGGGAAAGATTTGCTTTTTTGCAGaaaa > 1:339559/1‑138 (MQ=255) ataataAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc < 2:124714/139‑1 (MQ=255) aataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAT‑TTTTTTTTTCCTCATC‑TGTTTTAAGATGAGAAACATGAAATTTAACAGAGGGAAAAAAAATGGGGATGATTGGGTATTTTTCAGAAATTTACt > 2:522643/1‑138 (MQ=255) taaAAATTTCTCGAGGCGGGTCCCGAAAAAAACACATGCACAT‑TTTTTTTTTTCTCCTC‑TTTTTTTAGATGGGAAACATGAAATTTAACGGGGGGAAAAAAAATTGGGATGATTTGCTATTTTTTAGaaaattaat > 2:567674/1‑130 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAT‑TTTTTTTTTCCTCATC‑TGTTTTAAGATGCGAAAAATGAAATTTAACGGGGGGAAAAAAAATGGGGAAGATTGGGTATTTTTCAGAAAATGACt > 1:128769/1‑136 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa < 2:339559/138‑1 (MQ=255) gAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc < 1:223353/139‑1 (MQ=255) ggCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa < 2:128769/139‑1 (MQ=255) gTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCt < 1:522643/139‑1 (MQ=255) aGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 1:343290/1‑139 (MQ=255) gCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 2:302138/1‑138 (MQ=255) | TAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACAT‑TTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380205‑2380478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |