Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 45.6% | 90.2 / 15.6 | 57 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (27/4); new base C (0/26); total (27/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.11e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984291‑2984591 | cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca < 1:482730/151‑1 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg > 2:133249/1‑151 (MQ=255) cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccc < 1:465583/151‑1 (MQ=255) gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccc < 1:700507/150‑1 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg < 1:947703/151‑1 (MQ=255) ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGAt > 1:629207/1‑150 (MQ=255) tGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCg > 2:823737/1‑151 (MQ=255) ccaaaaaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGCCGGCAGA‑TATTTATCGGAAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTAAAGCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCa < 2:348623/144‑1 (MQ=255) ccaaaaaaaTGTATGTCGACTTGCCAGCCATCATTATCCAACAGATTAATGACGGCATA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGt < 1:39146/146‑1 (MQ=255) cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc > 2:408625/1‑151 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc < 1:802755/150‑1 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc > 2:93734/1‑150 (MQ=255) aaaaTGTTTGTCGACTGGCCAGCCATCATTATTCAACAAATTAATGCGGGCAGA‑TATTTATCGGGAAGGAGGAAATGCAGAGCTATTCGTCACGTTCATTGGTCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTc < 2:298866/147‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 1:498890/1‑151 (MQ=255) aaaaTGTATGGCGACTGGCCAGCCATAATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCGGGGCTTTTCGTCCCTTTCAATGGTCAACGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCtt < 2:392077/150‑1 (MQ=255) aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc > 2:641724/1‑151 (MQ=255) aaTGTATGTCGACTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAGGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc < 1:660415/151‑1 (MQ=255) atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg > 2:276028/1‑151 (MQ=255) tatgtCGTCTTGCCAGCCATAATTTTTAAACAAATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAAGGCAGAGCTATTCGTCCCGTAAAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg < 2:769145/150‑1 (MQ=255) tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt > 1:785081/1‑150 (MQ=255) taCTTGCCACCCATCTTTATTCAACAGATTAATGCCGGCATA‑TTTTTATGGGTAAGGAGAAAATCCAAACTTATTCGTCACGTTCAGTGGTCATTCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa < 1:561932/148‑1 (MQ=255) ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt > 1:556349/1‑151 (MQ=255) ccagccaTCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc < 1:595039/151‑1 (MQ=255) ccagccaTCATTATTAAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc < 1:824497/151‑1 (MQ=255) caccattcatTATCCAACAAATGAA‑AAGGGGAGATTATTTAGCGGAAAGGGGAAAATACGGAGCTTTTTGTCCGGTTCAATGGTCATTGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa < 2:371094/143‑1 (MQ=255) catTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt > 2:577528/1‑151 (MQ=255) attattCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 1:781403/1‑151 (MQ=255) ttattCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg < 1:52832/150‑1 (MQ=255) ttattCAAAAGATTAATGCCGCCAGA‑TTTTTAGGGGTAAGGAGAAAATGAAGAGCTTTTCGTCGCGTTAAAGGGTCAAGGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa < 1:17412/148‑1 (MQ=255) taattAAACAGATTAAGGCCGGAAAA‑TATTTATCGGTAAGGAGAAAATGCAAAGTTATCCGTCACTTACAATGATAACGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc < 2:650754/145‑1 (MQ=255) attCAACAGATTAATGGCGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAAGGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc < 2:629144/151‑1 (MQ=255) attCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc > 1:77699/1‑151 (MQ=255) aaCAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag < 2:351207/150‑1 (MQ=255) aaCAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 1:853791/149‑1 (MQ=255) aTTAATGCCGCCAGA‑TTTTTGTCGGAAAGAAGGAAATGAAGAGCTATTCCTCCCGTTCAATGGTCATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga < 2:12821/148‑1 (MQ=255) aTTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc > 1:253195/1‑151 (MQ=255) tGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg > 1:240999/1‑150 (MQ=255) tGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg > 1:848466/1‑150 (MQ=255) ggCAGA‑TTTTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc < 1:595059/150‑1 (MQ=255) a‑tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc > 2:782447/1‑151 (MQ=255) tatTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa < 2:173788/148‑1 (MQ=255) tttATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa < 2:701427/147‑1 (MQ=255) tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCACTGGTCAATGCCACACTCTGTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGAGGACATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg > 2:478153/1‑150 (MQ=255) tttAGCGGTAAGGGGAAAAGGCAGAGCTTTTCGTCCCTTTCAATGGTCAAGGCCATCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa < 2:881890/146‑1 (MQ=255) tATGGGTAAGGAGAAAATGAAGAGCTTTTCGTCACGTTCAGGGGTAAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa < 1:408666/150‑1 (MQ=255) tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 1:212632/1‑151 (MQ=255) ggggaaaaTGCAAAGCTTTTCGTCAGGTAAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt < 1:42374/145‑1 (MQ=255) aaaTGCAGAGTTTTTCGTCACGTTCAATGGTCAAGCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg < 2:254767/151‑1 (MQ=255) gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:617106/1‑151 (MQ=255) gagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt > 1:371144/1‑151 (MQ=255) agCTTTTCGTCACGTTCAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg < 1:823820/151‑1 (MQ=255) gCTATTCGTCACGTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc < 1:57582/151‑1 (MQ=255) cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 2:860167/1‑151 (MQ=255) tATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg > 1:392393/1‑150 (MQ=255) cGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc > 1:1608/1‑151 (MQ=255) ttCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAGACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc > 2:145633/1‑150 (MQ=255) cAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt > 2:777500/1‑151 (MQ=255) aaTGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg > 1:28880/1‑151 (MQ=255) acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc > 1:882320/1‑151 (MQ=255) | CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGA‑TATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984291‑2984591 |
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: ‑ |