Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,934,995 | +C | intergenic (‑308/‑173) | yieP ← / → BW25113_RS19520 | FadR/GntR family transcriptional regulator/16S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,934,991 | 1 | . | C | 93.2% | 123.2 / 3.3 | 44 | intergenic (‑304/‑177) | yieP/BW25113_RS19520 | FadR/GntR family transcriptional regulator/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base . (0/3); new base C (24/17); total (24/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.61e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT > NZ_CP009273/3934912‑3935079 | aGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGcgcc > 1:1438301/1‑90 (MQ=255) gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTgag > 1:1175514/1‑90 (MQ=255) aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGTGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAaaagcaaagc > 1:263560/1‑90 (MQ=11) aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAaaagcaaagc < 1:513220/90‑1 (MQ=255) aCTCCGGCAGAGAAAGCAAAAATAAATGCNTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGACAAaaagcaaagc < 2:546192/90‑1 (MQ=14) cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCg > 2:362498/1‑90 (MQ=255) cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCg > 1:769/1‑90 (MQ=255) ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc > 2:1116623/1‑90 (MQ=255) ggCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt < 1:1519299/90‑1 (MQ=39) agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc < 1:1083780/90‑1 (MQ=255) gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGctct > 1:742302/1‑90 (MQ=255) aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt > 1:1117749/1‑90 (MQ=255) aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt > 1:622408/1‑90 (MQ=255) gCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTaa > 1:229772/1‑90 (MQ=255) aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa < 2:1405707/90‑1 (MQ=255) aaaaTAAATGCTAGACTCTGTAGCGGGAAATGGTATTATGC‑A‑CACCCCCGTGCCGCTGAGAAAAAGCAAAGCGGCGCTGCTCTTTAACaa < 2:804465/90‑1 (MQ=11) aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc < 1:407886/90‑1 (MQ=255) aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc < 1:1314090/90‑1 (MQ=255) tGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg < 1:1229032/90‑1 (MQ=255) tGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg < 2:224869/90‑1 (MQ=255) ttGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACa > 2:537319/1‑90 (MQ=255) tGACTCTGTAGCGGGAAAGGGTATTATGC‑ACCC‑CCCCGCGCCGCCGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa < 2:1320846/90‑1 (MQ=14) tGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa > 1:599780/1‑90 (MQ=255) ctGTAGCGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt > 2:115736/1‑90 (MQ=255) cGGGAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc < 2:1422293/90‑1 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg > 1:138793/1‑74 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg < 2:138793/74‑1 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc < 2:1213696/80‑1 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc > 1:1213696/1‑80 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa > 2:1419379/1‑90 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa > 2:802260/1‑90 (MQ=255) gggAAAGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGCGTGGGCa > 1:839649/1‑90 (MQ=11) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc > 1:931688/1‑90 (MQ=255) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc > 1:1315625/1‑90 (MQ=255) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc > 2:1558873/1‑90 (MQ=255) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc > 1:654303/1‑90 (MQ=255) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc > 2:357853/1‑90 (MQ=255) aaaGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc < 1:537319/90‑1 (MQ=255) aGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGaaaaa < 2:292856/36‑1 (MQ=11) aGCGTATTATGC‑A‑CACCCCCGCGCCGCTGAGaaaaa > 1:292856/1‑36 (MQ=11) gCGTATTATGCAA‑CC‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGac < 1:293550/90‑2 (MQ=11) tatGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg < 1:1419379/90‑1 (MQ=255) tGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa < 1:237035/90‑1 (MQ=255) tGC‑A‑CACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa < 2:1035197/90‑1 (MQ=255) ca‑cCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCGTTAACAATTTATCAGACAATCTATGTGGGCACTCGAAGATACGGATTCtt > 2:1421758/1‑90 (MQ=255) | AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGC‑A‑CA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT > NZ_CP009273/3934912‑3935079 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |