Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 41.4% | 85.2 / 22.2 | 58 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (33/1); new base C (0/24); total (33/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.95e-15 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG > NC_000913/2984291‑2984587 | cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca < 2:288202/151‑1 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCt < 2:1137199/150‑1 (MQ=255) aaaaTGAGGGAGCCCAGAAAATGCATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa > 1:1079942/1‑150 (MQ=255) aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCACTGCCACACCCTGGCCCCCGCCGATACaa > 2:697818/1‑149 (MQ=255) gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt > 1:658010/1‑151 (MQ=255) ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc > 1:530125/1‑151 (MQ=255) gAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc > 2:666170/1‑151 (MQ=255) cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTAACGTTCAATGGTCAATGCCACAACCTGGCCCCCGCAGATACAAAGCGTTGCCAGTcc > 1:3318/1‑151 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt > 2:643329/1‑151 (MQ=255) aaatgttgttCGACTTGCCAGCCATCATTATTAAACAGATAAATGGGGGCAGATTTTTAGCGGTAGGGAGGAAAGCGAGAGCTATTCGTAACTTTAAATGTTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg < 1:1028608/142‑1 (MQ=255) gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg > 2:1036859/1‑151 (MQ=255) tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGAGGGCAGTTATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGCCACTTTCAATGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa < 2:956244/150‑1 (MQ=255) gtCGCCTTCCCACCCATTATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGGAAATGCAGAGTTTTTCGTCACGTTCAATGGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc < 1:335278/148‑1 (MQ=255) gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAATCGCTGCCAGTCCTTTGCGGGCATTACg > 1:267881/1‑151 (MQ=255) tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc > 1:123803/1‑151 (MQ=255) gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt > 2:1156980/1‑148 (MQ=255) gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc > 2:985845/1‑150 (MQ=255) cacccaTTATTATTCAAAAGATTAATGACGGAAGATTTTTAGCGGAAAGGAGAAAATGAAGATTTTTTCGTCCCGTTAAGTGGACAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc < 1:762004/147‑1 (MQ=255) acccaTCATTATAAAAAAAATAAGGGCCGGAAAATATTTATCGGTAGGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa < 2:693502/149‑1 (MQ=255) tcatTATTCAAAAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt < 2:41586/148‑1 (MQ=255) ttattCAACAGATTAATGCCGGCAGATATTTTTGGGTAAGGAGGAAATGCAGGGCTTTTCGTCACTTTCAATGTTAACTGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc < 1:46194/147‑1 (MQ=255) attacaaaaattattGACGGCGAATATTTGGCGGTAAGGAGAAAACGCAGAGTTTTTGGTCACTTTAAAGGGACAACCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc < 2:29966/142‑1 (MQ=255) attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc > 1:1055422/1‑151 (MQ=255) attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc > 1:573095/1‑151 (MQ=255) ttCAACAAATTAAGGACGGCAAATTTTTATCGGTAGGGAGAAAATGCAGAGTTTTTCGTCACTTTAAAGGTTCAACCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa < 2:793738/149‑1 (MQ=255) aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc > 2:403391/1‑148 (MQ=255) aaaTTAATGACGGAAGATTTTTATCGGAAAGGAGAAAAGGAAGAGTTTTTCGTCCCTTTAAATGTTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa < 2:1085430/149‑1 (MQ=255) gATTAATGAGGACAAATTTTTATCGGAAAGGAGAAAATGAAGACCTTTTCGTCCCTTTAAATGGTCACCGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga < 1:626483/147‑1 (MQ=255) ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 2:409366/1‑151 (MQ=255) taaaTGCCGGAAAATTTTTACCGGAAAGAAGAAAATGAAGAGTTTTTCGTCACGTCAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag < 1:542720/144‑1 (MQ=255) aaTGGCGGCAGATATTTATCGGAAAGGAAGAAATGCAGACCTATTCTTCACTTTCAATGATAAACCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc < 2:202027/149‑1 (MQ=255) gACGGCAGATATTTATCGGGAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGACCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt > 1:63897/1‑151 (MQ=255) cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt > 1:215593/1‑150 (MQ=255) ggCAAATATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGTTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg < 2:1120643/151‑1 (MQ=255) gCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg > 2:392011/1‑151 (MQ=255) aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc > 2:624807/1‑150 (MQ=255) atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg > 1:627710/1‑151 (MQ=255) ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCCTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc > 1:476974/1‑151 (MQ=255) gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg > 1:1018626/1‑151 (MQ=255) gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg > 1:554995/1‑151 (MQ=255) gagaaAATGCACTGTTTTTCTTCACATTCAATGGCCAACCCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAATCCAGGATTCGGCAACCGGAAg < 2:523272/142‑1 (MQ=255) ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg > 2:705655/1‑151 (MQ=255) ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCCTTGCGGGCATTAAGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGGCTATCggg > 2:856169/1‑150 (MQ=255) gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt > 1:212948/1‑151 (MQ=255) aggttttttCGTCCCGTTAAATGTTCAAGCCCCCCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg < 1:1161189/142‑1 (MQ=255) agagCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 1:1184674/1‑149 (MQ=255) agagCTATTCGTCACGCTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg > 2:809712/1‑149 (MQ=255) ttttCGTCACTTTAAATGCTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc < 1:444995/147‑1 (MQ=255) tttCGTCACTTTAATTGTTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc < 1:791831/147‑1 (MQ=255) ttCGTCCCTTTCAATGGCAAACCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt < 1:291350/150‑1 (MQ=255) ttCGTCACGTTCAATGGTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt < 2:594391/145‑1 (MQ=255) cACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc > 2:971083/1‑151 (MQ=255) gTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc > 2:217416/1‑151 (MQ=255) atgtgcaactcccACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg < 1:1156926/140‑1 (MQ=255) aaTGGCCAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg < 1:546246/151‑1 (MQ=255) agggtaaagcccACCCCCTGCCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat < 2:985657/138‑1 (MQ=255) ggccagccccACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga < 2:449502/141‑1 (MQ=255) aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg > 2:900211/1‑150 (MQ=255) ccACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg > 1:1022522/1‑151 (MQ=255) | CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG > NC_000913/2984291‑2984587 |
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: ‑ |