Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,094,581 | A→G | 13.6% | intergenic (‑309/+142) | rhaD ← / ← rhaA | rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,094,581 | 0 | A | G | 13.6% | 60.6 / 4.5 | 33 | intergenic (‑309/+142) | rhaD/rhaA | rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase |
Reads supporting (aligned to +/- strand): ref base A (13/15); new base G (3/2); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01 |
TCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGT > NC_000913/4094444‑4094722 | tcagtGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc < 1:310020/143‑1 (MQ=21) tcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc < 1:257180/143‑1 (MQ=35) aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGa > 1:94379/1‑143 (MQ=35) gCCGGATGCGACGCTGTCGCGGCTTATCAGGCCGACAACTGTGCCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCt < 2:485361/143‑1 (MQ=11) ggggCGACGCTGGCGCGGCGTATCGGGCCTACAACGGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAg < 1:26760/140‑1 (MQ=11) cGCTGTCGCGTCGGATCAGGCCTACAACGGTGGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGAt < 2:63053/143‑1 (MQ=255) ttGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatcaat < 1:371898/143‑1 (MQ=34) ccGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc > 1:458325/1‑143 (MQ=35) gCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCCCTTGTGCCGCATCCGGGAATCAATATCTTATGCGCCGCTGACGCGTcc > 1:222776/1‑142 (MQ=11) attgCCGGATGCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa > 1:288860/3‑143 (MQ=12) cacccGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc < 1:83949/140‑1 (MQ=17) aatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc < 2:333637/83‑1 (MQ=33) ccTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCCTCCTATCAGGCCCACAACTGTTGCCGCCCCg > 2:443095/1‑143 (MQ=11) cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGt < 2:369583/143‑1 (MQ=34) gCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCCACAACTGTTGCCGACCCGTCGGcc > 2:131845/1‑143 (MQ=25) cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTCGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCg < 2:168581/143‑1 (MQ=25) gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCgg > 1:104685/1‑143 (MQ=35) ccGCTGTCCCGTCTTATCAGGCCTCCAACCATAGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGa < 2:109069/142‑1 (MQ=255) cTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATaa < 1:475095/143‑1 (MQ=35) gTCGCGTCGTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAgg < 1:353019/143‑1 (MQ=21) cgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGCCTTCTCAGGCCCACCACTGTTCCCCACCCGTCGGCCGGCTGAGgcccc > 2:449203/1‑140 (MQ=11) agTCTTCCCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCt < 1:70280/142‑1 (MQ=255) tCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCGGACAACTGTTGCCGACCCGTCGGCCGGATAAGgcgcctg > 2:111599/1‑140 (MQ=11) tATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATACGGCGCTCgcgc > 1:112976/1‑143 (MQ=255) aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTCCAACTGTTGCCGACCCGTGGGCCGGATAAGGCTGTCgcgcg > 2:405319/1‑142 (MQ=11) aTCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgcc > 1:359073/1‑143 (MQ=255) tCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg < 2:419209/143‑1 (MQ=255) tCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGcgccg < 2:419177/143‑1 (MQ=255) cTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATCAGGCGCTCGCGCCGCCTCCg > 2:123053/1‑143 (MQ=255) cTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCACTCGTCAATCAATGCATGATGGGGAGCGGTTGCGTCGTATCAGGCCTACAACTGTAGCCCAGCCGCAAACCGGATAAGGGGCTCGCGCCGCCTCCg > 2:63951/1‑143 (MQ=14) cAACTGTTGCCGCGCTGTAGGCCGGATAAGGCACTCGTGCCGCATCCGGCAAGCAATGCCGGATGGGACGCTGTCGCGTCTTATCAGGCCTACAACGGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCa < 2:275873/143‑1 (MQ=14) aaCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAg > 1:166232/1‑143 (MQ=18) tATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGAGCCGCATCCCGCAgtgc > 2:57742/1‑142 (MQ=255) | TCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAGTGT > NC_000913/4094444‑4094722 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |