Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 321,484 | 0 | A | T | 21.2% | 65.6 / 14.3 | 33 | intergenic (+403/‑124) | rclR/ykgE | reactive chlorine species (RCS)‑specific activator of the rcl genes/cysteine‑rich LutA family protein; putative electron transport chain YkgEFG component |
Reads supporting (aligned to +/- strand): ref base A (12/14); new base T (5/2); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.98e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATAGGTATAAGGGATGTAGCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTT > NC_000913/321340‑321625 | aaTAGGTATAAGGGATGTAGCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATAc > 1:124938/1‑148 (MQ=255) aTGTAGCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTa < 2:155447/149‑1 (MQ=255) tttAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGc < 2:124938/148‑1 (MQ=255) ttAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGcc > 6:69273/1‑148 (MQ=255) aaTGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAAt < 5:69273/148‑1 (MQ=255) gATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTaaa < 2:49304/148‑1 (MQ=255) aTTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTaaaaa < 8:22357/149‑1 (MQ=255) ttATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTATGACAAATGataat < 6:50528/119‑5 (MQ=255) ttATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTATGACAAATGataat > 5:50528/1‑115 (MQ=255) ccAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTAt < 6:170056/120‑1 (MQ=255) ccAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTAt > 5:170056/1‑120 (MQ=255) aaaaTGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCAtatt > 1:167593/1‑149 (MQ=255) aaTGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCAtattta > 3:122076/1‑149 (MQ=255) gtcgATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATtgatga < 4:122076/149‑1 (MQ=255) cATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCg > 2:140642/1‑149 (MQ=255) cATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTATGACAAATGata > 8:31774/1‑87 (MQ=255) cATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTATGACAAATGata < 7:31774/89‑3 (MQ=255) tgAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTgggggg < 2:167593/148‑1 (MQ=255) aGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCgg < 8:39571/139‑1 (MQ=255) aGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCgg > 7:39571/1‑139 (MQ=255) gaaagaaaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCggg > 6:129810/1‑148 (MQ=255) aaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAAt > 8:147568/1‑149 (MQ=255) aTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATg > 6:93612/1‑149 (MQ=255) aTTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCgg < 5:151752/128‑1 (MQ=255) aTTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCgg > 6:151752/1‑128 (MQ=255) aTTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa > 7:70062/1‑149 (MQ=255) cATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGa < 7:147568/149‑1 (MQ=255) tttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGtt < 8:70062/149‑1 (MQ=255) gATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa < 5:129810/149‑1 (MQ=255) cctCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGCACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAAttt > 3:101610/3‑149 (MQ=255) cctCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAAtt > 8:129925/3‑148 (MQ=255) ctCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTc > 8:135027/2‑149 (MQ=255) aTTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCttt < 5:93612/149‑1 (MQ=255) aTACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa < 6:128804/128‑1 (MQ=255) aTACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa > 5:128804/1‑128 (MQ=255) | AATAGGTATAAGGGATGTAGCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTT > NC_000913/321340‑321625 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |