Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 53.7% | 46.9 / 87.7 | 54 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (24/1); new base C (0/29); total (24/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.78e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-10 |
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC > NC_000913/2984291‑2984589 | cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCcaca < 2:721894/151‑1 (MQ=255) tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCt > 1:582141/1‑151 (MQ=255) tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑cc > 2:755068/1‑151 (MQ=255) aTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATTGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑cccc > 1:45684/1‑151 (MQ=255) aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCccgccg > 2:793213/1‑150 (MQ=255) aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGtt > 2:36291/1‑151 (MQ=255) aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGt > 2:869602/1‑150 (MQ=255) aGCCAAGAAAATTTATGTCGACTTGCCAGCCATCATTTTTCAAAAGATTAATGCCGGCAAATTTTTCTCGTTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTAATT‑CCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGcc < 2:37922/150‑1 (MQ=255) aaaaatTGTTTGTCGACTTGCCACCCATCTTTATTAAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGGAAATGCAGGGCTTTTCGTCACTTTCATTGGTAAAT‑CCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCt < 2:891110/149‑1 (MQ=255) aaaaTGTATGTCGCCTGGCCACCAATAATTATCCAACAGATTAATGACGGCAAATATTTATGGGAAAGGAGAAATGGGAGAGTTTTTCTTCATGTTAAAGGGTAAAT‑CCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCtt < 1:659572/150‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 1:807895/1‑151 (MQ=255) aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAAT‑GCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc < 2:1117423/151‑1 (MQ=255) tatgtCGACTGGCCCGCCATCTTTATCAAAAAAATAAATGCCGGAAGATTTTTCGCGGTAGGGAGAAAATGCGAGGCTTTTCTTCACTTTCACTGGTAAAT‑GCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg < 2:76261/148‑1 (MQ=255) cgtcTTGCCACCCATCATTATTCAAAAAATTAAGGGCGGAAGATTTTTATGGGAAGGGAGAAAATGCAGAGTTTTTCTTCCCTTTAATTGGCAAAT‑CCCCCCCCGTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt < 1:230563/146‑1 (MQ=255) cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc > 1:662407/1‑151 (MQ=255) gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg > 1:714906/1‑151 (MQ=255) gACTTGCCAGCCATCATTATCCAACAGATAAATGGCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGACTTATACTTCACTTTCATGGGTAAAC‑CCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc < 1:789249/150‑1 (MQ=255) ttcccagccaTAATTATTCAACAGATTAATGCGGAAAGATTTTTACCGGAAAGGAGAAAATGCAGAGCTTTTCTTCCCGTTCAATGGTCATT‑CCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc < 1:65552/144‑1 (MQ=255) ttcccacccaTCTTTATTCAAAAAATAAATACCGGAATATATTTATGGGAAAGGAGAAAACGAAGACTTTTTCGTCAGGTTAAGGGGCCATT‑GCCCCACCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt < 1:1040644/146‑1 (MQ=255) tgcccacccaTTATTATCCAACAGATTAATGAGGGCATATTTTTAGGGGTAAGGAGGAAAGGCAGAGCTTTTCGTCGCTTTCAGGGGTAAAT‑GCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt < 1:1107692/146‑1 (MQ=255) gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc > 1:766737/1‑150 (MQ=255) ccacccaTAATTATTCAACAAATAAATGCCGGAAGATATTTATGGGTAAGGAGGAAATGCAGATCTTTTCGTCCCGTTCAATGGTCACC‑GCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc < 1:458955/149‑1 (MQ=255) cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc > 1:265047/1‑151 (MQ=255) ttattCAACAGATTAATGCCGGCAGATATTTACCGGTAAGGAGAAAATGCAGACCTTTTCGTCACGTTAAATGGTCAAT‑CCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt < 1:1161262/144‑1 (MQ=255) cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 2:903900/1‑149 (MQ=255) aaaagtttAATGACGGCAAATATTTATCGGGAGGGAGGAAATGCAGACTTTTTCGTCACGTTCAGTGGCTAAT‑CCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga < 2:625078/145‑1 (MQ=255) aTTAAGGACGGCAGTTATTTATGGGAAGGGAGAAAATGCAGACTTATTCGTCTCGATAATTGGTCAAT‑CCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc < 1:858157/151‑1 (MQ=255) aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 1:1156304/1‑149 (MQ=255) aaTGACGGCAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGCTCAAGGGTAAAT‑GCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc < 2:215280/149‑1 (MQ=255) aaTGACGGAAGATATTTATGGGTGAGGAGAAAATGAAGAGTTTTTCGTCTTTTTCAAGGTTAAAC‑GCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc < 2:38330/149‑1 (MQ=255) ggaaaatatTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAAT‑CCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc < 2:824119/145‑1 (MQ=255) ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTAAAT‑CCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg < 2:534162/151‑1 (MQ=255) agttatTTAGCGGTAGGGAGGAAATGCAGGCCTTTTCGTCTCGTCAAAGGGTAAAT‑CCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg < 1:692476/146‑1 (MQ=255) tttATCGGAAAGGAGAAAATGCAGATTTTTTCGTCACTTTAATTGGTCAAC‑GCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc < 1:1153662/148‑1 (MQ=255) ggaagggagaaAATGCAGATCTTTTCGTCCCTTTCAAGGGTCAAT‑CCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc < 1:137746/144‑1 (MQ=255) ggaaaGGAGAAAATGAAGAGCTATTCGTCACGTTAAATGGTCAAT‑CCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc < 2:738952/146‑1 (MQ=255) ggTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAAT‑CCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc < 2:652789/149‑1 (MQ=255) tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 2:1141580/1‑151 (MQ=255) aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 2:1084429/1‑150 (MQ=255) aaGGAGAAATTGCAGAGTTTTTCGTCAC‑TTTAATGGGTAAA‑GCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt < 2:347236/150‑1 (MQ=255) ggaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc > 1:416298/1‑151 (MQ=255) gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg > 2:136723/1‑150 (MQ=255) aaaTGCAGAGCTTTTCGTCACGTTAAAGGGTAAAT‑CCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg < 2:681290/149‑1 (MQ=255) aaaGGCAGAGTTATTCGTC‑CCTTTAATGGGTAATCCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg < 2:244221/151‑1 (MQ=255) aaTGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg > 2:575840/1‑150 (MQ=255) aaTGAAGAGCTTTTCGTCACGTTAAATGGTCAAT‑CCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa < 2:995287/151‑1 (MQ=255) aTGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCTCCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc > 2:78480/1‑151 (MQ=255) gagtttttCGTCCCTTTCAAGGGTCAAT‑CCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc < 2:1097219/141‑1 (MQ=255) gagCTTTTCGTCACGTTCAAGGGTCAAT‑GCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt < 1:182664/151‑1 (MQ=255) cGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac > 1:117586/1‑151 (MQ=255) gTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGACAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc > 1:77405/1‑151 (MQ=255) aagggTCAAT‑GCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg < 1:17665/148‑1 (MQ=255) aTGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg > 2:260986/1‑150 (MQ=255) cAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg > 2:162588/1‑151 (MQ=255) ccACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg > 2:648641/1‑151 (MQ=255) acacCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc > 1:995164/1‑151 (MQ=255) acacCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc > 1:890873/1‑151 (MQ=255) | CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC > NC_000913/2984291‑2984589 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |