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.6% | 90.2 / 15.0 | 57 | F110V (TTC→GTC) | inaA | acid‑inducible Kdo/WaaP family putative kinase |
Reads supporting (aligned to +/- strand): ref base A (28/3); new base C (0/26); total (28/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.30e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA > NC_000913/2348995‑2349291 | aGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGaa < 2:834098/151‑1 (MQ=255) gATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAc > 1:849274/1‑151 (MQ=255) tatGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAccc > 2:377399/1‑151 (MQ=255) atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc < 1:84394/151‑1 (MQ=255) cATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATc > 2:350882/1‑151 (MQ=255) aTAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCt > 2:998438/1‑151 (MQ=255) aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt > 1:915863/1‑151 (MQ=255) aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt > 1:885001/1‑151 (MQ=255) acaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAg > 2:3104/1‑149 (MQ=255) ccATGCTGACAATTAATGATATCAATTTTCTTAAACGCCGGCCCCACTGCTTCCAA‑CATGGTTTGCCTTATTTCGTAAAAAAAAGGCGATCCTGAATGCTGGGCATACCAGTAAGCAATGCTGATAACCCCCGCCATATCTTCAGTCa < 2:440817/148‑1 (MQ=255) ctgaataTTAATGCTATGCTTTTTCTAAAACGCCAGCCCCACTGTTTTAAA‑CAGGGCTTCCCGTACTTCGTCAGAAAAGGGCAATACGGCATGCGGGGCATACCAGCAAGAAAGGCTAATAAACCCCGCCATATCTTCAGTCACCAGCa < 2:594526/143‑1 (MQ=255) cTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAAc > 1:638479/1‑151 (MQ=255) gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg > 2:206883/1‑151 (MQ=255) gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg > 1:422134/1‑151 (MQ=255) ttattccTATGTATTTTCTAAACCGCCCGCCCCCGTGTTTACAC‑CTTGGCTTGCCGTCCTTCTTCAAAATAGGGCGTTACTGTTTGCTGGGCAAACCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCa < 2:640969/136‑1 (MQ=255) tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc > 2:1004346/1‑151 (MQ=255) aTGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAccc > 1:287742/1‑151 (MQ=255) tattttttCTAAAACCCCAGCCCCCTTTTTTTAAA‑CATGGTTTGCCGTCCTTCGAAAAAAAAGGGCGATCCGGCGTGCGGGCCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt < 2:25364/141‑1 (MQ=255) ttttCTAAACCGCCCGCCCCACTGCTTTCAA‑CGGGGCTTGCCGTACTTCGACAAAAAAAGGCGATACTGCATGCTGGGCATCCCAGCCACCAAGGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCaa < 2:566849/150‑1 (MQ=255) aaccGCCAGCGCCCCTGTTTTCAA‑CATGGCTTGCCGTACTTCGTCAAAATAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc < 1:357265/145‑1 (MQ=255) aaCGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCg > 1:793070/1‑151 (MQ=255) gcCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa > 1:1014085/1‑150 (MQ=255) cACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa > 2:232702/1‑150 (MQ=255) ttttttCAC‑CAGGGTTTCCCTTTCTTCGAAAAAAAAAGGCGATCCTGCTTGCGGGACATACCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaaa < 1:76538/146‑1 (MQ=255) tttCAA‑CATGGTTGCCCTTACTTCGTAAAAAAAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt < 1:978064/150‑1 (MQ=255) tttCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATc > 1:333680/1‑151 (MQ=255) ttCAACCATGGCTTCCCGTCCTTCGCAAAAAAAGGGCGATCCTGCAGGCGGGCCACCCCAGCCAACAAGGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACg < 2:56870/148‑1 (MQ=255) tcaaa‑CATGTCTGGCCGTACTTCGTAAAAATAAGGCGATACTGCATGCTGGACACACCAGTCATCAATGCTGAAAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt < 2:652063/146‑1 (MQ=255) tGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCAc > 2:973663/1‑151 (MQ=255) tcccGTCCTTCGTAAAAAAAAGGCAATTCTGCATGCTGGCCATCCCGGCCAAAACTGTTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATg < 2:385947/148‑1 (MQ=255) ccGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGac > 1:931195/1‑150 (MQ=255) gTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACAcc > 1:874071/1‑151 (MQ=255) cataaGAAAAGGGCGATTCTGCTTGCTGGGAATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCg < 2:287675/142‑1 (MQ=255) aGAAAAAGGCGATCCTGCATGCGGGCCATCCCAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGtt < 1:350969/149‑1 (MQ=255) aaaaagtcatcgcTGCAGGCGGGGCATCCAAGCCAAAAACGCGAAAGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTc < 2:745359/136‑1 (MQ=255) aaTAGGGAGATTCTGCAGGCGGGCAACCCCAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAg < 2:547972/151‑1 (MQ=255) aaTAAGGCGATCCTGCATGCGGGGCATCCCAGCCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAg < 1:45288/151‑1 (MQ=255) aaTAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc > 2:723105/1‑149 (MQ=255) ggCTATACTGCATGCTGGACATCCCAGTCAGCAATGCTAATACCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCt < 2:921492/150‑1 (MQ=255) gCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCttt > 2:1048669/1‑151 (MQ=255) gCGATACTGCATGCTGGCCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCt < 1:165689/149‑1 (MQ=255) gATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTa > 2:1032953/1‑150 (MQ=255) tgcgtgcgGGGCATCCCAGTCACCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTa < 1:232748/149‑1 (MQ=255) cgttttGGGCATACCCGTCGAAAACGATTATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTa < 2:895085/142‑1 (MQ=255) cATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGc > 1:133833/1‑151 (MQ=255) ttgctgcGCACCCCAGCCAACAATGCTAATGCCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAtt < 2:859010/144‑1 (MQ=255) ggccATACCAGTCAGCAATGCTAATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAAcc < 1:766806/146‑1 (MQ=255) tACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACgg > 1:608667/1‑150 (MQ=255) ccAGTAAGAAATGCGGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGa < 1:206926/149‑1 (MQ=255) gtcaaaaaTGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATc < 1:647429/145‑1 (MQ=255) aGCAGTCCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAAt < 2:767415/147‑1 (MQ=255) gCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTg > 2:777756/1‑151 (MQ=255) aaggcagatgaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTgg < 1:560208/144‑1 (MQ=255) cttcTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGt < 2:666970/147‑1 (MQ=255) tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg > 1:984516/1‑151 (MQ=255) tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg > 1:346496/1‑151 (MQ=255) tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg > 2:891308/1‑151 (MQ=255) atgaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGa < 1:998761/151‑1 (MQ=255) | AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA > NC_000913/2348995‑2349291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |