Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,161 | 0 | C | T | 60.6% | ‑3.3 / 32.1 | 33 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (1/12); new base T (20/0); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.66e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCG > NZ_CP009273/4594076‑4594239 | aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt > 1:337489/1‑90 (MQ=255) cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 1:156614/90‑1 (MQ=255) cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc < 1:712961/90‑1 (MQ=255) cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg < 2:417029/90‑1 (MQ=255) cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg < 1:675892/90‑1 (MQ=255) ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCC‑a > 2:97896/1‑90 (MQ=255) aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTg < 1:315201/90‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑ATTGGAAGGGAgggcg > 1:393303/1‑89 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑ATGAGAAGGGAGGGCa > 2:660347/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCC‑CTGAGAAGGGAGGGCa > 2:233193/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCCTTTTTTTCC‑ATGAGAAGGGAGGGCa > 1:475241/1‑90 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCaga < 1:97896/90‑1 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCaga < 2:265665/90‑1 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCaga < 1:301246/90‑1 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCagaga < 2:287861/90‑1 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑CTGAGAAGGGAGCGCagggag > 2:82555/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCC‑CTGGGAAGAGGGCGCagagag > 1:707310/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCC‑CTGAGGAAGGAGGGGagagag > 2:477104/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCAATGGGAA‑TGAGAGCCGAGAgg > 2:531305/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑CTGAGAAGGGAGGGCCGAGAgg > 2:636760/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCC‑CTGGGAAGGGAGCGCAGAGAgg > 1:717177/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCC‑CTGAGGAGTGAGCGCAGCGAgg > 2:786783/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑CTGGGAAGGGGGGGCAGgggggg > 2:299138/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑ATGAGAAGGGAGCGCAGAGAggg > 1:133094/1‑90 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑ATGAGAAGGGAGGGCCGAGAgggg > 2:526946/1‑89 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCC‑CTGAGAAGGGAGGGCAGGGGGGGGGTATCtt > 2:215569/1‑90 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑CTGAGAAGGGGGCGCAGAGGGGGGTTATCTTGc > 1:94747/1‑90 (MQ=255) ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc < 1:721196/90‑1 (MQ=255) cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt < 1:67686/90‑1 (MQ=255) tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCC‑CTGGGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCg > 1:47115/1‑90 (MQ=255) tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCC‑ATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCg > 1:773832/1‑90 (MQ=255) tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt < 2:740453/90‑1 (MQ=255) gcggcgGTCATTTTTTTTTTTTTTGCC‑CTGAGAAGTGAGGGCGGAGAGGGGTTATTTTTCCTGACTGAACTGGGGATTAAGTTTTCCCGt > 2:766568/1‑90 (MQ=255) ggTCATTTTTTTTTTTTTTTCC‑ATGAGAAGTGAGCGCCGAGGGGGGTTTTTTTTCCTGTCTTAACTGGCGATTAAGTTTTCGCGGACCCg > 2:423766/1‑90 (MQ=255) ggTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 2:90815/90‑1 (MQ=255) ggTCATTTTTTCCTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 2:475241/90‑1 (MQ=255) | AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCC‑ATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCG > NZ_CP009273/4594076‑4594239 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |