Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,071,939 | C→G | 10.8% | L214F (TTG→TTC) | pstB‑II ← | phosphate ABC transporter ATP‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,071,939 | 0 | C | G | 10.8% | 106.7 / 5.5 | 37 | L214F (TTG→TTC) | pstB‑II | phosphate ABC transporter ATP‑binding protein |
Reads supporting (aligned to +/- strand): ref base C (14/19); new base G (2/2); total (16/21) | |||||||||||
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.43e-01 |
AGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACC > NC_002947/6071793‑6072083 | aGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAg > 1:321818/1‑150 (MQ=255) aGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAg < 2:463208/150‑1 (MQ=255) tcggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCgag < 2:62979/150‑1 (MQ=255) cggtatcCCCGAATTCGACCAGTTTGCCCATGTAAATCGACTCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCgaag < 1:29576/150‑3 (MQ=255) atcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg < 1:57909/150‑1 (MQ=255) aTTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGGCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAgcg < 2:159555/150‑1 (MQ=255) aCCAGTTTGCCCATGTACATGAACGCGGCGCAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGTGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAg > 1:438442/1‑150 (MQ=255) cAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGAGACGTTATGACCAGGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCa < 2:438442/150‑1 (MQ=255) ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCgaagaa > 2:224706/1‑121 (MQ=255) ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCgaagaa < 1:224706/125‑5 (MQ=255) aCATGAACGCGGTGTAAACCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTCGATTTCAATTCATAGATCAGTTCTTCGACCGTCAAGGTCGAGATCGAGTCCATCGCCCCGCACTGTTCGTCAagcg > 1:249165/1‑149 (MQ=255) tGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCg > 1:121416/1‑117 (MQ=255) tGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCg < 2:121416/117‑1 (MQ=255) aCGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAgcagcagca > 1:75552/1‑150 (MQ=255) ggTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTc > 2:261413/1‑150 (MQ=255) tAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGc < 1:377545/150‑1 (MQ=255) tAATCCGAAACACGCGCAGCTTGCTGCATGTTGTGGGTCACCATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTCCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGc < 2:386082/150‑1 (MQ=255) tCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCAcc > 1:205421/1‑141 (MQ=255) tCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCAcc < 2:205421/141‑1 (MQ=255) cGAAACACGCGCCGCTTGCTGCATGGTGTGGGTCACGATGAAGATGGGGTACTTGGATTTCAATTCATAGATCAGTGCTTCGACCTTCAAGGTCGAGATCGGGTCCAGTGCCGAGCACGGGTCGTCAAGCAGCAGCACCTCCGGCTCCAc < 1:205442/150‑1 (MQ=255) tgtgGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGAc < 1:261413/150‑1 (MQ=255) gggTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgct < 2:75552/150‑1 (MQ=255) tCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgctgct > 1:481878/1‑150 (MQ=255) gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg > 2:160074/1‑150 (MQ=255) gtgtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGAc < 1:160074/150‑1 (MQ=255) aCTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGc < 1:121924/150‑1 (MQ=255) cTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc > 2:300344/1‑150 (MQ=255) ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgc > 2:151064/1‑126 (MQ=255) ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgc < 1:151064/126‑1 (MQ=255) ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc > 2:382884/1‑149 (MQ=255) ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGcc < 1:382884/149‑1 (MQ=255) gATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAg < 2:321818/150‑1 (MQ=255) gATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAg < 1:300344/150‑1 (MQ=255) ctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt > 2:369041/2‑150 (MQ=255) cGACATTCAACGTCGAGATCGGGTCCAGCGCCAAGCACTGTTCGCCAAGCGGCAGCACCTCCGGCTCCACGGCAATTGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCTGACAGGCACAGCGCCGACTCGGGCAGGCGGTCCTTGa < 2:249165/150‑1 (MQ=255) gACCTTCAAGGTCGAGATCGGATCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGAGAGGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGAc > 2:95622/1‑150 (MQ=255) aCCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGAcc > 1:382847/1‑150 (MQ=255) | AGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACC > NC_002947/6071793‑6072083 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |