Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,443,277 | 0 | A | C | 42.4% | 104.2 / 23.6 | 59 | V274G (GTA→GGA) | secY | preprotein translocase membrane subunit |
Reads supporting (aligned to +/- strand): ref base A (33/1); new base C (0/25); total (33/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.59e-16 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT > NC_000913/3443133‑3443423 | gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCg < 1:927025/151‑1 (MQ=255) agTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCatat > 2:148127/1‑151 (MQ=255) tAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTc > 2:283486/1‑151 (MQ=255) acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc > 1:252048/1‑151 (MQ=255) cacaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACt > 1:326957/1‑151 (MQ=255) aagggggTTGCCCGGGCTGAAAATACAGCAAATTTTTTGTCAGCCAGTCCCAACCAGTCCCCCCCCCAAACCAT‑GACGCGATGGTCGCGGGAAACAAAATAATCCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg < 2:137407/146‑1 (MQ=255) aggCGTTTGCCCAGCCTAAAAATACAGAAAATTTGTTGCCCGCCAGTCCAAACCAGCCCCCCCCCCG‑ACCATGGCCGGGGTGGTCGCGGGGAAAAAAATAATACTGGAAGCGAATTTGGCGGGGATCCCCCCCGCCATATTCACTTTCAGc < 1:580659/149‑1 (MQ=255) cGGTTCCCCGGGCTCCAAAAACAGAAAAATTGTGGTCGCCCAGTCCCACCCAGCCCCCCCCCCGAACCAT‑GCCGCGATGGTCGCCGGAAAAAAAAAAATACTGGAAGCGAAGATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCg < 1:837108/149‑1 (MQ=255) ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt > 2:614021/1‑150 (MQ=255) gTTCCCCAGGCTGAAAATAGAGAGAAATTTTTCCCAGCCAGTCCCACCCATCCCCCCCCCCGACCCAG‑GCCGCGTGGGCCGCCGGAAAAAAAATAATACTGGAAGCGAAGATTGCGGGAATCCCCCCCGCCATATTCACTTTCAGCGGt < 1:285444/149‑1 (MQ=255) cccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgt > 2:880206/1‑151 (MQ=255) cccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgt > 2:461442/1‑151 (MQ=255) cccAGGCTGAAAATACAGCAAATTTTTTGTCACCCATTCCCAACCAGTCCCCCCCCCAACCCAG‑GCCGGGATGGCCGCGGGAAAAAAAATAATCCTGGAAGCAAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAAt < 2:231040/149‑1 (MQ=255) ccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg > 2:89098/1‑151 (MQ=255) tGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct > 1:983035/1‑151 (MQ=255) gCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg > 2:297259/1‑151 (MQ=255) aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg > 2:853879/1‑151 (MQ=255) aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg > 2:65768/1‑151 (MQ=255) aatccagaaaaaTTTTTGCCCGCCAGTTCCAACCAGTACCCCCCCCGAACCAT‑GACGTGATGGTCCCCGGAAAAAAAAAAATATGGGAGGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc < 2:318974/142‑1 (MQ=255) aaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc > 1:21563/1‑151 (MQ=255) aGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt > 2:452366/1‑151 (MQ=255) gCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt > 2:922401/1‑150 (MQ=255) cGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt > 2:618981/1‑149 (MQ=255) aatttgttgTCACCCAGTCCCAACCAGTCCCGCCCCCGAACCAG‑GACGCGATGGTCGCGGGGAACAAAAAAATACTGGAAGGGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac < 2:776640/148‑1 (MQ=255) aaaTTGTTGCCAGCCTATTCCAACCAGTCCCCCCCCCGACCCAG‑GACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAATTGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac < 2:217328/151‑1 (MQ=255) aaaTTGTTGCCACCCATTCCCAACCAGTCCCCCCCCCGAACCAT‑GACGCGATGGCCGCGGGAAACAAAAAAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgca < 2:524958/143‑1 (MQ=255) gttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac > 2:389621/1‑151 (MQ=255) tgccacccattcccAACCAGACCCCCCCCCAAACCAG‑GCCGGGGGGGCCGCGGGGAAAAAAAAAATACGGGAAGAGAAGATTGCCGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac < 1:157240/137‑1 (MQ=255) aGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt > 1:11593/1‑151 (MQ=255) aGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTA‑CCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg > 1:614717/1‑151 (MQ=255) caaGTTCCAACCAGTCCCCCCCCCGACCCAT‑GCCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg < 1:247917/148‑1 (MQ=255) ccAGTCCCCCCCCCGAACCAG‑GCCGCGATGGTCGCCGGGACCAAAATAATACTGAAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc < 2:247613/150‑1 (MQ=255) cAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc > 1:705523/1‑151 (MQ=255) cccGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc > 1:648766/1‑151 (MQ=255) ccagggcgggATGGTGGCGGGAAACAGAATAATACTGGAAGCGAATGTGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac < 2:886603/138‑1 (MQ=255) t‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcgg > 2:742849/1‑151 (MQ=255) aCGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg > 2:950535/1‑151 (MQ=255) cgggATGGCCGCCGGAAAAAAAAAAATACGGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc < 2:926915/146‑1 (MQ=255) gggtggtccgcgGGAAACAAAAAAATCCGGAAGGCGAAGATGGCGGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt < 1:586508/143‑1 (MQ=255) gcggtggccgcggaaaaCAAAAAAATACTAGAAGCGAATATTGGGGTTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt < 1:851581/136‑1 (MQ=255) gcgAGGGCCGCCGGAAACAAAAAAATACTGGAAGGAAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACTACAATgcggcg < 2:921365/149‑1 (MQ=255) gagggTCGCCGGAAAAAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa < 2:31819/148‑1 (MQ=255) cGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc > 1:299291/1‑151 (MQ=255) ggaaaCAAAAAAATACTGGAAGCGAAGATTCCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGa < 2:707359/138‑1 (MQ=255) aGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt > 1:708250/1‑150 (MQ=255) aGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt > 2:859660/1‑150 (MQ=255) aGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt > 1:887890/1‑150 (MQ=255) aataatCCTGAAAGGGAAGTTTGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa < 2:433782/149‑1 (MQ=255) aataTTGTAGGCGAAGATTGCATGAGTACCCCCCGCCATATTCTCTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca < 2:498536/148‑1 (MQ=255) acgggaagggaagATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAg < 1:801503/148‑1 (MQ=255) ggAGGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa < 2:299307/150‑1 (MQ=255) gaaaGGGAAGATTCCCGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACaaa < 2:992319/145‑1 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg > 2:873244/1‑151 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg > 1:704349/1‑151 (MQ=255) gAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt > 1:530438/1‑151 (MQ=255) aTTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa > 1:241418/1‑150 (MQ=255) gCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaat > 2:628803/1‑150 (MQ=255) gggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact < 2:527431/150‑1 (MQ=255) gATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat > 1:496106/1‑151 (MQ=255) | GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCAT‑GACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT > NC_000913/3443133‑3443423 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |