Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,984,441 | 0 | A | C | 43.1% | 103.5 / 13.9 | 58 | G384G (GGT→GGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (30/3); new base C (0/25); total (30/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.88e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC > NC_000913/2984292‑2984585 | gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCacac > 1:393000/1‑151 (MQ=255) atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTg > 2:296365/1‑151 (MQ=255) gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCccg < 1:680446/151‑1 (MQ=255) agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCccgc < 2:611378/151‑1 (MQ=255) agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCccgc > 2:481501/1‑151 (MQ=255) aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCccgccg > 1:692239/1‑150 (MQ=255) aaaaaTGAGGGAGCCAAAAAAATATATGTCGACTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGATATTTACCGGTAAGGAGAAAATGCAGCGCTATTCGTCACTTTA‑AATGGTAATTGCCACCCCCGGG‑CCCCCGCCGATa < 1:115397/149‑1 (MQ=255) ccaaaaaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATTTTTATCGGTAGGGAGAAAAGGCAGAGCTTTTCGTCCCGTTA‑AATGGTCAAGGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCa < 2:737289/144‑1 (MQ=255) cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTcc > 2:311801/1‑151 (MQ=255) aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCt > 1:516449/1‑151 (MQ=255) aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTC‑AATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt < 2:799830/151‑1 (MQ=255) tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg > 1:828569/1‑151 (MQ=255) tatgtCGTCTTGCCACCCATCATTATTCAACAGATTAATGACGGAAGATATTTATGGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTA‑AATGGTAAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg < 2:966206/149‑1 (MQ=255) atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc > 2:69650/1‑150 (MQ=255) gACTTGCCAGCCATCATTATTCAAAAGATTAAGGACGGCAGATATTTATGGGAAGGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa < 2:717673/149‑1 (MQ=255) atttgcacgccaTCTTTTTTAAACAGATTAATGACGGAAATTTTTTATCGGTAGGAAGAAAATGCAGAGCTTTTC‑TCCCTTTCAAATGGTCAATCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCTTTTGCGGGCATTa < 2:330891/140‑1 (MQ=255) gccagccaTCATTATTCAACAAATAAATGACGGCAAATATTTATCGGTAAGGAGAAAATGCAAAGCTTTTCGTCACGTAA‑AATGGACAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa < 2:76756/151‑1 (MQ=255) cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc > 2:468793/1‑151 (MQ=255) catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt > 1:114631/1‑151 (MQ=255) atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc > 2:280543/1‑151 (MQ=255) attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 2:517579/1‑151 (MQ=255) attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 2:67928/1‑151 (MQ=255) attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg > 2:168581/1‑151 (MQ=255) attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc > 2:624095/1‑151 (MQ=255) ttCAACAAATTAATGACGGCAGATTTTTATCGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTA‑AATGGTAAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt < 2:957710/143‑1 (MQ=255) aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTA‑AATGGTCAATGCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag < 2:426380/150‑1 (MQ=255) aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga > 2:458919/1‑150 (MQ=255) aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc > 2:749975/1‑151 (MQ=255) aTTAAGGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTC‑AATGGTCAATGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa < 1:158519/150‑1 (MQ=255) tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc > 1:393593/1‑150 (MQ=255) aTGACGAGAGATTTTTATCGGGAAGGAGGAACGGCACAGTTTTTCTCCCCTTTC‑ATGGGTAAATCCCTCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg < 2:711491/151‑1 (MQ=255) gACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt > 1:638491/1‑151 (MQ=255) aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt > 2:6746/1‑151 (MQ=255) ggCAGATATTTACGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTC‑AATGGTCAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa < 2:49065/147‑1 (MQ=255) aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCAACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAGTCGGc > 1:378174/1‑150 (MQ=255) aGATAGTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTTG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCCTGAACCACAGAAACCAGGATTCg > 1:883735/1‑148 (MQ=255) gATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa > 2:676504/1‑150 (MQ=255) gATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa > 2:732016/1‑150 (MQ=255) atatTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa < 2:326469/150‑1 (MQ=255) ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTA‑AATGGTCAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg < 2:393174/150‑1 (MQ=255) tATGGGTAAGGAGAAAATGCAGACCTATTCGTCACGTTC‑AATGGCCAGTCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa < 2:955884/150‑1 (MQ=255) tATCGGAAAGAAGGAAAGGCAGAGTTTTTCGTCACGTTA‑AATGTTCATTGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa < 1:295239/150‑1 (MQ=255) ggTAAGGAGGAAAGGAAGAGCTTTTCTTCCAGTTA‑AAGGTTAATTGCCCCCCCCCGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg < 2:639393/150‑1 (MQ=255) aaGGAGAAAATGCAGAGCTTTTCGTCACGTTA‑AATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc < 2:969715/149‑1 (MQ=255) ggagaaAATGCAGAGCTTTTCGTAACGTTC‑AAGGGCCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg < 1:927814/143‑1 (MQ=255) gagaaAATGCAGAGCTTTTCGTCACGTTA‑AATGGTAAATGCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg < 2:875665/151‑1 (MQ=255) aggaAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg > 2:498751/1‑151 (MQ=255) ggaAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg > 1:891672/1‑151 (MQ=255) aagtccagactgtttCGTCACTTTC‑AAT‑GTAAACCCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg < 1:904298/138‑1 (MQ=255) aaaTGCAGATCTTTTCGTCCCGTTA‑AATGGTCAACCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt < 2:148068/150‑1 (MQ=255) aaTGCAGATCTTTTCGTCAGGTTC‑AAGGGTCAGTGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg < 1:839526/146‑1 (MQ=255) agCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg > 1:1003057/1‑151 (MQ=255) actttttCGCCAGGTTC‑AATGTTCAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt < 1:511930/145‑1 (MQ=255) ttCGTCACGTTA‑AAGGGTCAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt < 2:998343/150‑1 (MQ=255) aCGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACCAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCCTGAACCAGAGAAACCAGGAGTCGGCAACCGGAAGCGCCTATCGGGCGACCGAGTGCGATCGCGcc > 1:693040/1‑150 (MQ=255) cGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac > 2:526148/1‑151 (MQ=255) c‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt > 1:128943/1‑151 (MQ=255) gggTCACTACAACACCCCGG‑CCCCCGCCGATACAAAGAGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAACCCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCTAGTGCGATCGCGCCACcattgac < 1:964875/150‑1 (MQ=255) tGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc > 2:902469/1‑151 (MQ=255) | GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTC‑AATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC > NC_000913/2984292‑2984585 |
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: ‑ |