Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 4,266,090 | C→A | 16.5% | S231* (TCA→TAA) | yjcC → | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,266,090 | 0 | C | A | 16.5% | 72.6 / 2.1 | 25 | S231* (TCA→TAA) | yjcC | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
Reads supporting (aligned to +/- strand): ref base C (9/11); new base A (2/2); total (12/13) | |||||||||||
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.83e-01 |
TGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCG > CP009273/4265859‑4266310 | tgGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTAGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTAATATCAACGTCTTATAACCCATTTTTATAATCATCTTAt > 2:207378/1‑271 (MQ=255) tGGAGTGTAGTGATGTCTGATGACCCGACATTGTAATGGGGTGTTGATGATTCGGTGACGAAAACCTTTTTCTCGTTAAACAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACGCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTagg < 1:466297/275‑1 (MQ=255) tGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGCCAACACGCCCAATTGCAGCCATTGTTGCGACTTAATATCAACGTCATATAACCCATTTTTATAAACATCATATTTTTGAGTTGCCCGACGTTAtttccgg > 2:230590/1‑271 (MQ=255) tGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTAGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTAATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg < 1:207384/275‑1 (MQ=255) aCATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCt > 1:181373/1‑275 (MQ=255) tGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTGATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAttt > 1:253429/1‑240 (MQ=255) cAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCAATTTTTATAATCATCTTATTTTTGCGTTGCCAGCCGGTATTTTGGGGAGTCTTGTTCTTCTTTTATTNTGGCttcttag > 2:47337/1‑269 (MQ=255) tattatACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGTGGCCTCGGCATCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCAAATTGCAGCCATTGTTGCGACTTCATATCAAc < 2:240866/172‑1 (MQ=255) aCGAAAACCTTTTTCTCGTGAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTCACCGTACAAAGAAATGTCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCt < 1:317566/245‑1 (MQ=255) aGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAtt > 1:390291/1‑152 (MQ=255) aTGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTActc < 2:390278/152‑1 (MQ=255) tGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTtata > 1:193208/1‑79 (MQ=255) tATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAAccc < 2:193202/79‑1 (MQ=255) aTTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCg > 2:221792/1‑160 (MQ=255) tttATATGCGACAGTTTTTTCGACAATACGCCCAATTGCAGCCGTTGTTGCGACTTAATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGTGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAg < 1:230597/275‑1 (MQ=255) tgtttgcacncgTTTATTCCCCCAAACCCCACTTTGCAGCCATTTTTGCGATTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTGTCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCg < 2:181367/264‑1 (MQ=255) atatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaa < 1:221799/160‑1 (MQ=255) aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATc > 1:287010/1‑52 (MQ=255) cAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCaa > 1:28143/1‑203 (MQ=255) tttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATCTCAACg < 2:287002/52‑1 (MQ=255) ttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACttt < 2:28142/203‑1 (MQ=255) tCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt > 1:147198/1‑121 (MQ=255) cAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTAct < 2:147192/121‑1 (MQ=255) aaaaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCaa > 1:43499/1‑169 (MQ=255) cGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAAcgcg < 2:43497/169‑1 (MQ=255) | TGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCG > CP009273/4265859‑4266310 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |