Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,963 | 0 | T | G | 57.5% | 21.4 / 62.8 | 47 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/17); new base G (27/0); total (30/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.16e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT > NZ_CP009273/822876‑823052 | 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 < 2:1209785/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 1:1369464/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCGGAAGGGGGGTCGGCGAGGAAAAGAAAGTTCGGGTTATGCATCAg > 1:217585/1‑90 (MQ=255) aTGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGaa < 2:1584178/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1329612/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:832007/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1702031/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:174260/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:197496/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:345718/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) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 2:695406/1‑90 (MQ=255) tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt > 2:1042480/1‑90 (MQ=255) | TCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT > NZ_CP009273/822876‑823052 |
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: ‑ |