Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 58.3% | 38.7 / 60.4 | 48 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (2/18); new base G (28/0); total (30/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.60e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT > NZ_CP009273/822870‑823038 | gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg < 2:586622/90‑1 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggg > 2:1640364/1‑90 (MQ=255) cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAggg > 2:886725/1‑90 (MQ=255) tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCa < 1:1241905/90‑1 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAg > 1:613807/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCggggg > 2:1737594/1‑87 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCCGAAg > 1:678804/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAg > 2:356574/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGAAg > 2:492934/1‑90 (MQ=255) ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCGGGAGTCgg > 2:1274204/1‑90 (MQ=255) gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTGGGt > 1:1667789/1‑90 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGACGGgtcgtc > 1:569387/1‑90 (MQ=255) tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAGGCGGGTCGGCga > 1:832007/1‑90 (MQ=255) ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTGAACGCCCGAAGGCGGTTCGTCCAGaa > 2:1304770/1‑90 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAc > 2:1122567/1‑90 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCGGAAGGCCGGTCGTCGGGaaaaaa > 1:913031/1‑87 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTGGGGTTGTCCGGaaaaaa > 1:1197682/1‑87 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGaaaaaa > 1:228811/1‑87 (MQ=255) cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGt < 2:881484/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:54053/90‑1 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:356574/90‑1 (MQ=255) tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCg < 2:391613/90‑1 (MQ=255) cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATg > 2:694386/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGCAGGCGGGTCGTCCGGAAAAAGAAAGGCCGGGTCCTgg > 1:1368995/1‑89 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAACCGGGAAACGGGTCCTCGAGAAAAAGAAAGTCCCGGTCAttc > 2:1760915/1‑88 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTTGTCGAGAAACAGAAAGTCCGGTTCATGc > 2:729318/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTTATGc > 2:1106695/1‑90 (MQ=255) gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgt < 1:287838/62‑1 (MQ=255) gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgt > 2:287838/1‑62 (MQ=255) cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCGACGCCCGAAGTCGGGTTGTGGGGAAAAAAAAAATCCGGGTGATGCGTc > 1:1875930/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 1:1369464/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:1209785/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCGGAAGGGGGGTCGGCGAGGAAAAGAAAGTTCGGGTTATGCATCAg > 1:217585/1‑90 (MQ=255) aTGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGaa < 2:1584178/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:345718/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:197496/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:174260/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1702031/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:832007/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1329612/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGTCGGGTTGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc > 2:233567/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAAc > 2:542408/1‑90 (MQ=255) gtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACa < 2:573641/90‑1 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCGGAAGGGGGTTCGTGGAGAAAAAGAATGGCCGGGTCGTGCATAAGCGGACAGacc > 1:962947/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGGCGAGAAAAAGAATGTCCGGTTCATGCATCAGCGAACAAgcc > 2:718192/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGTCGGGAAACAGAAAGTCCGGGTTCTTCATCAGCGAAAAgggc > 2:1635507/1‑88 (MQ=255) aGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGTTCATGCATCAGCGAACAGgcc > 2:708116/1‑90 (MQ=255) cAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTTGTCGGGAAAAAGAAAGTCCGGGTCATTCATCAACGAACAAgccagc > 1:358694/1‑90 (MQ=255) ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt < 1:1762307/90‑1 (MQ=255) | GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCT > 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 < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |