Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,443,277 | 0 | A | C | 47.4% | 52.7 / 38.2 | 57 | V274G (GTA→GGA) | secY | preprotein translocase membrane subunit |
Reads supporting (aligned to +/- strand): ref base A (24/6); new base C (0/27); total (24/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.49e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATA > NC_000913/3443133‑3443424 | gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCg > 1:960307/1‑151 (MQ=255) aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGcc < 2:404844/151‑1 (MQ=255) gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCata > 2:996779/1‑151 (MQ=255) gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAACGAAGATTGCCGGGATT‑ACCCCCGCCata > 2:274920/1‑151 (MQ=255) tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCa > 1:1008692/1‑151 (MQ=255) ccccAGGCTGAAAATACAGCAAATTT‑G‑TTGCCCGCAAGTTCCAACCAGTACCCCCCCAGACCCAGG‑CCGCGGTGGCCGCCGGAAAAAAAATAATACTGGAAGCAAATATTGCCGGGATC‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtg < 2:340687/150‑1 (MQ=255) cccgggcagaaaaaccAGCAAATTT‑G‑TGGCCACCCAGTTCCACCCAGTCCCCCCCCCGAACCATG‑CGGCGAGGGTCCCGGGAAACAAAATAATACGGAAAGAGAAGATGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtg < 2:758450/135‑1 (MQ=255) ccTGAAAATCCAGCAAATTT‑T‑TTGTCACCCAGTCCCAACCATCCCCCCCCCCGAACCAGG‑CCGGGAGGGCCGCGGGGAACAAAAAAATCCTGGAGGCGAAGATTGCGGGTTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc < 2:698155/149‑1 (MQ=255) cTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct > 1:777462/1‑150 (MQ=255) cTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct > 2:692355/1‑150 (MQ=255) cTGAAAATACAGAGAAATT‑T‑TTGCCAGCCATTCCCAACCAGTCCCGCCCCCGACCCAGG‑CCGGGATGGCCCCGGGGAAAAGAATAATTCGGGAAGCTAAGATGGCGGGGTTT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:579705/151‑1 (MQ=255) gCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct > 2:845241/1‑148 (MQ=255) caaaaTG‑G‑TTGCCACCCAGTTCCAACCGGCCCCCCCCCCAAACCATG‑CCGGGAGGGTCGCGGGAAACAAAAAAATACGGGAAGCAAATATTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 2:671220/147‑1 (MQ=255) cGAAATT‑G‑TTGCCAGCCATTCCCAACAAGTCCCGCCCCCGACCCATG‑ACGCGATGGCCGCCGGGAACAAAATAATACTGAAAGCAAAGTTTGCGGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc < 2:995545/147‑1 (MQ=255) agAATTT‑T‑TTGCCACCCAGTCCCAACAAGACCCGCCCCCAACCCAGG‑AGGGGAGGGTCGCGGGGAAAAAAATAATACTGGAAGC‑AAAATTTCCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 2:923436/148‑1 (MQ=255) agAAGTT‑G‑TGGCCGGCCATTCCCACCCAGTCCCCCCCCCAAACAATG‑AGGCGTTGGCCGCCGGAAACAAAAAAATAGGGGAAGGGAAGTTTGCCGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 2:716025/148‑1 (MQ=255) aaaattt‑g‑ttgCCAGCCATTCCCACCCAGTCCCCCCCCCGACCCATG‑CCGGGATGGGCCCCGGGAACAAAATAAGACTGGAAGGGAAGATTGCGGGGTTT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 2:255496/144‑1 (MQ=255) aaaaaTT‑G‑TGCTCACCCAGTTCCAACCAGTCCCCCCCCCGACCCAGG‑CCGCGGTGGTCCCCGGGAACAAAATAATACGGGAGGAAAAATTGGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc < 1:48577/145‑1 (MQ=255) gAAATTGG‑TTGTCAGCCGGTTCAAACAAACCCCCCCCCGGAACCATG‑ACGCGAGGGCCGCGGGAAAAAAAAAAATTGTGGAAGAGAAGATGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 1:652666/149‑1 (MQ=255) ccaaattttttgcccgCCATTCCCACCCAGTCCCCCCCCAGACCCATG‑ACGGGAGGTCCGCGGGAAACAAAATAATCCTGGAAGCAAAGATTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 1:301738/135‑1 (MQ=255) aaaaTT‑GTTTGCCAGCAAGTCCCAACCAG‑ACCCCCCCCGGACCCGGAACGGGATGTGCCCCGGAAACAAAATAATCCGGAAGGCGAAGTTGGCCGGTATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 2:227234/148‑1 (MQ=255) aaaattttttgcCACCCAGTCCCAACCAGCACCGCCCCCAACCAATG‑ACGCGAGGGTCCCGGGAAAAAAAATAATTGGGGAAGAAAAGTTTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt < 1:421179/141‑1 (MQ=255) aaaTT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa > 2:534808/1‑150 (MQ=255) aaaTT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa > 2:162264/1‑150 (MQ=255) aaTT‑G‑TTGTCAGCCAGTCCCACCCAGTCCCCCCCCCGACCCAGG‑ACGGGTGGGTCGCGGGAAAAAAAATAATTCGGGAAGAGAAGTTGGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGaca < 1:605100/151‑1 (MQ=255) g‑ttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga > 1:366650/1‑150 (MQ=255) ttgtCAGCCAGTCCAAACCAGCCCCCCCCCAGACCCATG‑GCGCGTGGGTCGCCGGGAACAAAATAATAGGGGAAGAGAATTGGGCCGGTATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg < 1:512802/151‑1 (MQ=255) tggccAGCCAGTCCCACACAGCCCCGCCCCAGAACCAGG‑CCGGGTTGGTCGCGGAAAACAAAAAAATCCGGGAGGCGAATTTTGCGGGTATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac < 1:801651/146‑1 (MQ=255) tcccgCCAGTCCCACCCAACCCCGCCCCCAAACCAGG‑CCGGGGGGGTCGCCGGAACAAAAAAAATCCTGGAGGCGAAGTTTCCGGGAATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac < 2:517699/147‑1 (MQ=255) tCAGCCAGTTCCACCCAGCCCCGCCCCCGACCCATG‑CCGCGATGGCCCCCGGGAAAAAAAAAATACTGGAAGCAAAGATTGCCGGAATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc < 1:612049/151‑1 (MQ=255) gcggCCGGTTCCACCAGGTCCCCCCCCCAAACAATG‑CCGCAAGGGCCGCCGGAACCAAAATAATACTGAAAGGGAAGTTGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac < 2:97456/147‑1 (MQ=255) cAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg > 1:789257/1‑149 (MQ=255) gTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa > 1:252639/1‑151 (MQ=255) ccAACCAGTCCCCCCCCCAAACAATG‑AGGGGGGGGCCCCGGGAAACAAAATAATCCTGGAAGGGAAGATTGCGGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt < 1:938313/151‑1 (MQ=255) cAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt > 1:956765/1‑151 (MQ=255) aCCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTACTGACGttt > 1:690070/1‑150 (MQ=255) ccccccGACCCAGG‑CCGGGGTGCCCGCCGGAAAAAAAATAATATGGGAAGGGAAGTTTCCCGGAATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt < 2:228149/150‑1 (MQ=255) cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa > 1:50463/1‑151 (MQ=255) cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa > 2:441161/1‑151 (MQ=255) cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa > 2:547852/1‑151 (MQ=255) caaaCCAGG‑ACGCGAGGGTCGCCGGAAAAAAAAAAATACTGGAAGCAAAGATTGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac < 2:868452/149‑1 (MQ=255) cATG‑AGGGGATGTGCGCGGGAAACAAAATAATACGGAAAGAGAATTTTGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc < 1:534574/151‑1 (MQ=255) aCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGTTTGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc < 1:363199/150‑1 (MQ=255) tgctcccggggAACAAAATAATACTGAACGAAAATTTGCCCGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc < 1:171315/142‑1 (MQ=255) aacaaaaaaatcggggAAGGGAAATTTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa < 2:218429/137‑1 (MQ=255) gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa > 2:95915/1‑150 (MQ=255) gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa > 2:631381/1‑150 (MQ=255) gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa > 1:307731/1‑150 (MQ=255) gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCCCAACAAATa > 2:172404/1‑150 (MQ=255) ataCTGGAGGGAAAGAGTGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa < 2:794381/151‑1 (MQ=255) tggagggaaaGATGGCGGGATTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa < 2:899456/143‑1 (MQ=255) gCGAAGTTGGCCGGTTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc < 1:928711/150‑1 (MQ=255) aGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg > 1:718043/1‑150 (MQ=255) ttgcccGGTATC‑ACCCCCGCCATATTCNNTTNCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc < 1:800127/145‑1 (MQ=255) aTTGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc < 2:22840/149‑1 (MQ=255) gCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata > 1:883645/1‑151 (MQ=255) ggATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaa > 1:65131/1‑151 (MQ=255) aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaata > 2:619893/1‑148 (MQ=255) aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaata > 1:979357/1‑148 (MQ=255) | GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATA > NC_000913/3443133‑3443424 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |