Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 293,763 | 0 | C | A | 25.7% | 174.4 / 11.8 | 70 | L225F (TTG→TTT) | paoD | molybdenum cofactor insertion chaperone PaoD |
Reads supporting (aligned to +/- strand): ref base C (32/20); new base A (0/18); total (32/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.39e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCT > NZ_CP009273/293675‑293847 | gATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCa > 2:144128/1‑90 (MQ=255) cccAATGCGCCGAGAAAAAAGGGTTTTGCTTCCCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGAcc < 2:164670/90‑1 (MQ=255) aaTGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACGGGCAGCTCCCGGTTGAGATCATGGCACAGAAAAATGAccgcc < 2:235324/90‑1 (MQ=255) aaTGCGCCGAGATAAAAGGGTTTTGCTTCCCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGAccgcc < 1:111980/90‑1 (MQ=255) cgcCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCCCCCGGTTGAGATCATGGCACA‑AAAAATGACCGCCgtatc < 1:451580/90‑1 (MQ=255) gcggagaTAAAGGGGTTTTGCTTCCCGCGCGGCCGGCAACCCTGGCACCCCCAGGTTGATATCATGGCCCAGCAAAATGACCGCCgtatc < 2:364280/87‑1 (MQ=255) gcCGAGATAAAGGGTTTTTGCTTCTCGCGGGGCCTGCACCCCGGGCAGCCCCGGGTTGAGATCAGGGCACAGCAAAATGACCGCCgtatc < 2:487984/90‑1 (MQ=255) gcCGAGATAAAGGGGTTTTGCTTCCCGCGCGGCCTGCAACACTGGCAGCCCCCGGTTGAGATCATGGCACAGAAAAATGACCGCCgtatc < 2:384953/90‑1 (MQ=255) gcCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCCCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc < 1:196992/90‑1 (MQ=255) gcCGAGATAAAAGGGTTTTGCTTCCCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc < 1:14340/90‑1 (MQ=255) gcCGAAAAAAAAGGGTTTTGCTTCCCGCGCGGCCTGCACCACGGGCAGCTCCCGGTTGAGATCATGGCACAGAAAAATGACCGCCgtatc < 2:379411/90‑1 (MQ=255) ccGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatcg > 2:399908/1‑90 (MQ=255) ccGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatcg > 2:439417/1‑90 (MQ=255) aaaaaaaGGGTTTTCCTCCTCGCGCGGCCTGCAACACGGGCGGCCCCCGGTTGAGATAATGGCACAGAAAAATGCCCGCCgtatcggtat < 2:380897/87‑1 (MQ=255) aTAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACGGGCAGCTCCCGGTTGAGATCATGGCACAGAAAAATGACCGCCgtatcggtatc < 2:430640/90‑1 (MQ=255) ggggTTTTGCTTCCCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCt < 1:353612/89‑1 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGAcc > 1:263406/1‑67 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGAcc < 2:263406/67‑1 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc > 1:109702/1‑75 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc < 2:109702/75‑1 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 1:210042/1‑90 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 1:380897/1‑90 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 2:43627/1‑90 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 1:128746/1‑90 (MQ=255) ttttGCTTCTCGCGCGGCCTGCAACACTGGCAGCCCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 1:232267/1‑90 (MQ=255) tttGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgcg > 1:425705/1‑90 (MQ=255) ctcGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatcggtat > 2:282699/1‑72 (MQ=255) ctcGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatcggtat < 1:282699/72‑1 (MQ=255) tcGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGa > 1:412015/1‑90 (MQ=255) ccTGGCAGCTCCGGGTTGAGATCAGGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCg < 1:17740/89‑1 (MQ=255) ccTGGCAGCCCCCGGTTGAGATCATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCg < 1:121152/89‑1 (MQ=255) cTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGa > 1:73916/1‑90 (MQ=255) cTGGCAGCTCCCGGTTGAGATCACGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGa > 1:366625/1‑90 (MQ=255) cTGGCAGCTCACGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGa > 1:139449/1‑90 (MQ=255) tGGCAGCCCCCGGTTGAGATCAGGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGaa < 2:15310/90‑1 (MQ=255) ggCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAg > 1:438255/1‑90 (MQ=255) ggCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAg > 1:275859/1‑90 (MQ=255) ggCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAg > 1:293927/1‑90 (MQ=255) cAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGat > 2:30674/1‑90 (MQ=255) cTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc > 2:465401/1‑42 (MQ=255) cTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCgtatc < 1:465401/42‑1 (MQ=255) cTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATg < 1:197367/90‑1 (MQ=255) ccgggTTGAGATCATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 1:213828/87‑1 (MQ=255) cccGTTTGAGATCATGGCACAGAAAAAGGCCCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 2:436274/90‑1 (MQ=255) cccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 1:30674/90‑1 (MQ=255) cccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc > 2:56568/1‑90 (MQ=255) cccGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 1:400281/90‑1 (MQ=255) cccGGTTGAGATCATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 2:73916/90‑1 (MQ=255) cccGGTTGAGATCATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 1:442989/90‑1 (MQ=255) cccGGTTGAGATCAGGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 2:371551/90‑1 (MQ=255) cccGGTTGAGATAATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 1:443973/90‑1 (MQ=255) ggTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc > 2:97790/1‑52 (MQ=255) ggTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgc < 1:97790/52‑1 (MQ=255) ttGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCa < 1:12809/90‑1 (MQ=255) tGAGATCATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCAt < 2:135122/90‑1 (MQ=255) tGAGATCAGGGCACAGAAAAAGGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCAt < 2:210042/90‑1 (MQ=255) gagaTCATGGCACAGCAAAATGACCGCCgtatcggtatcg > 2:30560/1‑40 (MQ=255) gagaTCATGGCACAGCAAAATGACCGCCgtatcggtatcg < 1:30560/40‑1 (MQ=255) gaTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGc > 2:269962/1‑90 (MQ=255) gaTCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGc > 1:297587/1‑90 (MQ=255) cATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc > 1:193901/1‑78 (MQ=255) cATGGCACAGAAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc < 2:193901/78‑1 (MQ=255) aTGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCtgt > 1:83943/1‑90 (MQ=255) tGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgcgc > 2:411328/1‑43 (MQ=255) tGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAgcgc < 1:411328/43‑1 (MQ=255) tGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGGCATAGCCtgtg > 2:238812/1‑90 (MQ=255) cacaGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGgct < 1:56568/90‑1 (MQ=255) caGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGaa < 2:192453/63‑1 (MQ=255) caGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGaa > 1:192453/1‑63 (MQ=255) gCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGc > 1:1587/1‑69 (MQ=255) | GATAGCTGCCCAATGCGCCGAGATAAAAGGGTTTTGCTTCTCGCGCGGCCTGCAACACTGGCAGCTCCCGGTTGAGATCATGGCACAGCAAAATGACCGCCGTATCGGTATCGATCTGAGCGCTGGCTGAGGCCGGAAAAAGATCGAAGATATGGCTGTCATAGCCTGTGGCT > NZ_CP009273/293675‑293847 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |