Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,349,145 | 0 | A | C | 45.2% | 59.6 / 13.3 | 42 | F110V (TTC→GTC) | inaA | acid‑inducible Kdo/WaaP family putative kinase |
Reads supporting (aligned to +/- strand): ref base A (20/3); new base C (0/19); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.49e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG > NC_000913/2349002‑2349287 | gCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGcc > 2:38570/1‑151 (MQ=255) aTGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGc > 2:855740/1‑151 (MQ=255) gACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGc > 1:336215/1‑151 (MQ=255) cGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTc > 2:606217/1‑151 (MQ=255) ttAATGTTTTGATTTTTCTAAACCGCC‑GCCCCCCTGGTTTTAAACTGGGTTTCCCGTCCTTCTTAAAAAAAGGGCGT‑TCCTGCATGCTGGCCATCCCAGACAGCATTGCTAAAGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCg < 1:850393/149‑1 (MQ=255) ttAATGCTATGCTTTTTCTTAACCCCCAGCGCCACTGTTTTCAACTTGGTTTG‑CCGTTCTTCGTCAGAATAAGGCGA‑TCCTGCATGCTGGGCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCaa < 2:37207/143‑1 (MQ=255) aaTGTTATGCTTTTTCTAAACCGCCGGCCCCCCTGTTTTAAACGTGCCTTC‑CCGTACTTCGTCAAAATAGGGCGA‑TCCTGCAGGCGGGGAATCCCAGTCAGCAATGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGcc < 2:542968/149‑1 (MQ=255) ctatttaTTTTCTAAACCCCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTTCTCCGTAAAAATAAGGCGA‑TCCTGCATGCGGGGCATCCCAGACAGAATTGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt < 2:789035/142‑1 (MQ=255) tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc > 1:476362/1‑151 (MQ=255) aacccccAGCGCCCTTCTTTTAAACATGGCTTC‑CCTTACTTCG‑CAAAAAAAGGGAATTCCTGCTTGCGGGGCAACCCAGCCAACAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCa < 1:19323/144‑1 (MQ=255) cGCCAGCCCCCCTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAAAATAAGGCGA‑TACTCCTTGCGGGGCATCCCAGTCACCATTTCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc < 1:850820/145‑1 (MQ=255) gCCATCGCCACTGCTTTCAACATGGCTTG‑CCGGACTTCGTCAGAATAAGGCGA‑TACTGCACGCTGGGCATCCCAGTTAGCAATCCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc > 1:119167/1‑147 (MQ=255) ccggccccccTGCTTCAAACAGGGTTTC‑CCTTACTTCTAAAAAAAAAGGCGA‑TACTGCATGCTGGGCATCCCAGCCAGAAATGTTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc < 2:272482/142‑1 (MQ=255) cAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc > 1:806311/1‑150 (MQ=255) cgcCACTGTTTTCAACATGGCTTC‑CCGTACTTCGTAAAAATAGGGGGA‑TACTGCATGCGGGCCATCCCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa < 2:582460/151‑1 (MQ=255) cACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa > 2:921731/1‑150 (MQ=255) ttCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt > 2:766969/1‑151 (MQ=255) aaCATGGCTTG‑CCGTACTTCGTCAAAATAAGGCGA‑TACTGCATGCGGGGCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt < 2:702064/148‑1 (MQ=255) aCATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCgg > 1:629203/1‑151 (MQ=255) cATGGTTCG‑CCGTACTTCCACAAAAGAGGGGTA‑TCTCGCAGGCGGGTCATCCCAGAAAGAAATCCAGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt < 2:227622/146‑1 (MQ=255) ttG‑CCTTTCTTCGTAAAAATAAGGCGA‑TCCTGCATGCGGGGCATACCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACg < 2:600723/148‑1 (MQ=255) ttG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATg > 2:76639/1‑151 (MQ=255) cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac > 2:726949/1‑151 (MQ=255) cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac > 1:280816/1‑151 (MQ=255) cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac > 1:357948/1‑151 (MQ=255) gTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACAcc > 1:234327/1‑151 (MQ=255) ttcttaaaaaTAAGGCTA‑TACTGCGTGCTGGGCATACCAGCCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCgg < 2:337310/141‑1 (MQ=255) cAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTc > 2:131718/1‑151 (MQ=255) aaaaaGGCTA‑TCCGGCAGCCTGGCCATACCAGTCAGCAGTACACATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCg < 1:166347/142‑1 (MQ=255) aaTAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc > 2:249504/1‑149 (MQ=255) aTAAGGCTA‑TACTGCTGGCTGGGCTTCCCAGTCAGAATTGTGGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCa < 2:254120/149‑1 (MQ=255) ggCGA‑TGCTGCGGGCGGGGCATCCCAGACAGCAATGCGGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTc < 1:411671/149‑1 (MQ=255) ggCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt > 1:2433/1‑151 (MQ=255) gccGA‑TCCTGCATGCGGGGCAACCCAGCCAGAATTGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt < 2:292267/149‑1 (MQ=255) a‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAt > 2:752251/1‑151 (MQ=255) tGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCt > 1:176300/1‑151 (MQ=255) gggCATCCCAGCCAGCATGGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGc < 2:339198/145‑1 (MQ=255) gggCATCCCAGCCAGCAATGCGGAACACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTc < 1:594882/151‑1 (MQ=255) cATCCCAGCCAGCAATGTTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCt < 1:809867/147‑1 (MQ=255) aGTCAGCAATGCTAATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAAt < 1:542886/151‑1 (MQ=255) gacAGCATTGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATc < 1:580657/149‑1 (MQ=255) atgcACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCgg < 1:825578/147‑1 (MQ=255) | GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG > NC_000913/2349002‑2349287 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |