Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,434,787 | G→A | 26.7% | intergenic (‑563/+171) | ynaE ← / ← ttcC | cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,434,787 | 0 | G | A | 26.7% | 21.2 / 6.4 | 15 | intergenic (‑563/+171) | ynaE/ttcC | cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related |
Reads supporting (aligned to +/- strand): ref base G (5/6); new base A (2/2); total (7/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 |
AAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATG > NC_000913/1434643‑1434929 | aagaaAAATATCAATGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTTTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTg < 2:198773/150‑1 (MQ=255) agaaAAATATCAAGGTGATATAGCCTATATGCCTGTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTgg < 2:231675/150‑1 (MQ=255) gaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACATTTTGGa > 2:160335/1‑150 (MQ=255) gaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACATTTTGGa > 2:34441/1‑150 (MQ=255) tGATATAGCGTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTgg < 2:231136/150‑1 (MQ=255) tGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTgg < 1:184788/150‑1 (MQ=255) taGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaaa > 2:123791/1‑150 (MQ=255) gATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACATTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCaaaaaaa < 1:160335/150‑1 (MQ=255) gATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACATTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCaaaaaaa < 1:34441/150‑1 (MQ=255) ttGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTa < 1:19492/150‑1 (MQ=255) aaaaaCCGCAGTCACGTCGTATGCAAGAACGTGCAGCCGTTGGAATGACGTTGATTCAGACGTGTTCGAATTTTTATCGGAGGCAAAAAAGGTGCAAAACGCTGCAAAACGGGCAAAAAAGGGGCAAAAAAAGAGAGGATTAACGTAGCt > 2:185644/1‑150 (MQ=255) cgtcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtg > 2:238821/1‑150 (MQ=255) gAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAac > 2:19449/1‑150 (MQ=255) ggttggttTGACTTTGATTGAGACGTTTTGGGATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGaaaaa < 2:247052/150‑1 (MQ=255) ttgaGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATg > 2:219027/1‑150 (MQ=255) | AAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATG > NC_000913/1434643‑1434929 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |