Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,963 | 0 | T | G | 55.6% | 9.5 / 42.6 | 36 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (1/15); new base G (20/0); total (21/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 = 2.31e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822886‑823045 | 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 < 2:1515725/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑Tcc < 2:1376956/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑Tcc < 1:1251657/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:368665/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGc < 2:493857/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) gAGGGGGTCACCC‑CCGGAAGCCGGTCCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCCAACAGGCCAGCGCCCGCCGCTGTTTaaa < 2:1033632/90‑1 (MQ=255) gAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:1443741/90‑1 (MQ=255) | GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACG‑CCGGAAGTCGGTTCGTCGAGAAACAGAATG‑TCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822886‑823045 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |