Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,528,480 | 0 | A | T | 20.0% | 75.4 / 14.6 | 35 | Y671F (TAC→TTC) | yrfF | putative RcsCDB‑response attenuator, inner membrane protein |
Reads supporting (aligned to +/- strand): ref base A (14/14); new base T (3/4); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTGCG > NC_000913/3528332‑3528614 | tttCACGCCGAGGGGATCGTCACCAAAATCCTCACTGCCGCCAATGGTACGCAGCATATTTGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTa < 6:13299/149‑1 (MQ=255) cGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCg > 6:72550/1‑149 (MQ=255) gAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGCTACTCAGCATATTGTCCTTCATCCGATCGCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGCTGCTGGGTAGCGCCATTGTAAATGGCgc < 7:31039/149‑2 (MQ=255) aGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTa < 3:92929/149‑1 (MQ=255) aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGcgtc > 6:108481/1‑149 (MQ=255) aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGcgtc < 5:108481/149‑1 (MQ=255) aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGc < 6:126606/146‑1 (MQ=255) aaaaaTCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGc > 5:126606/1‑146 (MQ=255) cACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCa < 5:72550/147‑1 (MQ=255) cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAAt > 1:146095/1‑90 (MQ=255) cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAAt < 2:146095/90‑1 (MQ=255) cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTa > 6:41302/1‑113 (MQ=255) cTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTa < 5:41302/113‑1 (MQ=255) tCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCtat > 1:39989/1‑148 (MQ=255) tCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGa > 8:14134/1‑148 (MQ=255) ttCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGTGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGc < 3:32281/145‑1 (MQ=255) ttCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGAGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGc > 4:32281/1‑145 (MQ=255) gcgcTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc < 8:120974/146‑1 (MQ=255) gcgcTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc > 7:120974/1‑146 (MQ=255) aGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTc < 2:39989/149‑1 (MQ=255) cacaTTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGcc < 2:4576/119‑1 (MQ=255) cacaTTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGcc > 1:4576/1‑119 (MQ=255) aCTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGaaa > 8:20670/1‑137 (MQ=255) aCTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGaaa < 7:20670/137‑1 (MQ=255) cTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAAt < 7:14134/148‑1 (MQ=255) tGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc > 8:50437/1‑108 (MQ=255) tGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAAc < 7:50437/108‑1 (MQ=255) gggTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTgcg > 6:72225/1‑149 (MQ=255) ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa > 4:60031/6‑135 (MQ=255) ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa < 3:60031/130‑1 (MQ=255) ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa < 3:62122/130‑1 (MQ=255) ctgtacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa > 4:62122/6‑135 (MQ=255) tacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa < 1:51461/130‑1 (MQ=255) tacgcCATTGTAAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGaa > 2:51461/3‑132 (MQ=255) cgcCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAgg < 7:118801/71‑1 (MQ=255) cgcCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAgg > 8:118801/1‑71 (MQ=255) | TTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGCGTCATCGCACTCGCATGATGGAGATTCAGGCCTATTATGAAAGCTGCCTGAACCCGCAACTGATCACCCCTTCAGAAAGCCTTATCGAATAACACGTTTGCG > NC_000913/3528332‑3528614 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |