| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 15,395 | 0 | G | A | 38.9% | 43.3 / 12.5 | 36 | intergenic (+97/‑50) | dnaJ/insL1 | chaperone Hsp40, DnaK co‑chaperone/IS186 transposase |
| Reads supporting (aligned to +/- strand): ref base G (5/17); new base A (14/0); total (19/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.24e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT > NC_000913/15257‑15536 | tttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTTACTCCCCAAAGCGATAAAt‑t > 2:884208/1‑148 (MQ=255) tGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCATAAACGCTAACT‑TAAggg > 1:187238/1‑148 (MQ=255) gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTtgtg < 2:915674/151‑1 (MQ=255) aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GGCACTCCCATCAAAGCTGACT‑TACAGGGTGTGGTAtt > 2:244890/1‑151 (MQ=255) gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAAAGCTAAAT‑TAAAGGGTTTGGTATTACGcc > 1:1003860/1‑151 (MQ=255) tttcGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GATACACCCACAAGCGCTAACA‑TTAGGGGTggggggttataccttc > 1:591308/1‑139 (MQ=255) gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATg < 2:639880/151‑1 (MQ=255) acaacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTGCGCCTGATATGa < 2:563098/148‑1 (MQ=255) cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAACGCCAACT‑TAAAGGGTGTGGGATTTTGCCTGat > 1:1000195/1‑146 (MQ=255) gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACACCCCCAAAAGCTAACT‑TTAGGGGTGTGGGGTTACGCCCGGTattatt > 1:116122/1‑151 (MQ=255) cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg < 2:711582/151‑1 (MQ=255) tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAACGCTAACT‑TTAGGGGTGTGGGGTTACGCCTGGTATGATTTTAAGTg > 1:980011/1‑146 (MQ=255) ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAATCA‑TCCCATTAACCCTAACT‑TAAGGGGTGTGGGAATAAGCCCGATATGATTTaaattggccacga > 1:208676/1‑139 (MQ=255) ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAAAGCTAACT‑TTAGGGGTGGGGTATTATGCCTGATATGATTTAAAGTGc > 2:716746/1‑144 (MQ=255) cTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCAATCCCCTAAACGCCAACT‑TAAGGGGTGTGGTATTACGCCCGATATTATTTTAAGTGCCGATGaa > 1:936193/1‑150 (MQ=255) cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc < 2:1005733/151‑1 (MQ=255) aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct < 2:32142/151‑1 (MQ=255) cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAACGCTACCT‑TCAGGGGTGTGGGTTTATGCCTGATATTATTTTAAGGGCCGTCGAAATACTCTCAc > 2:865156/1‑148 (MQ=255) cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑AGCAATCCCCAAAGAGAGAACTATAAGGG‑TGTGGAGTTATTGATGATATGATTTAACGCGCCgccggaataat > 2:601127/1‑133 (MQ=255) gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa < 1:284882/151‑1 (MQ=255) cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑AACAATCCCATAAACGCTAACC‑TAAAGGGTGGGGGATTTAGCCCGATATGATTTAACGTGCCCCAGAAATAATTTCACCATAAATGGTCAGCAAtt > 2:325452/1‑150 (MQ=255) gTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc < 2:540877/151‑1 (MQ=255) tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCAAAACGCCAAAT‑TAAGGGGTGGGGGGTTTAGACCGATATGATTTAAAGGGCCCGGGAATTACTCTCCCGATAAATGGTCAACAAtttttggc > 1:173178/1‑145 (MQ=255) tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑GTCACTCCCCTAAACGCTCACT‑TTAGGGGTGTGGGATTTCGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAAATCTTGCCCCTAtt > 1:424324/1‑151 (MQ=18) tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATGAACCCCCCAAAGGCTAACT‑TAAGGGGGGGGGGATTAAGCCTCCTATGCTTTAAAGTGTCGATGAAGTTCTATCACGATAAATTGTCAGAAACTTTGGCCCCTAtt > 2:992424/1‑151 (MQ=255) aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg < 2:57135/151‑1 (MQ=255) gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg < 1:112177/151‑1 (MQ=255) gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa < 1:717026/150‑1 (MQ=255) gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAg < 1:858374/151‑1 (MQ=255) atatGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagag < 2:106098/151‑2 (MQ=255) aaGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc > 1:205657/1‑150 (MQ=38) aaGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg < 1:568613/151‑1 (MQ=38) gTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa < 1:245058/151‑1 (MQ=37) gTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa < 1:884533/151‑1 (MQ=37) gCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg < 1:563745/151‑1 (MQ=34) cTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct > 1:629991/1‑151 (MQ=11) | TTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGG‑ATCACTCCCATAAGCGCTAACT‑TAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT > NC_000913/15257‑15536 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |