Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 701,009 | 0 | T | A | 45.8% | 32.1 / 34.7 | 48 | I537N (ATC→AAC) | nagE | PTS N‑acetyl glucosamine transporter subunit IIABC |
Reads supporting (aligned to +/- strand): ref base T (8/18); new base A (22/0); total (30/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.08e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG > NZ_CP009273/700921‑701098 | gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtc < 1:330080/90‑1 (MQ=255) aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCTACAGATAAAAtcg < 2:1893414/90‑1 (MQ=255) aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtcg < 1:2592921/90‑1 (MQ=255) tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAtcgt < 2:2169937/90‑1 (MQ=255) tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAacgtcg > 2:745779/1‑90 (MQ=255) gTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTa < 1:1207096/90‑1 (MQ=255) tGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTAt < 2:1867270/90‑1 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCa > 2:2803741/1‑90 (MQ=255) cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGTATCACCa > 2:1538755/1‑90 (MQ=255) tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTTGAAACCG‑CATACAAAAATGTCGGATCACCac > 1:2249144/1‑89 (MQ=255) cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa < 1:401760/90‑1 (MQ=255) ggTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCCgg > 2:474680/1‑90 (MQ=255) ttCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGGAAACGACAGATAAAATTGTCGTATTACCAGCCGCAgggg > 2:2801293/1‑89 (MQ=255) tCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGAc < 1:622075/90‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGGAAACGACAGATAAAAACGTCGTATTACCAACCGCCGGGGCa > 1:1004861/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAg < 1:2361799/79‑1 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAg > 2:2361799/1‑79 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAACCCCAGGGGCa > 2:1790609/1‑90 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGAAAAAAAAGTCGTATCACCAGCCCCCggggga > 1:1227900/1‑87 (MQ=255) ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACAACAGGTAAAAAAGGCGGATTACCCGCCGCCggggga > 1:2127695/1‑87 (MQ=255) cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGGAAACCACAAAAAAAAAAATCGGATCACCCGCCGCAGGGACaa > 1:1955342/1‑90 (MQ=255) tGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAt > 2:1008273/1‑90 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGGATCACCAGCCGCAGGGGCAATCGTg > 2:2614456/1‑90 (MQ=255) gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGGATCACCAGCCGCCGGGGCAAACggg > 2:890162/1‑88 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGa > 2:726194/1‑90 (MQ=255) aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACCACAGATAAAAACGTGGTTTCACCAACCCCAGGGGGAATCgaga > 1:1221656/1‑90 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:2564332/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:79516/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:819936/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 1:2019625/90‑1 (MQ=255) cATTCGCCAGCAAAGCGGTCGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAAAGTCGTATCACCAACCGCAGGGGGAAACGTGgaaa > 1:590388/1‑90 (MQ=255) acagcGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa < 2:2320275/86‑1 (MQ=255) ttCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAAAGTCGTATCACCAGCCGCAGGGGCAAACGTGaaaaac > 1:1016937/1‑88 (MQ=255) tCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAAACCACAGAAAAAAAAGTCGTATCACCAGCCGCAGGGGCAAAAGTGAAAAACt > 1:462663/1‑90 (MQ=255) cGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑GGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCAGGGGGAAACGTGGAAAACtt > 1:525822/1‑90 (MQ=255) gCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTc < 1:726194/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 2:2512950/90‑1 (MQ=255) ccAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa < 1:2562412/90‑1 (MQ=255) aGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa < 2:2237764/90‑1 (MQ=255) aGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa < 2:2606353/90‑1 (MQ=255) gCGGTGGGTGACGGTGTGGCGG‑GGAAAACGGCAAAAAAAAAAGTCGTATTACCAGCCGCCGGGGCAAACGTGGAAAAATTCAAAacccgg > 1:2320275/1‑87 (MQ=255) gtgggtgACGGTGTGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGAAAAACTTCAACACCAACCAc > 2:2612965/1‑90 (MQ=255) tgACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGtt < 2:65601/90‑1 (MQ=255) aCGGTGTGGCGG‑TGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAAACTTCAACACCAACCACGCGTTCt > 2:481698/1‑90 (MQ=255) gtgtGGCGG‑TGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAATCGTGAAAATCTTCAACACCAACCACGCGGTCTGcc > 2:2312593/1‑90 (MQ=255) gtgGCGG‑TGAAAACCACAAATAAAAAAGTCGTATCACCAGCCGCAGGGGGAAACGTGAAAAAATTTAAAAACCACCACCCGCTCTGCCTg > 1:1569865/1‑90 (MQ=255) ggCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa < 2:1291368/90‑1 (MQ=255) ggCGG‑TGAAAACGACAGATAAAAAAGTCGTATCACCAGCCCCAGGGGCAAACGGGAAAAAATTCAAAACCAACCCCCCGGTCTGCCTggg > 1:2216166/1‑89 (MQ=255) gAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa > 2:1891294/1‑90 (MQ=255) gAAACCGACAGATAAAATCGGCGTGTCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa > 2:644936/1‑90 (MQ=255) gAAAACGACAGAAAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAAAATTTAACACCCACCACCCGTTCTGCCTGGGAACCGa > 1:1688698/1‑90 (MQ=255) cGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAg > 1:95754/1‑90 (MQ=255) aTCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtc < 1:926190/90‑1 (MQ=255) tcgtcgTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg < 1:2624220/90‑1 (MQ=255) tcgtcgTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg < 1:1166662/90‑1 (MQ=255) | GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGG‑TGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG > NZ_CP009273/700921‑701098 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |