Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 42.1% | 96.3 / 20.1 | 57 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (32/1); new base C (0/24); total (32/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.32e-15 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC > NC_000913/2984292‑2984585 | gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCacac > 2:978742/1‑151 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCGATGGTCAATG‑CCACAcc > 2:1086597/1‑151 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc > 2:1232977/1‑151 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc > 2:1129929/1‑151 (MQ=255) atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAcc > 1:511604/1‑150 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCt > 2:846074/1‑150 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACCCCCTGGCCCccgccg < 1:862474/151‑1 (MQ=255) aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACa > 2:655139/1‑150 (MQ=255) aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGCTACaa > 1:989496/1‑149 (MQ=255) aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGtt > 2:796002/1‑151 (MQ=255) aGGGAGCCAAGAAAATGTATGTCGACTTGCCAGCCATAATTATTCAACAGATTAATTACGGCAGATATTTATCGGTGAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACACTAGCCCCCGCCGATACAAAGCg > 1:512065/1‑149 (MQ=255) gggccccaaaaaaaTTTATGTCCACTTCCCACCCATTATTATTCAAAAGATAAATGCCGGAAAATTTTTAGCGGAAAGGAGGAAATGCAGATCTTTTCACCACGTTAAATGTTAAACC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGtt < 1:889891/140‑1 (MQ=255) gggAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGTTATTCGTCCCGTTAAATGGTCAATC‑CCACCCCCGGGCCCCCGCCGATACAAAGCGTTg < 2:1093220/151‑1 (MQ=255) taGCCCAGAAAATGTATGTCCACCTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGc > 1:578277/2‑150 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt > 1:369792/1‑151 (MQ=255) aaaaTGTATGTCGATTTCCCAGCAATTATTATTAAAAAGATTAATGCGGGCAGATTTTTATGGGAAGGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTAAAGG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTc < 2:734372/147‑1 (MQ=255) aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc > 2:1009642/1‑151 (MQ=255) gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg > 1:156786/1‑151 (MQ=255) ttGCCAGCCATCTTTATACAACAGATAAAGGACGGCAGATTTTTAGGGGTAGGGAGAAAAGGGAGAGCTTTTCGCCACGTTAAATGGTAAATG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa < 1:194366/146‑1 (MQ=255) tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc > 2:392824/1‑151 (MQ=255) gccagcctTCATTATCCAACAGATTAATGACGGCAGATATTTTTCGGTAGGGAGGAAAGGCAGAGTTTTTCGTCACGTT‑AATGGGTAATGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc < 2:576358/150‑1 (MQ=255) ccagcaaTCATTATAAAACAAATTAATGGCGGGAGATATTTATCGGAAGGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAGTC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc < 2:334255/149‑1 (MQ=255) ccacccaTAATTATCCAACAGATTAAGGACGGCATATTTTTATCGGTAGGGAAGAAAGCCAGCGCTATTCGTCCCGTTCAATGGTCATTC‑CCCCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa < 2:25043/150‑1 (MQ=255) ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt > 2:131485/1‑151 (MQ=255) attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 2:794048/1‑151 (MQ=255) ttattaAACAGATAAATGAGGGAAGATTTTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAATGGTCAACG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc < 2:703999/147‑1 (MQ=255) attAAACAGATTAATGACGGAAAATTTTTAGGGGTAAGGAGAAAATGCAGAGTTTTTCG‑CCCCTTTAAAGGGAAATGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa < 2:1030670/150‑1 (MQ=255) aaaagttaaaggagggCAAATTTTTAGCGGAAAGGAGAAAAGGAAGAGTTTTTCTTCACTTTAAATGGTAAACC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa < 1:899407/131‑1 (MQ=255) aCAGATAAAGGACGGAAATTATTTACGTGTAGGGAGGAAATGCAGAGCTTTTCGTCCCGTTAAATGGTCAATG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt < 1:1301399/142‑1 (MQ=255) aTTAATGCCGGAAAATTTTTATCGGAAGGGAGGAAAGGCAGTGCTTTTCGTCACGTTAAATGGTCATTC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa < 2:149064/150‑1 (MQ=255) ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 1:1236285/1‑151 (MQ=255) tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg > 1:223945/1‑150 (MQ=255) gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt > 1:239441/1‑151 (MQ=255) cGGCAGTTATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGAAAAGC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc < 2:263278/151‑1 (MQ=255) cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 2:343215/1‑151 (MQ=255) atatttttATAGGAAAGGAGAAAAGCCAGATTTTTTCGTCCCTTCCATTGGTCATGC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg < 1:1159899/143‑1 (MQ=255) aGATTTTTATGGGAAAGGAGAAAATGCAGAGTTTTTCGTCACGTTCGATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg < 2:1005767/149‑1 (MQ=255) tatTTAGCGAAAGGGAGGAAATGAGGACTTTTTCTCCACGTTAAATGGTAAAGG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg < 1:75787/146‑1 (MQ=255) ttATCGGTAAGGAGAAAATGCAGAGCTTTTCGCCACTTTAAATGGTAAATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg < 2:875736/150‑1 (MQ=255) tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa > 2:1052659/1‑150 (MQ=255) aaGGAGGAAATGGATATTTTTTCGTCACGTTCAATGTTCAATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa < 1:1214760/151‑1 (MQ=255) aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 2:1147736/1‑150 (MQ=255) ggaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCTGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc < 1:639200/151‑1 (MQ=255) gagaaaggcaagtgCTTTTCGTCACGTTAAAGGGTCAAGC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg < 1:585544/138‑1 (MQ=255) aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg > 2:569433/1‑151 (MQ=255) cAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:1688/1‑150 (MQ=255) agggCTTTTCGTCCCTTTAAATGGTCAACC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa < 2:476770/147‑1 (MQ=255) agCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg > 2:1110838/1‑151 (MQ=255) tATTCTTCCCTTCCAATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg < 1:795972/148‑1 (MQ=255) aTTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa > 2:202883/1‑150 (MQ=255) cGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc > 2:618805/1‑151 (MQ=255) cACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 2:597617/1‑151 (MQ=255) cACGTTCAATGGTCAATC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc < 2:1033228/148‑1 (MQ=255) aCGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 1:743130/1‑150 (MQ=255) gTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc > 1:216441/1‑151 (MQ=255) tGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac > 2:405982/1‑151 (MQ=255) tGGTCAATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga < 2:592331/150‑1 (MQ=255) tG‑CCACACCCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc > 1:363170/1‑151 (MQ=255) | GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC > NC_000913/2984292‑2984585 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |