Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,014,023 | 0 | T | G | 43.1% | 85.5 / 14.3 | 58 | G255G (GGT→GGG) | pqiB | paraquat‑inducible, SoxRS‑regulated MCE domain protein |
Reads supporting (aligned to +/- strand): ref base T (10/23); new base G (25/0); total (35/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.13e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGG > NC_000913/1013884‑1014163 | aCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGG‑GCTGGATTTAgg > 2:409526/1‑151 (MQ=255) gACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGATTTTggg > 2:922901/1‑149 (MQ=255) gACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGGCAGCTTTTTTGGGCCGGAAGG‑GCTGGATTTTggg > 1:85442/1‑149 (MQ=255) ggTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCa < 1:596041/151‑1 (MQ=255) gACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGt < 1:1031903/150‑1 (MQ=255) gACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCCTTTATGTGGCGGGAGG‑GCTGGGTTTAGGGCAGCCAGTg > 2:839951/1‑151 (MQ=255) aCCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGTAGCTTTGATGTGGCGGGAGG‑GCTGGTTTTTGGGCAGCCAg > 1:293730/1‑148 (MQ=255) caacaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGTGAGCTTTTATGTGGCGGGAGG‑GCTGGGGTTTAGGCAGCCAGTgg > 1:941063/1‑149 (MQ=255) aacaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTGAACTTTTATGTGCCGGGAAG‑TGTGGGTTTTTGGCAGCCCGTggggc > 2:594165/1‑148 (MQ=255) acaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCCGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGCAGCCCGTGGGAcc > 2:824348/1‑151 (MQ=255) caaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGGAGCCAGTGGCAcc > 2:44688/1‑150 (MQ=255) caaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCg < 2:258034/151‑1 (MQ=255) aaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGGAGCCCGTGGCAcc > 2:498789/1‑149 (MQ=255) gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGaaaa < 2:186015/151‑1 (MQ=255) ttCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGATGGGCCGGGAGG‑GCTGGGTTTAGGGCAGCCAGGGGGACCGaaaa > 2:841804/1‑150 (MQ=255) cGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACa < 1:524277/150‑1 (MQ=255) cGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAg < 2:976445/151‑1 (MQ=255) tCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCCTTTGTGTGCCGGGAGGTTTTGGGTTTTGGGGGGCC‑GGGGCCCCCAAAACAACttt > 2:517056/1‑151 (MQ=255) ggAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGGAGG‑GCTGGGTTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCg > 2:727552/1‑150 (MQ=255) gATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGGCGGGAGG‑GGTGGGGTTAGGGGGGCCAGTGGCACCGGAAAAAGCTTTCGtttt > 2:607937/1‑149 (MQ=255) aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGGCGGAAGG‑GGTGGATTTTGGGGAGCCAGGGGGACCGAAAACAACTTTTGtttt > 1:152347/1‑148 (MQ=255) tAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGGAGGCAGGGGGACCGGAAAAAGCTTTTGtttttt > 2:588969/1‑147 (MQ=255) tGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTGCGCTTTTTTGGGCCGGAAGG‑GCTGGGTTTTGGGGAGCCAGGGGCACCGCAAAAAACTTTCGTTTTTTat > 1:10346/1‑148 (MQ=255) gCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGGCTGCTTTTATGTGCCGGTAGG‑TCTGGCATTAGGGGAGGCAGCGGGCCCGAAAAAAGCATTTGTtttgtttgtt > 2:670779/1‑147 (MQ=255) cATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATgatgat < 1:301021/151‑1 (MQ=255) cATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCCGAAGG‑TCTGGATTTAGGGCAGCCCGGGGCACCGAAAACAGCCTTCGtttttt > 2:337647/1‑141 (MQ=255) cGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGGAGCCCGGGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCag > 2:904472/1‑149 (MQ=255) tGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCCGCCTTTATGTGCCGGAAGG‑GCTGGGTTTAAGGCAGCCCGGGGGCCCGAAAAAAGCTTTTGTTTTTTATGATGATTagaa > 2:753651/1‑147 (MQ=255) gACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCa < 2:830370/146‑1 (MQ=255) cGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTTAGCTTTTTTGTGCCGGGAGG‑GCTGGGTTTAGGGGAGCCCGGGGGACCGGAAACAGCTTTCGTTTTTTTTGATTATCAGAAGA‑GCATTCAgg > 2:700365/1‑150 (MQ=255) gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGTCCGCGTTTATGTGCCGGAAGG‑GCTGGATTTAGGGGAGCCCGAGGCGCCGAAAACAGCTTTCGTTTTGTATGATGTTAAGCAGA‑GCATTCAgg > 1:780962/1‑149 (MQ=255) gcagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑GCTGGATTTAGGGCAGCCCGTGGGACCGGAAACAACTTTTGTTTTTTATGATGATCAGAAGA‑GCATTCAGGATTCg > 2:162187/1‑150 (MQ=255) cagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGtt < 2:394633/151‑1 (MQ=255) agGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTg < 1:1029751/151‑1 (MQ=255) aTGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGTGGCTTTTATGGGCCGGGAGG‑TCTGGGTTTTTGGGCGCCCGGGGGACCGAAAAAAGCCTTCGTTTTTTATGATGATTAAAAAA‑GGACTCAGGGATCGtt > 1:159167/1‑146 (MQ=255) gCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACAcc < 1:953918/151‑1 (MQ=255) cGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTTGCGGTGTCAGCTCTGATGTGCCGGAAGG‑TGTGGATTTAGAGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTAGGATGATCAGAAGA‑GCATTCAGGATTCGTTGTAcccc < 1:449183/150‑1 (MQ=255) gtgGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGAt < 2:143040/150‑1 (MQ=255) ggAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATcac < 1:387753/151‑1 (MQ=255) ggAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTTTATGATGATCAGAAGA‑GCATTCCGGGTTCGTTGTACACCCATca > 2:271144/1‑150 (MQ=255) gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGAt < 2:921850/151‑1 (MQ=255) ggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATTTGGCCGGAGG‑GCTGGGTTTTGGGCAGCCAGGGGCACCGAAAAAAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGAtt > 2:571939/1‑151 (MQ=255) ggCTCATTGACAACGCTGCTGAGTGGCGGGGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTAGGGGAGGCAGGGGGACCGGAAAAAGCTTTTGTTTTGTATGATGATCAGAAGAGGCTTCCAGGA‑TCGGTGTACACCGactcccattttt > 1:623750/1‑140 (MQ=255) tCATTGACAACGCTGCTGAGTTGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATc > 2:65445/1‑151 (MQ=255) tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGAAGG‑GCTGGGTTTTGGGGAGCCCGTGGCACCGAAAAACACTTTCGGTTTTTTTGATGATCAGAAGA‑GCATTCAAGGTTCGTTGTTCAACGGTCACATTGTTTAt > 2:509980/1‑150 (MQ=255) gACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGATTTAGGGCAGCCAGGGGGACCGAAAACAGCTTTCGTTTTTTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCCATCACATTTATTATCTGAt > 2:540328/1‑150 (MQ=255) gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt < 2:189045/151‑1 (MQ=255) gCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGt < 1:32136/151‑1 (MQ=255) cGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGt > 1:772396/1‑150 (MQ=255) cGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTa < 2:335369/151‑1 (MQ=255) ggTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAc > 2:272859/1‑151 (MQ=255) gCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCt < 1:867535/151‑1 (MQ=255) tGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCa > 2:152755/1‑150 (MQ=255) gATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAc > 1:986055/1‑151 (MQ=255) gATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAc < 1:922541/151‑1 (MQ=255) aTGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAcc < 2:816272/151‑1 (MQ=255) gCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg < 1:36789/150‑1 (MQ=255) gCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg < 1:286512/150‑1 (MQ=255) | ACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGG > NC_000913/1013884‑1014163 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |