Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,705,661 | 0 | T | C | 51.7% | 62.4 / 86.9 | 58 | I44V (ATC→GTC) | dppB | dipeptide/heme ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base T (22/6); new base C (0/30); total (22/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.70e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-11 |
AACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTA > NC_000913/3705516‑3705805 | aaCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAc > 2:532104/1‑151 (MQ=255) aCTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACg > 1:297580/1‑151 (MQ=255) tctTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGtt > 2:140662/1‑151 (MQ=255) tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTc < 1:444830/149‑1 (MQ=255) cGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCAAAATGTAATGGGGATACGGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCAGGTTcgcccgcc < 2:365925/148‑1 (MQ=255) ggATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGat > 1:342838/1‑151 (MQ=255) ggATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGat > 2:129649/1‑151 (MQ=255) cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGTCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc > 2:933691/1‑150 (MQ=255) ttttCATTAAAATGCCTAGCCCGCCATGCATACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCTAGGCCGAGTTCACCCCGCCGCTGCGCGTGCCGCTCGGGGGAGCCCCCACGTTCGCCCGCCATGATCATCAc < 2:158276/148‑1 (MQ=255) ttttataaaaaTGCCAAGACCCCCATGCATACCGCCCAAAATGGAAGGGAGATACTGCTGCCACAGCGGTTTATCAAACCCGAGTTCACCCAGCAGCTGCGCTTGCCGCTCGGGGGAGCTCCCACGTTCGCCCGCCATGAtcatca < 2:730989/138‑1 (MQ=255) cATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCg > 1:34502/1‑151 (MQ=255) tGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGccc > 1:630716/1‑151 (MQ=255) aaaTGCCTAGATCGCCATGCATACCCCCCAAAATGTAATGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGccc < 2:658974/149‑1 (MQ=255) aaTGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCggg > 2:469487/1‑151 (MQ=255) ccTAGTTCGACATCAAAACCGCCCAAAAAGTAAGGGCGATAGTGCTCCCACTTGGTTTTACCGACGCCGTGCCCAACGACCCGCTGGGCGGGACGCAGGGGGGTGTCCCCGTGTTCGTCCGCCATGATCGTCACCGGATCGCCCGGGATc < 2:548639/150‑1 (MQ=255) ccTAGACCGCCAGCCATACCCCCCAAAATGTAAGGGAGATACTGCTCCCACAGCGTTTTATCTAAGCCGATTCCAGCAAGCAGCGGCGCGTGCCGCTGTGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGAt < 2:109740/149‑1 (MQ=255) aGATCCCCATGAATACCCCCCAAAATGAAAGGGAGATCCTGCGCCCACATCGTTTTATCAAAGCCGATTTCAGCCAGCGGCTGCGCGTGCCGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGAt < 1:446838/146‑1 (MQ=255) aGACCGCCATGCAAACCCCCCAAAATGTAAGGGAGATCCTGCCGCCACATCGTTTTATCAAAGCCGAGTTCACCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt < 2:15531/151‑1 (MQ=255) cGCCAGGCATACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCTAAGCCGATTGCAGCCCCCAGCTGCGCGGGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgtg < 2:305256/148‑1 (MQ=255) ccATGCATAACGCCCAAAATGTATGGGAGATCCTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGGGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa < 1:305167/151‑1 (MQ=255) aTGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa > 2:521469/1‑149 (MQ=255) aTGCATAACCCCCAAAATGTAGGGGAGATCCTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCCGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAg < 2:510243/151‑1 (MQ=255) aTGAATACCGCCCAAAATGTAAGGGGGATCCTGCGCCCACCTGGTTTTATCAAAGCCGAGTCCGGCCAGCAGCTGGGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa < 1:944739/150‑1 (MQ=255) gCATAACCCCCAAAATGTAATGGAGATACTGCTGCCACATGGTTTAATCTAAGCCAAGTCCACCCGGCGGCGGCCCGTGCCGCTCTGGGGAGCCCCCACGTCGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa < 2:548975/147‑1 (MQ=255) accccccaAAATTTAAGGGAGATCCTGCTGCCACATGGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa < 1:647264/146‑1 (MQ=255) aaCGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAt > 1:35200/1‑151 (MQ=255) ccccccAATGGTAGTGGAGATATTGCCGCCACGGCTTTTTATAAACGCGTATTCCGCCCGGCAGCGGCGAGTGCCGCTTTGGGGATTTCCAGCGTGGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa < 1:467516/146‑1 (MQ=255) cGCCCCAAATGTAATGGAGATCCTGCTGCCACAGCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCTGCGGGTGACGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtg < 1:259750/150‑1 (MQ=255) cAAATGTAAGGGAGATCCGGCGCCCACTTGTTTTTATCAAAGCCGAGTTCAGCCCCCAGCTGGGCGGGACGCTGGGGGGAGATCCCAGGTGCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga < 2:65768/150‑1 (MQ=255) aaaTGTAAGGGAGATCCTGCGGCCACATCGGTTTATCTAAGCCGAGTTCCGCCAGCCGCTGCGCGTGACGCTTGGGGGAGACCCCGCGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg < 2:524322/151‑1 (MQ=255) aagggAGATCCTGCTCCCACTTCGGTTTATCAAAGCCGATTTCAGCCAGCACCGGCGCGTGCCGCTTGGGGGAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt < 2:766858/145‑1 (MQ=255) aaTGGAGATTCTGTTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg < 1:774613/147‑1 (MQ=255) aTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt > 2:79701/1‑150 (MQ=255) tgcgcccACATCGTTTTTTCTAGGCCGATTCCAGCCAGCAGCGGCGCGTGCCGCTGTGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa < 1:684809/140‑1 (MQ=255) tgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg > 2:528806/1‑151 (MQ=255) tgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGCAATACCGATCAACGTgg > 2:513560/1‑151 (MQ=255) ccacagggttttATCAAAGCCGAGTTCGCCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg < 2:654016/139‑1 (MQ=255) ccACATCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGt < 1:25075/147‑1 (MQ=255) ccACATCGGTTTATCTAGGCGGATTTCAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGATCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg < 1:663696/148‑1 (MQ=255) gttttATCAAAGCCGATTCCAGCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGAc < 1:260746/147‑1 (MQ=255) tAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc > 2:343415/1‑147 (MQ=255) aggcggAGTTCACCCAGCAGCTGCGCCTGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc < 1:88728/146‑1 (MQ=255) aggCCGGGTTCAGCCGGCCGCGGCGCGGGCCGCCTTGGGGGCCCCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc < 1:11053/144‑1 (MQ=255) ccGAGTCCAGCCAGCAGCTGCGCGGGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCa < 2:672916/145‑1 (MQ=255) aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg > 2:934924/1‑150 (MQ=255) gtggcgctggACGCTGTGGGGCGCCCCCAGGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTg < 1:848142/139‑1 (MQ=255) tGCGCGTGACACTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc > 2:533303/1‑150 (MQ=255) gcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATa > 1:739441/1‑150 (MQ=255) gggccgctggggggAGACCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAc < 1:361501/140‑1 (MQ=255) gTGCCGCTCGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCaa < 1:538118/144‑1 (MQ=255) tGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc > 2:70839/1‑151 (MQ=255) aCGCTGTGGGGGGGCCCCCTGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGTGTCCCAAACGTCGGAGAATAAACTGCAACAt < 2:216329/144‑1 (MQ=255) cgctggggggAGCCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc < 2:375015/142‑1 (MQ=255) cGCTCTGGGGAGCCCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCgg < 1:159683/150‑1 (MQ=255) tctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtc > 2:906155/1‑151 (MQ=255) tGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtctc > 2:914680/1‑151 (MQ=255) ggggAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt < 2:566514/146‑1 (MQ=255) agaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGGGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta > 2:490482/1‑148 (MQ=255) | AACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTA > NC_000913/3705516‑3705805 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |