| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,380,341 | 0 | G | A | 43.6% | 42.1 / 25.0 | 39 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
| Reads supporting (aligned to +/- strand): ref base G (6/16); new base A (17/0); total (23/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.06e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380220‑2380478 | cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCA‑GAAAAAACGACAAGCACaatttat > 1:759195/1‑139 (MQ=255) tGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGAGACATGCACattttattt > 1:199268/1‑139 (MQ=255) gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAACGACACTGAAattttttttt > 2:209359/1‑139 (MQ=255) ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTcc < 1:655290/139‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt < 1:878169/138‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt < 1:777924/138‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt < 2:296602/138‑1 (MQ=255) aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt < 1:290624/138‑1 (MQ=255) gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGACAAGAACATTTTTTTTTTTCT‑CATCTGTTTTa > 2:518702/1‑138 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAAAGACATGCACATTTTATTTTTCCT‑CATATGTTTTTAGAttc > 2:213475/1‑137 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CAt < 2:395881/125‑1 (MQ=255) aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCAAATTTTATTTTTCCT‑CAt > 1:395881/1‑125 (MQ=255) ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACa < 1:72753/139‑1 (MQ=255) aGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGaaa < 2:598949/139‑1 (MQ=255) gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACg < 1:518702/138‑1 (MQ=255) aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACACGCACATTTTATTTTTCCT‑CATCTGTTTTTAGAAACGAAACATGAAAATCAACggagg > 1:47943/1‑138 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAACGACATGCACAATTTATTTTTCCT‑CATCTGTTTTAAGAAACGAAACATAAAATTCCACGGAGGAAaaaaaaattg > 1:270828/1‑136 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGCTCCAGGAAAAAGGAAAATGC‑CCTTTTATTTTTCCTCCACTTTTTTTAAGTAGGGAAAAAGGAAA‑TCAACGGAAGAAAAAAAAATggg > 2:530469/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAACGGCATGCACATTTTATTTTTCCT‑CCTCTTTTTTTAGATGCGGAAAATGAAATTTAACCGGGGAAAAAAAAATAGGGATGATTTGCt > 1:34799/1‑139 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGACAAGCACATTTTTTTTTTCCT‑CATCTGTTTTTAGAAGCGAAAAATGAAAATTAACGGGGGAAAAATAAATTGGAATGAtttggt > 1:533725/1‑135 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGAAATGCACATTTTATTTTTCCT‑CATCTGTTTTTAGAAGCGAAAAAAGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc > 2:501294/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAAAGAAAAGAACATTTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAAAATGAAAATCAAAGGGGGAAAAAAAAATGGGAATGATTGGc > 2:186932/1‑138 (MQ=255) ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCA‑GAAAAAAAGAAATGAAAAATTTATTTTTCCT‑CACCTGTTTTTAGAAGGGAAAAAAAAAAATTAAAAAGGGGAAAATAAATGGGGATGAtttgg > 2:953469/1‑135 (MQ=255) tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACattttat < 2:504460/68‑1 (MQ=255) tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAGCGACATGCACatttttt > 1:504460/1‑66 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAAAGACAAGCACATTTTTTTTTCTCT‑CATCTGTTTTTAAAAGGGGAACAAGAAAATCAAAGGGGGGAAAAAAAATGGGGATGATTGGGTAttttt > 2:612492/1‑136 (MQ=255) aacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCC‑GAAAAAACGACATGCACAATTTTTTTTTCCT‑CCTCTGTTTTAAGATGCGAAACATGAAAATCAACGGGGGGAAAATAAAATGGAATGATTGGCTAttttt > 1:661717/1‑135 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt < 2:199268/139‑1 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt < 1:938245/139‑1 (MQ=255) tAAATAATAAAAATTTCTCGAGGCGGGGCCA‑GAAAAAACGACATGCACATTTTATTTTTCCT‑CATCTTTTTTTAGATGCGGAAAATGAAAATTAAAGGGGGGAAAAAAAATTGGGATGATTGGGTATTTTTTAGaaaat > 2:938245/1‑137 (MQ=255) aataatAAAAATTTCTCGAGGCGGGTGCA‑GAAAAAACGACATATACAATTTTTTTTTCCT‑CATCTGTTTTAAGATGAGAAACAAGGAATTTAAAGGGGGGAAAATAAATAGGGATTATTTGGTATTTTGCAGAAAATg > 1:796043/1‑138 (MQ=255) taaAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACATGCACATTTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAACAAGAAATTTAACGGAGGGAAAATAAATGGGAAAGATTGGGTATTTTGCAGAAAATGACt > 1:626291/1‑136 (MQ=255) taaAAATTTCTCGAGGCGGGTCCA‑GAAAAAACGACATGCACAATTTTTTTTTCCT‑CATCTGTTTTTAGATGCGAAAAAAGAAATTTAACGGGGGAAAAAAAAATGGGAATGATTGGCTATTTTTCAGAAAATGACt > 1:640669/1‑136 (MQ=255) aaaaTTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa < 1:381686/138‑1 (MQ=255) ctcGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAtt < 1:213475/138‑1 (MQ=255) cA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa < 1:501294/139‑1 (MQ=255) aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc > 1:220434/1‑135 (MQ=255) aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc < 2:220434/135‑1 (MQ=255) aGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 2:226319/1‑139 (MQ=255) gCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt > 1:693134/1‑138 (MQ=255) | CTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCA‑GAAAAAGCGACATGCACATTTTATTTTCCCT‑CATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACT > NZ_CP009273/2380220‑2380478 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |