Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 599,576 | 0 | T | G | 47.4% | 82.7 / 15.2 | 57 | V288G (GTG→GGG) | cusA | copper/silver efflux system, membrane component |
Reads supporting (aligned to +/- strand): ref base T (2/28); new base G (27/0); total (29/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.09e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCA > NC_000913/599427‑599716 | aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg < 1:746569/151‑1 (MQ=255) gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATcc > 1:676078/1‑148 (MQ=255) gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc < 2:964305/151‑1 (MQ=255) gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCt > 1:373332/1‑149 (MQ=255) ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑a < 1:301628/149‑1 (MQ=255) tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGGGCCTGCgcgcc > 2:915890/1‑147 (MQ=255) caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCaaa < 1:87310/151‑1 (MQ=255) tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGGGGGGGGGTCCCGCGCACCCGGCaaa > 2:648117/1‑150 (MQ=255) tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCCGctggcccgggaaaa > 2:41864/1‑134 (MQ=255) tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCG‑GTTCGGGAAAAAAGccc > 1:206598/1‑150 (MQ=255) aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCG‑ATCCGGCAAAAACGccc > 1:301275/1‑148 (MQ=255) gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGAt < 2:691266/151‑1 (MQ=255) gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATcg < 2:767268/151‑1 (MQ=255) aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCCTCG‑CAACGGCAAAAAAACCCgaggagagggcg > 2:933261/1‑142 (MQ=255) aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 2:46478/151‑1 (MQ=255) ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc < 2:532447/150‑1 (MQ=255) gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt < 1:695011/151‑1 (MQ=255) cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGATCgccgcc > 1:372878/1‑148 (MQ=255) gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGccgccgcgg > 1:78309/1‑148 (MQ=255) gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa < 2:205714/149‑1 (MQ=255) gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa < 1:705391/151‑1 (MQ=255) ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa < 1:167308/150‑1 (MQ=255) tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCG‑AACCCGCAAAAAAGACCGAGAAGGGGAAGCCGCCGAGAAGGACaaa > 1:689282/1‑150 (MQ=255) tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg > 1:64160/1‑151 (MQ=255) cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGCCTGCG‑GTCCGGGAAAAAAGACCGCGGAGAGATCGCCGCCGCGAAAGaaaa > 1:214027/1‑146 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt < 1:752122/147‑1 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt < 2:249587/147‑1 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTgg > 1:774866/1‑149 (MQ=255) aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAAAATGCG‑ATCCGGCCAAAAAACACGACGCGCGATCACAGCCGGGGAGGACaaaagcggaacg > 2:903140/1‑140 (MQ=255) gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa < 2:1001433/150‑1 (MQ=255) gCGAAGCTCCAGAGTGGCCCGGAGATGCGCCGGGACATTGCCGAACTAAACGGCGAAGGCGAAGTAGCGGGCGGGGTGGTGATCCAGCG‑ATCCGGCAAAAACGCCCGAGTAGTGAGCGCCGCCGTGAAGGACAAACTGGAAACGGTGaaaa < 2:340597/151‑1 (MQ=255) aGGTCCAGATTGGCGCAGAGATGCGCCGTGGCATTGGCAAACTAAACGGCGATGGCGAAATGGCCGGGGGCGTGGTGATCCTGCG‑ATCCGGCAAAAACGCTCGAGAAGTGATGGGCGCCGTGAAGTACAAACTGGAAACGCTGAAAagtag < 2:875516/151‑1 (MQ=255) aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag < 1:831038/151‑1 (MQ=255) aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCG‑GTCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAACGCTGAAAAGTAGTCTGcc > 1:204740/1‑151 (MQ=255) cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc < 1:497314/151‑1 (MQ=255) cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAgg > 1:330811/1‑150 (MQ=255) cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGCGGCGGGCGGGGGGGTGATCCCGCC‑ATACGGCAAAAACGCCCGCGGAGTGATCGCCGCCGTGACGGACAAACAGGAAACCCTCAAAAGAAACCTGCCGGAAg > 2:555442/1‑149 (MQ=255) ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCG‑GTCCCGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGCAGTCTGCCCGAAgg > 1:845352/1‑149 (MQ=255) ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGGGTTCCTGCG‑CTACCGCAAAAACGGCCGCGGTGGGATCGGCGCCGCGGAGGAAAAACAGGAAAAGGTTAAACAAAGTCTGGCGGAAGGCggg > 2:656290/1‑149 (MQ=255) gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCG‑CTCCCGCCAAAACCCCCGAGAAGGGATCGCCCCCCGGAAGGACCAACTGGAAACGCCGAAAAGGAGTCTGCCGGGAGGCgggg > 2:434083/1‑151 (MQ=255) gaTGCGCCGGGGCATTGCCGAGCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga < 2:30376/150‑1 (MQ=255) gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt < 1:150276/151‑1 (MQ=255) tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa < 1:939691/150‑1 (MQ=255) cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa > 1:521236/1‑150 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt > 1:224319/1‑150 (MQ=255) tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCG‑ATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCCGGAGGGGTGGGGATAGTTACAACATACGATCg > 1:197886/1‑148 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca < 2:936229/149‑1 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCCGGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGGGGGGATAGTTACAACATACGATCGca > 1:992914/1‑149 (MQ=255) aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca < 2:475926/148‑1 (MQ=255) aCGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCCCG‑ATCCGGGAAAAAACCCCGAGAAGAGATCCCCGCCCGGGAGGACAAACTGGAAACGCCGAGAAGTAGGCTGCCGGAGAGCGGGGGGGTAGTTAAAACAAACGATCGca > 1:681810/1‑147 (MQ=255) cGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca < 1:685786/150‑1 (MQ=255) gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc < 2:971381/150‑1 (MQ=255) aggcgaagTGGCGGGCGGGGGGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGGACGCCGCCGTGAAGGGCAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTAAAACAAAAGATCGCAGCCCGCTCAtt > 2:81046/1‑151 (MQ=255) aagTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg < 2:685757/151‑1 (MQ=255) aagTGGCGGGCGGGGGGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAACGCAGCCAGCTCCTTGAcc > 2:402821/1‑150 (MQ=255) gTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg < 2:490162/151‑1 (MQ=255) ggcgggcgGGGTGGTGATCCTGCG‑ATACGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGCGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc > 2:807873/1‑151 (MQ=255) gcgggcggGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa < 1:430852/151‑1 (MQ=255) | AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCA > NC_000913/599427‑599716 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |