Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,819,318 | 0 | G | A | 49.1% | 37.6 / 41.2 | 57 | intergenic (‑124/+79) | ves/spy | environmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy |
Reads supporting (aligned to +/- strand): ref base G (27/2); new base A (0/28); total (27/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.09e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTC > NZ_CP009273/1819230‑1819402 | tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga > 2:692284/1‑90 (MQ=255) aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa > 2:2432667/1‑90 (MQ=255) aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa > 2:2930825/1‑90 (MQ=255) aaCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg > 1:730003/1‑90 (MQ=255) ctgTGGATGGTTATCGGCAACTCGGGTGATTTCCGTTAGTGGTGATCAGGCTTTTTTTACAATTGTTTGAGAAGAGAAAAAAAAAccgcc < 2:524461/88‑1 (MQ=255) ccGGGGATGTTTATCGGCAACTGGGGTGATTTACGTTGGTGGTGATCAGTTTTTTTTTACAATTGTTGGAAAAGAAAAAAAAAAAccgcc < 1:558386/90‑1 (MQ=255) gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa < 2:1125478/90‑1 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc > 2:1932180/1‑90 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc > 1:2941543/1‑90 (MQ=255) ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc > 2:2974494/1‑90 (MQ=255) ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc > 2:1271029/1‑90 (MQ=255) aTCGGCAACTCGGGTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:1357408/90‑1 (MQ=255) tCGGCAACTGGGGTGATTTACTTTAGGGGTGATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc < 1:358535/90‑1 (MQ=255) tCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc < 1:776412/90‑1 (MQ=255) cAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAt > 2:1030710/1‑90 (MQ=255) aaCTGGGGTTTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCAtt < 1:2362734/90‑1 (MQ=255) aaCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt < 1:824480/90‑1 (MQ=255) aaCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt < 2:2767955/90‑1 (MQ=255) aaCTGGGGTAATTTACGTGAGTGGTAATCGGGTTTTTTTAACAATTTTTTAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCAtt < 1:2418798/90‑1 (MQ=255) aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt < 2:757418/90‑1 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:526984/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:2362734/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:2838084/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:259730/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc > 1:867697/1‑90 (MQ=255) ctggggTGTTTTACGTTGGTGGTGATCGGGTTTTTTTTACAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTAc < 1:1061839/87‑1 (MQ=255) gtgtTTTATGTTGGGGGTGATCGGTTTTTTTTAACATTTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 1:2548795/90‑1 (MQ=255) gtgtTTTACTTTAGTGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 1:692284/90‑1 (MQ=255) gtgtTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:1182693/90‑1 (MQ=255) gtgtTTTACGTTAGGGGTGATCAGGTTTTTTTTAAAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 1:94105/90‑1 (MQ=255) gTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:1673838/90‑1 (MQ=255) gTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 1:815763/90‑1 (MQ=255) gTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:958974/90‑1 (MQ=255) aTTTACTTTAGGGGGGCTCGGTTTTTTTTTACATTTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt < 1:2785927/90‑1 (MQ=255) tttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa > 2:2452068/1‑90 (MQ=255) ttACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:623507/90‑1 (MQ=255) ttACGTTAGTGGTGATAAGTTTTTTTTTACAACTGTTTAAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:1534609/90‑1 (MQ=255) ttACGTAAGGGGTGATCGGTTTTTTTTTAAAACTTTTGGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 1:1797722/90‑1 (MQ=255) gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:3002695/1‑90 (MQ=255) tAGTGGTAATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa < 2:2997092/90‑1 (MQ=255) aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa > 1:1279706/1‑90 (MQ=255) aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa > 2:2304114/1‑90 (MQ=255) gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg > 2:1660491/1‑90 (MQ=255) gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg > 2:681536/1‑90 (MQ=255) gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg > 1:469151/1‑90 (MQ=255) aTCGGGTTTTTTTCACACTTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc < 1:1271029/90‑1 (MQ=255) aTCAGGTTTTTTTTACAACTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc < 2:2769883/90‑1 (MQ=255) tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 1:1114634/1‑90 (MQ=255) aGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCaa > 1:59347/1‑68 (MQ=255) aGGCTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:59347/68‑1 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggc > 1:2138021/1‑88 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggc > 1:2558753/1‑88 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc > 1:1796345/1‑90 (MQ=255) tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg > 2:1214927/1‑90 (MQ=255) aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt > 1:931225/1‑90 (MQ=255) gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg > 2:507277/1‑90 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa > 1:2403517/1‑90 (MQ=255) tGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa < 2:2168285/90‑1 (MQ=255) gaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg < 2:732024/87‑1 (MQ=255) aaagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat < 2:2200774/89‑1 (MQ=255) gaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat < 1:1263588/87‑1 (MQ=255) aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat < 1:2974494/87‑1 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt > 2:2070036/1‑90 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt > 2:821095/1‑90 (MQ=255) acaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt > 2:2349958/3‑90 (MQ=255) aaaaaaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc < 1:2991485/89‑1 (MQ=255) | TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTC > NZ_CP009273/1819230‑1819402 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |