Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 55.6% | 32.1 / 45.5 | 37 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (1/15); new base G (20/0); total (22/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.87e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT > NZ_CP009273/822870‑823038 | gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg < 2:696028/90‑1 (MQ=255) gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg < 1:539980/76‑1 (MQ=255) gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAgg > 2:539980/1‑76 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACG‑CCGGGAg > 2:288438/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAG‑CCGGGAg > 1:433664/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAG‑CCGGGAg > 1:42558/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGAGGGGGGGTGAAAG‑CCGGagg > 1:1355971/1‑88 (MQ=255) ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAAgcgg > 2:1860307/1‑87 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACG‑CCGGGAGTCGGGtc > 1:683438/1‑90 (MQ=255) tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAG‑GCGGAAGGCGGGtcgtcgg > 2:1502261/1‑89 (MQ=255) ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACG‑CCGGGAGTGGGGtcgtcgggga > 2:558172/1‑86 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACG‑CCGGAAGTCGGTTCGTCGAGAAAc > 1:676141/1‑90 (MQ=255) cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGaa < 1:1239524/90‑1 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACC‑CCCGGAGGCGGGTCGTCGAGAAAAAGAAAG‑Tc > 2:1254695/1‑90 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑Tcc < 1:1251657/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑Tcc < 2:1515725/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑Tcc < 2:1376956/90‑1 (MQ=255) ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACG‑CCGGGAGGCGGGTTGTCGGGAAAAAGAAAG‑TCCGGTTCAt > 2:659908/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGc < 1:659908/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGc < 2:493857/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGc < 2:368665/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAG‑CCGGGAGTCGGTTCGGCGGGAAAAAGAAAG‑TCCGGTTCATGc > 1:233184/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAG‑CCGGGAGACGGGTCGTCCGGAAACAGAAAG‑TCCGGTTCATGc > 1:430838/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAG‑CCCGGAGGCGGGTTGTCGGGAAAAAAAATGTTCCGGttgtgc > 2:620456/1‑86 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACG‑CCGGGAGGCGGGTCGTCGAGAAACAGAAAG‑TCCGGTTCATGc > 1:435827/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACC‑CCGGAAGTCGGGTCGGCGAGAAAAAGAATG‑TCCGGGTCATGc > 1:368665/1‑90 (MQ=255) cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACC‑CCGGGAGGCGGGTCGTCGAGAAACAAAAAG‑TCCGGGTCATTCATc > 1:1802031/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACG‑CCGGGAGGCGGTTCGTCGGGAAAAAAAAAG‑TCCGGGTCATGCATCAg > 1:1578486/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAAc < 1:1182507/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACG‑CCGGGAGGCGGGTCGTCGAGAAAAAGAAAG‑TCCGGGTCATGCATCAGCGAAc > 1:1534743/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACG‑CCGGGAGTGGGGTGGTCGAGAAAAAGAATG‑TCCGGTTCATGCAGCAGCGGACAgg > 1:1095854/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACC‑CCGGGAGGCGGGTCGTCGAGAAAAAAAAAG‑GCCGGGTCGTGCATCAGCGaaaaag > 2:1817480/1‑86 (MQ=255) cggcggGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATT‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:1101228/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:209586/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:435827/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:29053/90‑1 (MQ=255) ggcggGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt < 2:1534743/90‑1 (MQ=255) | GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT > NZ_CP009273/822870‑823038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |