Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,900,744 | G→A | 6.9% | G418D (GGC→GAC) | PP_1703 → | assimilatory nitrate reductase/sulfite reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,900,744 | 0 | G | A | 6.9% | 171.0 / 2.8 | 58 | G418D (GGC→GAC) | PP_1703 | assimilatory nitrate reductase/sulfite reductase |
Reads supporting (aligned to +/- strand): ref base G (21/33); new base A (2/2); total (23/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01 |
CAACTGCCACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGA > NC_002947/1900596‑1900881 | cAACTGCCACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACTGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGc < 1:305869/150‑1 (MQ=255) ccACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGTCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑Tgcgc > 1:1311650/1‑150 (MQ=255) ccACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑Tgcgc > 1:335545/1‑150 (MQ=255) aCCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCaa < 1:790829/150‑1 (MQ=255) cccGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAg < 2:1116810/150‑1 (MQ=255) cGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGc < 2:1350973/150‑1 (MQ=255) gaGCCGGCGAAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCCTTGCGC‑AGGCCgt < 2:1369865/149‑2 (MQ=255) aGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGaa < 2:593929/150‑1 (MQ=255) gCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAAc < 2:284286/150‑1 (MQ=255) gcgcAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGGTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGt < 1:800360/150‑1 (MQ=255) aGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGaa < 1:440580/141‑1 (MQ=255) aGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCAGaa > 2:440580/1‑141 (MQ=255) gCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTa < 2:382569/150‑1 (MQ=255) gAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGcc < 2:930182/150‑1 (MQ=255) gAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGcc < 2:929991/150‑1 (MQ=255) aCGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTc < 1:1377916/150‑1 (MQ=255) aCGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTc < 2:421547/150‑1 (MQ=255) aCGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTc < 1:188503/150‑1 (MQ=255) gcgcTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACAcc > 2:285781/1‑150 (MQ=255) gTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGt > 1:993265/1‑117 (MQ=255) gTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGt < 2:993265/117‑1 (MQ=255) gTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaac < 1:999544/150‑1 (MQ=255) ttCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAg < 1:658464/105‑1 (MQ=255) ttCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAg > 2:658464/1‑105 (MQ=255) gggCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaac < 2:543413/150‑1 (MQ=255) gCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCg > 2:804545/1‑150 (MQ=255) aGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGACC‑Tgcg > 1:73035/1‑87 (MQ=255) aGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGACC‑Tgcg > 1:73076/1‑87 (MQ=255) aGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGACC‑Tgcg < 2:73076/87‑1 (MQ=255) aGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGACC‑Tgcg < 2:73035/87‑1 (MQ=255) ggTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATc < 1:285781/150‑1 (MQ=255) tgaagtgaaAGACTACTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCct > 2:420106/1‑150 (MQ=255) aagtgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCctgc > 1:391374/1‑150 (MQ=255) aagtgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTAGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCctgc > 1:502321/1‑150 (MQ=255) tgaaAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCtgctgc < 1:194597/150‑1 (MQ=255) gTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTgg < 1:778824/150‑1 (MQ=255) aTCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaacc > 1:1297393/1‑112 (MQ=255) aTCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaacc < 2:1297393/112‑1 (MQ=255) tCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCg > 1:687692/1‑150 (MQ=255) cagcaACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAgg < 1:420106/150‑1 (MQ=255) aCCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAgg < 2:687692/150‑1 (MQ=255) ccGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGgt > 1:177692/1‑150 (MQ=255) ggCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGtgat < 2:1311650/150‑1 (MQ=255) cAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGt < 2:1326989/54‑1 (MQ=255) cAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGt > 1:1326989/1‑54 (MQ=255) gtgtCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCa > 1:315574/1‑150 (MQ=255) cGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaacc < 1:271851/86‑1 (MQ=255) cGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAaaccaacc > 2:271851/1‑86 (MQ=255) cAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagcga > 2:427584/1‑150 (MQ=255) cAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagcga < 2:177692/150‑1 (MQ=255) aCCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGTCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCaa < 2:1309718/142‑1 (MQ=255) aCCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGTCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCaa > 1:1309718/1‑142 (MQ=255) tCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCtgtg > 2:803524/1‑111 (MQ=255) tCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCtgtg < 1:803524/111‑1 (MQ=255) gCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCt > 1:423525/1‑150 (MQ=255) aaGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGa < 2:315574/150‑1 (MQ=255) ttcGGCCTTGCGC‑AGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagc > 1:987327/3‑129 (MQ=255) ttcGGCCTTGCGC‑AGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACagc < 2:987327/127‑1 (MQ=255) | CAACTGCCACCCGGCAGCCTGCGCAGCGAGGGGGGCGAGGGTACGGTGGCGCTGTTCGGGCAGATGAAAGCCGGTGAAGTGAAAGCCTGCTGGGTCATCTGCAGCAACCCGGTGGCCAGTGTCGCCAACCGTCAGCAAGTGATCGACGGCC‑TGCGCAAGGCCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGA > NC_002947/1900596‑1900881 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |