Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 51.8% | 7.4 / 36.3 | 27 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (4/9); new base A (14/0); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380231‑2380478 | aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt < 1:402174/138‑1 (MQ=255) cagacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGtttt < 2:357240/139‑1 (MQ=255) gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAtt > 2:320967/1‑109 (MQ=255) gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATg < 1:320967/110‑1 (MQ=255) acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACCGCAAGCACATTTTTTTTTTCCTCCACCGttttta > 2:402174/1‑136 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc < 1:49813/106‑1 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc > 2:49813/1‑106 (MQ=255) cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGGAACATGa > 1:380630/1‑138 (MQ=255) aGGTAGTCGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGTCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCGCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagg < 2:412728/138‑1 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAAAAAAAACAACATGCAAATTTTATTTTTCCTCATCTGTTTTTAGATGGGAAACATGAAAATTAAAGGAGGgaaaa > 2:317648/1‑138 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCACATTTTATTTTTCCTCATCTGTTTTTAGAAGCGGAACAATAAATTTAACGGAGGAAAAATAAATtgg > 1:342737/1‑139 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAACATAAAAATTAAAGGGGGGAAAAAAAATggg > 2:335007/1‑139 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCACATTTTATTTTTCCTCCTATTGTTTAAGAAGCGAAAAATAAAAATTAAAGGGGGAAAAAAAAATggg > 1:100015/1‑138 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTATTTTTCCTCATCTGTTTTAAGATGGGAAACATGAAAATCAAAGGGGGGAAAATAAATggg > 1:299808/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGGAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGGAACATTAAATTTAACGGAGGGAAAATAAATTGGGATGATTgggt > 2:80006/1‑137 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGGCATTCCCAATTTATTTTTCCTCACCTGTTTTTAAATGCGAAACCAGAAATTCAAACGGGGAAAAATAAAAGGGGATGGTTTGCt > 1:445839/1‑139 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt < 1:335007/139‑1 (MQ=255) tAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGAACATTTTTTTTTTCCTCAACTGTTTTTAGATGCGGAAAATGAAATTTAACGGAGGGAAAAAAAAAGGGGATGATTTGCTATTTTGCAGAAAtt > 1:312782/1‑139 (MQ=255) aataatAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACAATTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAAAGAAATTCAAAGGAGGAAAAATAAATGGGGGTGATTGGCTTTTTTTTAGAAAtag > 2:353274/1‑136 (MQ=255) aataatAAAAATTTCTCGAGGCGGGGCCCGAAAAAAACACAAGAAAATTTTTTTTTTTCCCATCTGGTTTTAAATTGGAAACAATAAAATACACGGAGGAAAAATAAATGGGGAAGATTTGCAATTTTTTAGAAAATg > 2:178805/1‑138 (MQ=255) aataaAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTAAGATGCGAAAAATGAAATTCAAAGGAGGGAAAAAAAATGGGGATGATTGGGTATTTTGCAGAAATTGaat > 1:349050/1‑136 (MQ=255) taaAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTATTTTTCCTCCTCTGTTTTTAGATGCGAAAAATGAAAATTAAAGGAGGGAAAAAAAATTGGAATGTTTTGCTATTTTTCAAAAAATTACt > 2:107015/1‑136 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa < 1:178805/138‑1 (MQ=255) aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa < 2:366617/138‑1 (MQ=255) ggCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa < 2:349050/139‑1 (MQ=255) ggTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGAAACATGAAATTCAACGGAGGGAAAATAAATGGGAATGATTGGCTATTTTTTAGAAATTTACTCCGAAAAAATTAAACAACTGc > 2:8226/1‑139 (MQ=255) aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 2:494280/1‑139 (MQ=255) gCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa > 2:590512/1‑131 (MQ=255) gCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa < 1:590512/131‑1 (MQ=255) | ATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380231‑2380478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |