| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 2,956,956 | C→A | 18.3% | T40T (ACG→ACT) | ppdB ← | putative prepilin peptidase‑dependent protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 2,956,956 | 0 | C | A | 18.3% | 57.6 / 2.8 | 26 | T40T (ACG→ACT) | ppdB | putative prepilin peptidase‑dependent protein |
| Reads supporting (aligned to +/- strand): ref base C (9/11); new base A (3/2); total (13/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.45e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCA > CP009273/2956724‑2957190 | gCACATGCTCCTTCAGACGCAATCCAATCTGGTCGGACTCTTTTACCGGTTGGCGATCCCAGATACCGTTACTGATCGCCCCCCACCGCACAATGACACAGTCACCCTGTCCGACAATTTACAGCCATTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCGGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCGCGCTGTCACGCAGGCAGAAAGCGTgccgcc < 1:176608/275‑1 (MQ=255) gCTCCTACAGACGAAATCCAATCTGTTCGTACTCTTTTACCGGTTCGCGATCCCATATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTCCGCCGTTACAGATGCCATGACAATACCCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTAGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGc < 1:429060/275‑1 (MQ=255) ttCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACACGCAGCCAGATTTCATCTTCCAGCGCCCGCTTACGGGTGCTAGTTAAACTTTAACGCTGTCACGCAGGCAGAAAGCGTGCCACCCccaccacaaa > 2:429045/1‑271 (MQ=255) aCGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACCCTGTAACGCAGGCAGAAAGCGGGCCGCCCCCACCNACAATAcac > 2:193072/1‑273 (MQ=255) aCGAAACCCAAGCTGGTCGTCCTCTGTTCCCGCTTCGCGATCCCAGATACCGTTACTGTTCGCAGCACCCTTCACAATGACACAGTCACCCTGTTCGACAAGTTCCAGCCCTGTGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGGGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGc < 2:436881/275‑1 (MQ=255) ggTTCGCGCTCCCAGTTGCCGGTACTGATCGCATCCCACTGCACAAGGACACAGTCCCCCTGTCCGACAATTTCCAGCACCTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCAATAGCAATCAACACTTCcag < 1:193079/275‑1 (MQ=255) aTACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTGGTTAAACTTTCACCCTGTACTGCAGGCAGAAAGCGTGCCGCCCACAGCAACAATACGCTACTGATCCCCATAGCAATCACAACATCAAGCAGACAACAACCtag > 1:316310/1‑273 (MQ=255) tGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGcaac > 2:86430/1‑191 (MQ=255) acaGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATa < 1:86431/191‑1 (MQ=255) tGTCCGCCAATTTCCAGCCCTTCTCCGGTACAGATGCCATTACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGGAGCCAGAGTTCGTCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAATACTTCCAGCCGAGAAAAACCTTGCTCTTTTAc < 1:442661/239‑1 (MQ=255) tCCGACAATTTCCAGACCTCCGCCTGTACCGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTAGTTaaa > 1:370119/1‑130 (MQ=255) aCAATTTCCAGCCCTCCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACttt < 2:370107/130‑1 (MQ=255) tACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGt > 2:59970/1‑134 (MQ=255) gATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTgccg < 1:59970/134‑1 (MQ=255) cAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGt < 2:304970/86‑1 (MQ=255) gCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCa > 2:249879/1‑150 (MQ=255) ccTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAAt > 2:68796/1‑127 (MQ=255) tctGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCa < 1:249888/150‑1 (MQ=255) tGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGc < 1:68796/127‑1 (MQ=255) aTGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGCGTCTTCCCCACGGCGAAACCACCAGCCACCACCCGCCCGCTGAGGTTTTGACGCGAACATg > 2:424139/1‑275 (MQ=255) cTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTAGTTAAACttt > 2:94906/1‑71 (MQ=255) tCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAg > 1:317488/1‑92 (MQ=255) gcgACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTAGTTAAACTTTCACg < 1:94907/71‑1 (MQ=255) gACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTg < 2:317476/92‑1 (MQ=255) cAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGCAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGCCAAATCCGCAGTCTTCCCCCCGGCGAAACCAACAGCCACCACCCGCCCGTTGCGTTTTTCAAGCTCATATGCCCCGCCCATGCCTTCTTgcgcc > 2:148741/1‑273 (MQ=255) cTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCtt > 1:275981/1‑112 (MQ=255) | GCACATGCTCCTTCAGACGAAATCCAATCTGGTCGGACTCTTTTACCGGTTCGCGATCCCAGATACCGTTACTGTTCGCATCCCACTGCACAATGACACAGTCACCCTGTCCGACAATTTCCAGCCCTTCGCCGGTACAGATGCCATGACAATAACCCGCCCTCTGGAGATGCTTCGCGACGGTAAATACCCGCAGCCAGATTTCATCTTCCAGCGCCAGCTTACGGGTGCTCGTTAAACTTTCACGCTGTAACGCAGGCAGAAAGCGTGCCGCCCCCAGCAACAATACGCTACTGATCGCCATAGCAATCAACACTTCCAGCAGAGAAAAACCTTGCTCTTTTACAGGCATCCTTCTGTTTCTCCTTGCTGACAAAGCCGGAGTCTTCCCCACGGCGAAACCACCAGCCACCACTCGCCCGTTGAGTTTTTGAAGCGAATATGCCCGGCCCATGCGGTATTGCGCA > CP009273/2956724‑2957190 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |