Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 46.3% | 93.3 / 11.2 | 55 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (27/2); new base C (0/25); total (27/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.17e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984291‑2984591 | cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca < 2:453151/151‑1 (MQ=255) cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca < 1:925586/151‑1 (MQ=255) gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac < 1:721347/151‑1 (MQ=255) tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATCTTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt > 2:549305/1‑151 (MQ=255) tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc > 2:470609/1‑151 (MQ=255) gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgcc < 2:13772/151‑1 (MQ=255) aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgcca < 2:241689/151‑2 (MQ=255) aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc > 1:249548/1‑150 (MQ=255) gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt > 1:1053534/1‑151 (MQ=255) ggggGCCAAAAAAATTTAGGTCGATTTGCCAGCCATCATTATTCAACAGATTAAGGATGGCAGATATTTTTTGGTAGGGAGAAAAGGAAGTGTTATTCGTCACGTTAAATGGTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTg < 1:591089/151‑1 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 2:386355/1‑151 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt > 1:635816/1‑151 (MQ=255) aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg > 1:94211/1‑151 (MQ=255) ttcccagcaaTCTTTATACAACAGATTATTGGCGCAATAAATTAAGAGGAAAGGAGAAAATGAAAGTTTTTTCGTCACTTGAAGGGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt < 2:397903/147‑1 (MQ=255) ccacccaTAATTATTAAACAGATTAATGACGGCAGATATTTATGGGGAAGGAGGAAATGAGGAGTTTTTCGTCACTTAAAATGGACAATGCCACCCCCGGACCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt < 1:194013/146‑1 (MQ=255) caagccaTCATTATCAAACAGATAAATGACGGCAATTTTTTATCGGAAGGAAGAAAATGCAGAGTTTTTCGTCTTTTAAAAGGTTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc < 2:1050470/147‑1 (MQ=255) atcttTATTAAAAAAATTAATGACGGCAGATTTTTATCGGTAGGGAGAAATTGCAGAGTTTTTGGTCGCGTACAATGTTAAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt < 1:249169/150‑1 (MQ=255) catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt > 2:485407/1‑151 (MQ=255) ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 1:869559/1‑151 (MQ=255) acAACAGATTAATGACGGAAGATATTTATCGGAAAGGGGGAAATGCAGAGTTTTTCGTCACTTTCAAGGGTCAGTGCCACCCCTTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 2:934943/150‑1 (MQ=255) aaaaaaTTAAGGCCGGAAAATTTTTACCGGTAAGGAGAAAATGCAGATTTTTTCTCCCCTTTCATGGGTAACTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga < 2:265637/146‑1 (MQ=255) aaaaGATTAATGACGGCAATTATTTATGGGAAAGGAGAAAAGGCAGAGCTTTTCGTCACTTTAAAGGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc < 2:90460/145‑1 (MQ=255) aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa < 1:767137/149‑1 (MQ=255) aaCAGATAAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga < 1:856198/151‑1 (MQ=255) aaCAGATAAATGACGGCAAATTTTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga < 1:687475/151‑1 (MQ=255) aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag > 1:297680/1‑149 (MQ=255) aaTGCCGGAATATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa < 1:719647/150‑1 (MQ=255) aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 2:114565/1‑149 (MQ=255) aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 1:529280/1‑149 (MQ=255) aaTGACGGCAGATATTTACGGGTAAGGAGAAAATGAAGAGTTTTTCGTCACGTTAAATGGTCAAGCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg < 1:612831/151‑1 (MQ=255) gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa > 2:1066266/1‑150 (MQ=255) ggaaaattttttTCGGTAAGGAGAAAAGGCAGAGTTTTTCGTCACGTTAAAGGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt < 2:1062330/141‑1 (MQ=255) ggCAGATATTTATGGGTAGGGAGAAAAGGCAGAGTTATTCGTCACTTTCAATGTTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc < 1:552128/150‑1 (MQ=255) gCATATTTTTACGGGAAAGGAGGAAATGCAACGTTTTTTGGCACTTTCGGTGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc < 2:522756/141‑1 (MQ=255) tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc > 2:596877/1‑151 (MQ=255) tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc > 1:570397/1‑151 (MQ=255) tttATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc < 1:397261/149‑1 (MQ=255) ggaaagaagAAAATGCAGAGCTTTTCTTCACGTTCAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg < 1:619515/145‑1 (MQ=255) ggTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc < 1:452721/151‑1 (MQ=255) gTAAGGAGGAATTAGAGGGATTTACGACACTTTCAGTGTTAAATCCCCCCCCGTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa < 1:1052297/146‑1 (MQ=255) aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt > 1:307823/1‑150 (MQ=255) ggggaaaaTGCAGAGCTTTTCGTCACGTTAAAGGTTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt < 1:971518/145‑1 (MQ=255) aggaAAGGCAGACTTATTCGTCGCGTTAAAGGGTCAACCCCCCCCCTCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg < 1:824261/150‑1 (MQ=255) gaaaaTGGAGAGTTATTCGGCACATTCAAGGGACAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg < 2:788835/149‑1 (MQ=255) aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg > 1:272326/1‑150 (MQ=255) aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc > 2:152536/1‑151 (MQ=255) agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:537825/1‑149 (MQ=255) agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:51592/1‑149 (MQ=255) agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 1:89329/1‑149 (MQ=255) cgatgggtttattGGTCAATGCCACCCCCTCCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc < 1:116433/137‑1 (MQ=255) cctgTTCATGGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc < 1:1030504/148‑1 (MQ=255) cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGCGACCGAGTGCGATCGCGcc > 1:193832/1‑151 (MQ=255) cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCACCTATCGGGTGACCGAGTGCGATCGCGcc > 1:591366/1‑151 (MQ=255) cACGTTCAATGGTCAAAGCCACACCCTGGCCCCCGCCGATACAAAGCGGTGCCAGTCCTTTGCGGGCATTACGTATCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGTGTGCGATCGCGcc > 2:987153/1‑151 (MQ=255) gTTCAATGGTCACTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGCGAAACCAGGATTCGGCAAACGGAAGCGCCTATCGGGAGACCGCGTGCGATCGCGccacc > 2:176555/1‑151 (MQ=255) tGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac > 2:9133/1‑151 (MQ=255) gTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat > 1:232582/1‑151 (MQ=255) acgCCACTCCCCTGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGACTATCGGGTGACCGAGTGCGATCGCGCCTCcattgacattga < 2:190463/148‑1 (MQ=255) acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc > 2:440158/1‑151 (MQ=255) | CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC > NC_000913/2984291‑2984591 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |