Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,435,718 | 0 | G | A | 41.9% | 59.9 / 19.5 | 43 | intergenic (‑99/+42) | uspF/ompN | stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific |
Reads supporting (aligned to +/- strand): ref base G (24/1); new base A (0/18); total (24/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.12e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT > NC_000913/1435570‑1435866 | gTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGgaagaa < 1:399428/151‑1 (MQ=255) ggAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGc > 1:574436/1‑151 (MQ=255) atatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCg > 2:22848/1‑150 (MQ=255) aTAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAAc > 2:73055/1‑151 (MQ=255) gggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGc > 1:719341/1‑151 (MQ=255) aCAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCggg > 2:565999/1‑151 (MQ=255) aaaaaTCTTTCTGTTAATAAACCCTCCTTTTTTAGTACCCGCAAAAAGTATAACCCCAGATGCCATTTTTTGTAAAAGGAGGAAGTTCCCATTTTTATTTA‑AAAACAAAGGGGAAAAAAA‑AACAAAAAAAAGCCCGCCATAACAGCgg < 2:429730/145‑1 (MQ=255) tttttttCAAAAACCCTCTTTTTTTAGTACCCGCAAAAAGTGAAACGCAAGATGACCTTTTTTGTGTAATGAGGGAGTTCAATTTTTTAATTA‑AGAAAAAAGGAGAAAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCaggagg < 1:615673/143‑1 (MQ=255) ttCTTTTAAAAAACCCTCTTTTTTTAGTACCCGAAAAAAGTGTAACGCAAGATGACATTTTTGGTAAAAGGAGGTATTACACATTTTAAATTA‑AGAACAAGGGAGAAAGAAT‑AACAAGAAAAAGCCCGCCATAACAGCGGGCaggag < 2:456096/147‑1 (MQ=255) ttCAAAAACCCTCCTTTTTTAGTATCCGAAAAAAGTGTAACGCAAGATGACTTTTTTTGTGTAATGACGGATTTCACTTTTTTAATTT‑AGAAAAAAGGAGGAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa < 1:617180/150‑1 (MQ=255) taaaaaacCCTCTTTTTTTATTACCCGAAAAAAGTGTAACGCAAAATGACTTTTTTTGTGAAAGGACGGATTCCACATTTTAAATTA‑AGAGAAAGGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa < 2:239504/145‑1 (MQ=255) tCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc > 1:497294/1‑151 (MQ=255) aaacCCTTCTGTTTTAGTATCCGAATAATGTGTAACGCAAGATGACTTTTTTGGTGTAATGAGGGATTTCACTTTTTTAATTT‑AAATCAAGGGAGAAAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt < 2:784747/148‑1 (MQ=255) aaacCCTCTTGTTTGAGTAACCACAGAAAGTTAAACCACAGATGACTTTTTTTGTGAAATGAGGGAGGCCACATTTTAAATTA‑AGATCAAAGG‑GAAAAAATGAAACAAAAAAAGTCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa < 2:577358/150‑1 (MQ=255) aaacCCTCCTTTTTTATTACCCTAAAAAAGTAAAACGCAAAATGACTTTTTTGGTGAAATGACGGAGTTCAC‑TTTTTAATTTAAGAAAAAAGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTa < 2:401517/143‑1 (MQ=255) aaaCCTCGTGTTTGATTACCCGAAAAAATTGAAACGCCAGAGGACATTTTTTGGGAAAGGACGGAGTTCGCATTTTAAATTA‑AGAAAAAGGGAGAAAAAAT‑AAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa < 2:643589/151‑1 (MQ=255) aaCCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATa > 1:642576/1‑150 (MQ=255) cctcttGTTTGAGTACCGGAATAAAGTGTAACGCCAGATGACATTTTTTGTGAAAGGCGGGAGTTCACATTTTTAATTT‑AGATAAAGGGAGAAAAAAA‑AAGAAGAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaaa < 1:356172/150‑1 (MQ=255) ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGGGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc > 2:467563/1‑151 (MQ=255) ctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAcc > 2:107763/1‑151 (MQ=255) tcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc > 1:402686/1‑149 (MQ=255) gTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc > 1:168376/1‑151 (MQ=255) aGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa > 1:398442/1‑149 (MQ=255) tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg > 1:855855/1‑151 (MQ=255) ccGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc > 2:96323/1‑150 (MQ=255) aaaaagggaaccGCCAAATCACTTTTTTTGTAAAAAGAGGGATTTCATATTTTAAATAA‑AAATGAAGGAAAAAAAAAA‑AAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACaa < 2:645183/137‑1 (MQ=255) aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat > 1:130944/1‑151 (MQ=255) aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat > 1:110600/1‑151 (MQ=255) accGCCAGATGACAATTTTGTTATAATGAGGGAGTTCACTTTTTTATTTT‑AGAACAAGGGAGAAAAATT‑AAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCgg < 2:771604/146‑1 (MQ=255) agagGACATTTTTTGTGTAAGGACGGAGTTCATTTTTTTAATTA‑AGAAAAAAGGAGAAAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATg < 2:551526/150‑1 (MQ=255) aTGACTCTTTTTGTGAAAGGAGGGGTGTCAAATTTTAAATTA‑AAAAAAAAGGAGAAAAAAA‑AAGAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag < 2:123515/145‑1 (MQ=255) tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCAt > 1:484924/1‑151 (MQ=255) actttttttGTGAAATGAGGGATTTCAATTTTTTAATTA‑AGAAAAAGGGAGAAAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGc < 1:784722/142‑1 (MQ=255) tttttttGTATAGGGGCGGAGTACACATTTTAAATTT‑AAAAAAAAGAAGAAAAAAA‑AAAAAAAAAAGTCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg < 1:421148/148‑1 (MQ=255) ttggtatAATGAGGGATTCCAATTTTTTAATAT‑AAATAAAAGGAGAAAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGc < 2:651500/144‑1 (MQ=255) gtgtAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt > 2:512508/1‑149 (MQ=255) atttCACATTTTAAATTG‑AGAAAAAAGGAAGAAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATc < 2:336481/148‑1 (MQ=255) gTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCtt > 1:165839/1‑151 (MQ=255) ttCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTc > 1:648091/1‑151 (MQ=255) caTTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATcc > 2:390252/1‑151 (MQ=255) aaTTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt > 2:167803/1‑151 (MQ=255) ttt‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg > 1:49292/1‑151 (MQ=255) tCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt > 2:640786/1‑150 (MQ=255) ggagaaAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt < 2:227005/150‑1 (MQ=255) gaagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg > 2:378158/1‑151 (MQ=255) aagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtata > 2:94060/1‑143 (MQ=255) aagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGt > 1:817862/1‑151 (MQ=255) | GTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTT‑AGATCAAAGGAGGAAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGT > NC_000913/1435570‑1435866 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |