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 | 96.3% | 77.7 / ‑2.8 | 27 | intergenic (‑304/‑177) | yieP/BW25113_RS19520 | FadR/GntR family transcriptional regulator/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (11/15); total (12/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01 |
GGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACG > NZ_CP009273/3934915‑3935072 | ggggTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt < 1:151829/90‑1 (MQ=255) tctcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGaa < 1:146958/90‑1 (MQ=255) tGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagc < 2:171052/90‑1 (MQ=255) gagaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagca > 1:116206/1‑90 (MQ=255) gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa > 2:372074/1‑90 (MQ=255) gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa > 2:309758/1‑90 (MQ=255) ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt < 1:190400/72‑1 (MQ=255) ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCt > 2:190400/1‑72 (MQ=255) gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt < 2:52962/90‑1 (MQ=255) agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc < 1:189044/90‑1 (MQ=255) agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc < 1:396573/90‑1 (MQ=255) agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc < 2:450127/90‑1 (MQ=255) aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTa < 2:434996/90‑1 (MQ=255) cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAc > 2:364991/1‑90 (MQ=255) cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAc < 2:438692/90‑1 (MQ=255) aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa < 1:434547/90‑1 (MQ=255) aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa < 2:506330/90‑1 (MQ=255) aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa < 2:288931/90‑1 (MQ=255) tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa > 2:173073/1‑90 (MQ=255) ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt > 2:417593/1‑90 (MQ=255) ctGTAGCGGGAAAGCGTATTATGCACA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtg > 1:303789/1‑90 (MQ=39) gCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTggg > 2:103769/1‑90 (MQ=255) gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt > 1:479689/1‑54 (MQ=14) gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt < 2:479689/54‑1 (MQ=14) gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa > 1:323376/1‑90 (MQ=255) gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa > 2:183883/1‑90 (MQ=255) tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg < 2:40361/90‑1 (MQ=255) | GGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACA‑CCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACG > NZ_CP009273/3934915‑3935072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |