Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,903,525 | G→C | 12.3% | intergenic (+50/‑221) | PP_1703 → / → nirB | assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,903,525 | 0 | G | C | 12.3% | 126.0 / 9.7 | 49 | intergenic (+50/‑221) | PP_1703/nirB | assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit |
Reads supporting (aligned to +/- strand): ref base G (25/16); new base C (3/3); total (28/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01 |
CGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGG > NC_002947/1903385‑1903666 | cGCCGAGCATGCCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCATGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgc < 7:365723/149‑1 (MQ=255) gCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAgt < 4:215512/149‑1 (MQ=255) ccGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAg > 2:188790/1‑147 (MQ=255) ccGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaat > 2:290099/1‑137 (MQ=255) ccGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaat < 1:290099/141‑5 (MQ=255) gatgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTgggcgggc > 8:196014/1‑125 (MQ=255) gatgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTgggcgggc < 7:196014/125‑1 (MQ=255) tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaatgc > 8:140271/1‑128 (MQ=255) tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaatgc < 7:140271/134‑7 (MQ=255) ggcggcTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCcgc > 1:94783/1‑149 (MQ=255) cggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTAAACCcgctcc < 2:1700/149‑4 (MQ=255) ctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAagag < 3:190294/149‑1 (MQ=255) ctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAagag > 5:227257/1‑149 (MQ=255) ggcgggtTCTGTCGACGGCTTGAGCACGGCCAGGCGTTATC‑GGCGCGATGTGTACTGCAGTGCGCCAAGCTGGGGTGTTGGGCACCTTGAATGGGCATTGGCGGGGCTGGCCCAGCGCAAGGGTAGGAGCGGGTTTACCCGCGAAGAgg < 6:327159/148‑1 (MQ=255) gcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAagagtc > 8:83464/1‑147 (MQ=255) cggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAgg < 2:94783/147‑1 (MQ=255) cTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAATGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccag > 1:134579/1‑149 (MQ=255) ttGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccagg > 3:361467/1‑149 (MQ=255) ggCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCAc < 1:188790/149‑1 (MQ=255) cAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGc > 1:378694/1‑148 (MQ=255) aGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCaa > 2:372960/1‑149 (MQ=255) aTC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaaccc > 4:362311/1‑127 (MQ=255) aTC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaaccc < 3:362311/127‑1 (MQ=255) tgtgTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTcc > 1:177360/1‑149 (MQ=255) aCTGAAGTGCTCCAAGATGGGGTGTTGTGCACCGTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCGGGCACAGGCAACCGAGGACACCCAGGCATCCCGATTGGCACGCTCCCTGCt > 8:147381/1‑149 (MQ=38) tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTNTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt > 1:376962/1‑148 (MQ=255) tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt > 4:231265/1‑148 (MQ=255) tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt > 1:371030/1‑148 (MQ=255) tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt > 8:122459/1‑148 (MQ=255) tGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTAACCGCGAAGAGGCCGACACAGGCAACCCAGGAAACACACGCAACCCGAATGGCACGCTCCCTGCtt > 8:23730/1‑148 (MQ=255) aaGTGCTNCAAGATNGGGTGTTGTGCACCTTGNATNTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt < 2:376962/146‑1 (MQ=255) aaGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCtt < 2:371030/146‑1 (MQ=255) aaGTGCTCCAAGATGGGGCGTTGTGCACCTTGAATGTGCATTCGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGGAACCCGATTGGCACGCTCCCTGCTTaa < 4:94552/148‑1 (MQ=255) cTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCgg < 4:361467/149‑1 (MQ=255) tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTACCCTCGGa > 4:74915/1‑149 (MQ=255) tgaaGATGGGGTGTTGTGCACCTTGAATGTGCATTGGGTGGGCTGGCCCAGCGCAAGTGTAGGAGCGTGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCGACTCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGa < 4:324174/146‑1 (MQ=39) ccAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGa > 3:324174/1‑148 (MQ=255) gATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATac < 6:227257/147‑1 (MQ=255) gtgtTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGc > 5:382885/1‑149 (MQ=255) gttgtGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaa > 2:65636/1‑95 (MQ=255) cACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGa < 2:177360/148‑1 (MQ=255) aTGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATac < 7:142279/125‑1 (MQ=255) aTGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATac > 8:142279/1‑125 (MQ=255) gtgCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCGCTGATCAACGGCGAt > 4:118349/1‑149 (MQ=255) gCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATc > 6:140296/1‑148 (MQ=255) ttgcgCTGGGCCCGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGCCACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAACGCCCACTGATCAACGGCGATCacc > 6:124295/5‑147 (MQ=38) aTTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAAc > 8:276251/1‑149 (MQ=255) ttGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAAcc > 7:310179/1‑149 (MQ=255) gCTGGGCGGGGCGAGCGCAAGTGTAGGAGCGGGATTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACCACAAAGCCCACTGATCAACGGCGATCAACCggg < 8:310179/149‑1 (MQ=255) gCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCggg < 2:134579/149‑1 (MQ=255) gCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCggg < 2:378694/149‑1 (MQ=255) | CGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATC‑GGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGG > NC_002947/1903385‑1903666 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |