Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,957 | 0 | A | G | 54.5% | 28.7 / 48.0 | 44 | L494P (CTC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (6/14); new base G (24/0); total (30/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.37e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822878‑823045 | cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAAc < 1:441531/90‑1 (MQ=255) gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCggg > 2:1125990/1‑89 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAg > 2:1249155/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAg > 2:608281/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAg > 1:746796/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAg > 2:98221/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAg > 2:620505/1‑90 (MQ=255) ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAgg > 2:895497/1‑89 (MQ=255) ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGACACGCCGGAAGt > 2:1457921/1‑90 (MQ=255) gtggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTc > 2:76419/1‑90 (MQ=255) tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGtt > 2:283380/1‑90 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCgggtt > 2:113012/1‑87 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGTCGGGtc > 2:471882/1‑90 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCGGgtcgtc > 2:1394570/1‑90 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAc > 2:228978/1‑90 (MQ=255) gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGAAGGCGGGTCGTTGAGAAACaaaaa > 1:601523/1‑89 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTc > 2:33214/1‑90 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:1164295/90‑1 (MQ=255) ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt < 2:36663/90‑1 (MQ=255) ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt < 1:1190414/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGCTTCATGc < 1:1571081/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGc > 2:1164295/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGAAGTCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGc > 2:1155912/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTTGTCGGGAAAAAGAAAGTTCGGGTTATGc > 1:36663/1‑90 (MQ=255) gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAAAAAGGCCGGGTCATGCAt > 2:1385084/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:758476/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 1:498072/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:1140547/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAAAAAAATGTCCGGGTCATGCATCAg > 2:1190414/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGGGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCac > 1:1516367/1‑89 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGCGGGGGGGGGAAAGCCCGGGGTGGGGTGGGTGAGAAAAAGAAAGGTCGGGTGATGCATCACCg > 1:1563221/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGGCGGGAAACAAAAAATCCGGTTTATGCAACCGCGGAc > 2:577908/1‑90 (MQ=255) gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACa < 1:176074/90‑1 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACaag > 2:879199/1‑88 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACCCCCGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAgg > 2:1112031/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGGCGAGAAAAAGAAAGTCCGGTTCATGCATCAGCGAAAAgg > 2:718102/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAAAAAAAAAAAGTGCGGGGCATGCAGCAGCGGACAGgcc > 1:734318/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGGGGGGTCGTCGAGAAAAAAAAAGTGCGGGTCATGCATCCCCGAACAgggc > 1:1205825/1‑88 (MQ=255) aaTTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTTGTCGGGAAAAAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcc > 2:346895/1‑90 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:419860/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 1:476325/90‑1 (MQ=255) gcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTg < 2:1516367/90‑1 (MQ=255) gcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTg < 1:1347840/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 1:620505/90‑1 (MQ=255) | CATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822878‑823045 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |