Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,791,931 | 0 | G | A | 29.6% | 36.3 / 11.9 | 27 | H637Y (CAC→TAC) | pvdI | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base G (10/9); new base A (8/0); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.58e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCATGAATGCCATTCGCTGCT > NC_002947/4791785‑4792074 | cATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGCTCACATAAGCCAGGTTGCCGGAGTGa > 5:74989/1‑148 (MQ=255) cATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCACGTTCCCGTAGTg > 2:84637/1‑147 (MQ=255) cAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGCCCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGcccc > 1:13520/1‑148 (MQ=255) aaGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGCCCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGAGCACCCCTGCCTCg > 3:27690/1‑148 (MQ=255) aaTCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACAGGTTAACATATTAAAAATCAACGGGCACCCCTGCCTCGCCACCt > 2:81961/1‑148 (MQ=255) tCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACTTAAAAAACAAACACGGGCACCACTGCCTCGCCACCTGCcgcg > 4:124460/1‑147 (MQ=255) aaGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAgcg > 7:156491/1‑148 (MQ=255) aGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAAACCCACGGACACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGa > 4:72434/1‑148 (MQ=255) gcCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAACGAAACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATgg > 2:89975/1‑148 (MQ=255) acacCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGGTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAgg < 5:95660/148‑1 (MQ=255) cgcgACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAAAAAAACGGGCAACCCTGCCTCGCCACCTGCCGCGTTCAGTGCCAGGCACTGTACGCCCCCGATTGCCGGGCCAGGCTCg > 2:99425/1‑148 (MQ=255) aCCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAAAAAAACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCt > 4:83059/1‑148 (MQ=255) ccTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCgg > 3:15954/1‑130 (MQ=255) cTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCggg < 4:15954/130‑1 (MQ=255) ggCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCt < 1:102512/146‑1 (MQ=255) gCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGa < 1:97716/147‑1 (MQ=255) gCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAAAAACACAGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCCCTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGa > 4:3726/1‑147 (MQ=255) gCCCGAGGTGTAGATCACATAAGACACCTTAACATAAAAAAAAAAAACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATc > 2:156705/1‑148 (MQ=255) cGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGTCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCtgt < 3:113101/148‑1 (MQ=255) gtgtAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGt < 4:168454/148‑1 (MQ=255) tgtAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGCTGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTc > 5:112048/1‑148 (MQ=255) aTCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACgca > 5:108242/1‑148 (MQ=255) caTAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCAt < 6:96162/148‑1 (MQ=255) aTAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCAt > 6:127698/1‑147 (MQ=255) aTAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCAt > 4:135710/1‑147 (MQ=255) gCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCATGAAt < 8:156491/147‑1 (MQ=255) gCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCCGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCATGAATGCCATTCgc > 8:122198/1‑148 (MQ=255) gAGTGAACCGCCACGAGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCATGAATGCCATTCgctgct < 4:27690/148‑1 (MQ=255) | CATGGCGAAGTGTTCCAGATCGCCAAGGGCAATCCGCTCGGACAAGCCTTCAAGGTAATTGGCCAGCGCCCCATGAGACACCGCGACCCCCTTGGGCATACCGGTGGTGCCCGAGGTGTAGATCACATAAGCCAGGTTGCCGGAGTGAACCGCCACGGGCACCCCTGCCTCGCCACCTGCCGCGTTCAGCGCCAGGCACTGTACGCCCCCGATGGCCGGGCCAGGCTCGTCCTCGGCCCTGAGCAAGACCTTGATCTGTGCGTCACGCAGCATGAATGCCATTCGCTGCT > NC_002947/4791785‑4792074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |