Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 15,410 | 0 | T | G | 41.8% | 92.6 / 11.7 | 56 | intergenic (+112/‑35) | dnaJ/insL1 | chaperone Hsp40, DnaK co‑chaperone/IS186 transposase |
Reads supporting (aligned to +/- strand): ref base T (4/28); new base G (23/0); total (28/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.40e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT > NC_000913/15261‑15536 | aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGtt < 1:759228/151‑1 (MQ=255) ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTtg < 2:180421/150‑1 (MQ=255) gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTtgtg < 2:754998/151‑1 (MQ=255) gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCT‑AAACG‑CCTACTT‑TAGGGGTGTgg > 2:171165/1‑148 (MQ=255) aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTAtt < 2:1077734/151‑1 (MQ=255) aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACg < 2:580354/151‑1 (MQ=255) agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGc < 1:490325/151‑1 (MQ=255) gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGcc < 1:663325/151‑1 (MQ=255) ttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGCTTACCCCCAT‑AAACG‑CCGACCT‑AAAGGGTGTGGTATTACGcc > 1:1115647/1‑146 (MQ=255) tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCC‑AAGCG‑CTAAATT‑AAAGGGTGTGGGATTATGCCTggt > 1:566365/1‑147 (MQ=255) gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATg < 2:794083/151‑1 (MQ=255) cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTAATCCCCT‑AAGCG‑CCCACCT‑TAGGGGTGGGGGAGTTCGACCGATATg > 2:898420/1‑149 (MQ=255) acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑CTAACTT‑AAAGGGTGTGGGATTTCGCCTGATATGAttt > 1:1118499/1‑151 (MQ=255) ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTa < 1:384032/151‑1 (MQ=255) aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCT‑AAGAG‑CTCAATT‑AAAGGGTGGGGGGTTTCGCCTGCTATGAatt > 1:182743/1‑146 (MQ=255) aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATAAGATTTAAc < 2:699596/149‑1 (MQ=255) cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCT‑AAGCG‑CGAAATT‑TAGGGGTGTGGGTTTTCGCCCGGTATTAttt > 2:990812/1‑145 (MQ=255) cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGt < 1:990839/150‑1 (MQ=255) cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca < 2:1000454/151‑2 (MQ=255) cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca < 2:1000456/151‑2 (MQ=255) cTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAACG‑CTGACTT‑TAAGGGTGTGGGGTTATCCCTGGTATGtttttaaatgcc > 2:1026181/1‑141 (MQ=255) aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAATCCCAT‑AAACG‑CTAACTT‑AAAGGGTGTGGTATTACGCCCGAAATGATTTTAAGTGcc > 2:561837/1‑147 (MQ=255) aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa < 1:373740/151‑1 (MQ=255) ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa < 1:807768/151‑1 (MQ=255) ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CGAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa < 2:654213/151‑1 (MQ=255) ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa < 1:752302/151‑1 (MQ=255) cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc < 1:840106/151‑1 (MQ=255) aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCT‑AAAAG‑CTGAAAT‑AAAGGGTGTGGTGTTAAGCCTGATATGGTTTTAAACGCCGATGaa > 1:506212/1‑143 (MQ=255) aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCC‑AAACG‑CTCAATT‑AAGGGGTGTGGGATTACGCCTGATATTTTTTAAACTGcc > 1:1072006/1‑137 (MQ=255) aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc < 2:408385/151‑1 (MQ=255) gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa < 2:266980/150‑1 (MQ=255) cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑CTCACTT‑AAGGGTTGTGGTGTTACGCCTGATATGATTTAACGTGCCGATGAATTAATCTCACGATaa > 1:837599/1‑148 (MQ=25) gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCAA‑AAAGC‑CTAACTT‑TAAGGGTGTGGGATTAAGCCTGCTATGATTTAACAAGTCGAAGAATTTCTCTCCCCATAACTg > 1:1073297/1‑150 (MQ=255) gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc < 1:614548/151‑1 (MQ=255) cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAAAG‑CTAAATT‑TAAGGGTGTGGGATTTAGCCTGGTATGATTTTAAGTGCCGATGAATTAATCTCCCGATAACTgg > 2:463900/1‑146 (MQ=255) cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGACTCCCCT‑AAAGG‑CTAACTTCAAGGGTT‑TGGTATTTCCCCCGCTTTGATTTTACGTGCCGGAGGATTAATCTCAccaaaaatagggcg > 1:14890/1‑137 (MQ=255) ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CCAACAT‑AAAGGGTGTGGGGGTATGCCCGGTTTTATTTAACCTGCCGGAGAAATAATCTCACGATAACAGGGCa > 1:617618/1‑147 (MQ=255) ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa < 1:540027/151‑1 (MQ=255) ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAACCCCCT‑AAACG‑CTAACTT‑TAAGGGTGTGGTGTTTCACCCGGTTTTATTTAATCTGCCGATGGATTACTCTCCCGCAAAATGGTCAGCa > 1:565225/1‑150 (MQ=255) gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCAT‑AAAAG‑CTCAATT‑AAGGGGTGTGGGGTTTCGCCTGATTTGATTTTAAGGGCCGCTGAATTAATTTCACGaaaaataggccacaaattt > 2:579739/1‑135 (MQ=255) gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAAAG‑CTAACTT‑AAAGGGTGTGGGGTTACGCCTGTTTTTATTTAACGTGCCGATGAAATACTCTCACCATAAATTGTCAGCAATTCt > 2:707816/1‑150 (MQ=255) ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg < 2:672461/151‑1 (MQ=255) aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCT‑AAACG‑CTAACTT‑AAGGGGTGTGGGATTACGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAAATCTGGccc > 1:1174788/1‑150 (MQ=21) aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCT‑AAAAG‑CCCACTT‑AAGGGGTGGGGGATTACGCCCGCTATGATTTAACGTGCCGGTGAATTACTCTCACCATAACTGGTCAGGAATTCTGGccc > 2:416263/1‑148 (MQ=255) aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACTCCCAT‑AAACG‑CTAAATT‑TAGGGGTGTGGGGTTACGCCTGATATGATTTAAAGTGCCGACGGATTAATTTCACAAAAACATGGCAGACATTCTTGGCCACAtt > 2:835729/1‑150 (MQ=255) cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCCCTCCCCT‑AAACG‑CCAACTT‑TAGGGGTGTGGGATTTCGCCCGTTATGATTTTAAATGCCGATGAATTACTCTCACCATAAATGGTCAGCAAATCTTGGCCCTAtt > 1:677184/1‑144 (MQ=255) tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGAGACTCCCCT‑AAGGG‑GCAACTT‑AAGGGGTGTGGGGTTACGCCTGATATGATTTTACATTGCGGTGAAATACTCTCCAGAGAAAATGTCAGCAAATCTTGGCCACATTGGT‑AAGc > 2:913716/1‑148 (MQ=255) gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCgaaga < 1:480022/151‑1 (MQ=255) atatGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCCC‑AAACGCCTACCTA‑AAGGGTT‑TGGGATTTCTCCCGGTATGATTTTAAGTTCCGATGAATTTCTCTCCCCATAAATGGTCCGCAATTCTTGCCCCCATTTGT‑AAGccc > 2:621493/1‑145 (MQ=255) atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACt < 1:368365/151‑1 (MQ=255) gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTgg < 1:407607/151‑1 (MQ=255) gCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CCAAAAT‑AAGGGGTGGGGGGTTAAGCCTGATATGATTTAACGTGCCGCTGGAATAATCTCACGATAACTGGTTCGCAATTTTGGCCCCCATTTGT‑AAAAACGAAGAAAAGGATa > 2:627366/1‑150 (MQ=9) aCCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACtt < 2:1161389/151‑1 (MQ=255) cTGTAAAGTGGCGGGGGTCACCCCCCTAAACCC‑CCAATTT‑AAGGG‑TGTGGGATTACGCCTGATATGATTTAAAGTGCCCCTGAATTAATCTTCCGATAAATGGGCAGCCAATTTGGGCCACATTTGTAAAACCCGAAAAACgggtaatt > 2:386378/1‑144 (MQ=2) gCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATg < 2:138699/151‑1 (MQ=34) aTCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg < 2:206196/151‑1 (MQ=33) cTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct > 1:1064235/1‑151 (MQ=11) | ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CTAACTT‑AAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGT‑AAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT > NC_000913/15261‑15536 |
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: ‑ |