Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,014,023 | 0 | T | G | 47.4% | 56.1 / 10.3 | 39 | G255G (GGT→GGG) | pqiB | paraquat‑inducible, SoxRS‑regulated MCE domain protein |
Reads supporting (aligned to +/- strand): ref base T (1/19); new base G (18/0); total (20/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.13e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG > NC_000913/1013874‑1014162 | gCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGt < 1:903365/150‑1 (MQ=255) gACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAg < 1:26331/151‑1 (MQ=255) aCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGA‑ttt > 2:257117/1‑148 (MQ=255) aCTGGTGACCAACAATGTTCGCTGCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAg < 2:501304/151‑1 (MQ=255) cTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGGCAGCTTTTATGTGGCGGAAGGGCTGGA‑TTTACGGc > 1:360175/1‑148 (MQ=255) acaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCAcc < 2:218743/151‑1 (MQ=255) tGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGGCCGCTTTGATGGGCCGGAAGGACGGGG‑TTTAGGGCAGCCAGGGGCACCGaa > 1:780950/1‑150 (MQ=255) gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGA‑TTTAGGGCAGCCAGTGGGACCGaaa > 1:875842/1‑150 (MQ=255) gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGaaaa < 2:494735/151‑1 (MQ=255) ttCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGA‑TTTAGGGCAGCCAGGGGCACCGAAAAAAGCtt > 2:322414/1‑151 (MQ=255) tCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGTTGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCATTGTATTGCCGGACAGGCTGGG‑TTTTGGGTAGCCTGGGGGACCGaaaaaaacgt > 1:375813/1‑144 (MQ=255) cTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGGTGGG‑GGTTGGGGAGCCAGTGGCACCGAAAAAAACtt > 2:462453/1‑149 (MQ=255) gAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGAAGGGCTGGG‑TTTAGGGGAGCCCGTGGCACCGAAAAAAACTTTTGtt > 1:881714/1‑151 (MQ=255) aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGGTTTGTTGTGCCGGGAGGGCTGGG‑TTTTGGGGAGCCAGTGGCACCGAAAACAGCCTTCGtttt > 2:497152/1‑151 (MQ=255) aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtttt < 1:579039/151‑1 (MQ=255) aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTAAGCCTTGTGGTGGCGGGAAGGCTGGG‑GTTTGGGCGGGCAGGGGCCCCCGAAACAGCTtttgttttgt > 2:105686/1‑150 (MQ=255) tGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTGATTTGGCGGAAGGGCTGGG‑TTTTGGGGAGCCAGGGGGACCGCAAACAGCTTTCGGTTTGTTTGATGATaagaag > 2:724130/1‑148 (MQ=255) gATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGc < 1:415916/150‑1 (MQ=255) gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGG‑TTTTGGGCAGCCAGGGGCACCGAAAACAGCTTTCGTTTTTTTTGATGATCAGAAGAGCATTTAgg > 1:817357/1‑149 (MQ=255) gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt < 2:415704/151‑1 (MQ=255) tCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTTATGTGGCGGGAGGGCTGGG‑TTTTGGGGAGCCCGTGGGACCGGCAAAAACCTTTGTTTTTTTTGATGATCAAAaaaggaattagggtt > 2:847285/1‑140 (MQ=255) gggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGt < 1:952573/151‑1 (MQ=255) ggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTCTCCGATTTTAGTTGGGGGGAGGGCCGGAGTTTTGGGGAGCC‑GGGGGACCGGAAACAACTTTCGGTTTGTTTTTTGGTTCGAAAAGCCGTCAGGAGTCGtt > 2:575298/1‑148 (MQ=255) gtgGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGG‑TTTAGGGGAGCCAGTGGCACCCAAAACAGCTTTTGTTTTGTATGATGATCAGAAGAGCATTCAGGAATCGTTGTACACCgct > 1:260052/1‑148 (MQ=255) gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGa < 2:115499/150‑1 (MQ=255) gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGAt < 1:220596/151‑1 (MQ=255) cTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGTTGTGCCGGAAGGGCTGGG‑TTTTGGGGAGGCAGTGGGACCGGAAAACGCTTTCGTTTTGTATTGTGATTAGAAGAGGATTCAAGATTCGTTGTTCACCGATTACATTGAt > 2:705421/1‑148 (MQ=255) tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGA‑TTTTTGGGAGCCCGGGGGACCGAAAAAAGCTTTCGTTTTGTGTGATGATTAGAAAAGCATTCCGGATTCGTTGTACACCGATCACAtttttt > 2:187539/1‑144 (MQ=255) tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGGCTGGG‑TTTAGGGGAGCCAGTGGGACCGAAAACAGCTTTTGTTTTTTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACAACGATCACATTGAtt > 2:492459/1‑148 (MQ=255) cAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGG‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGtt > 1:800331/1‑151 (MQ=255) cGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCtt < 2:505999/151‑1 (MQ=255) gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt < 1:588156/151‑1 (MQ=255) gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt < 2:759643/151‑1 (MQ=255) gAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGAt < 2:87576/151‑1 (MQ=255) gTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTc < 1:518104/151‑1 (MQ=255) gtCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgc < 1:226620/150‑1 (MQ=255) tCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgcg < 2:609436/150‑1 (MQ=255) gCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCTATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCt > 1:742007/1‑151 (MQ=255) aTGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAcc < 2:380890/151‑1 (MQ=255) gtgCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCgg < 1:422008/151‑1 (MQ=255) | GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG > NC_000913/1013874‑1014162 |
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: ‑ |