Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,718,633 | G→A | 8.9% | A426T (GCC→ACC) | PP_5020 → | methyl‑accepting chemotaxis protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,718,633 | 0 | G | A | 8.9% | 120.9 / 3.2 | 45 | A426T (GCC→ACC) | PP_5020 | methyl‑accepting chemotaxis protein |
Reads supporting (aligned to +/- strand): ref base G (20/21); new base A (2/2); total (22/23) | |||||||||||
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.09e-01 |
TCGGCGGTGACAGAACAGACCAGCGCCGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAGCGACAAGATCGGCAGCGTA > NC_002947/5718489‑5718770 | tCGGCGGTGACAGAACAGACCAGCGCGGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGAc > 2:147757/1‑150 (MQ=255) gACAGAACAGACCAGCGCCGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGc > 2:487818/1‑150 (MQ=255) agaCCAGCGCCGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAgg < 2:170169/150‑1 (MQ=255) aCCAGCGCCGGGGTCAACAATCAGAAGGTCGAAACCCTGTTTGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGt < 1:262234/150‑1 (MQ=39) aGCGCCGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGatgacatt < 2:443863/150‑8 (MQ=255) ggTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCgtggt > 1:364233/1‑150 (MQ=255) tCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGtggtgg > 1:562295/1‑150 (MQ=255) tCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGtggtgg < 2:627900/150‑1 (MQ=255) aacaaTCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGc > 1:635183/1‑150 (MQ=255) caaTCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGa < 2:756548/150‑1 (MQ=255) gAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGt < 1:504665/150‑1 (MQ=255) gTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGca > 2:75430/1‑150 (MQ=255) gAAACCGACCAGGTGGCGACTGCGCTGATCCATATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAg > 1:451615/1‑150 (MQ=255) aaaCCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGc < 2:449797/97‑1 (MQ=255) aaaCCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGc > 1:449797/1‑97 (MQ=255) aaaCCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGa < 1:75430/150‑1 (MQ=255) cGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGa < 1:220456/150‑1 (MQ=255) aGGTGGAGACGGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGACGACCTGATGGCCGACCAGCAGGCCCGCGAGGGTCATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGcc < 1:487818/150‑1 (MQ=39) gTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTg > 1:732525/1‑150 (MQ=255) tGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc < 2:21909/137‑1 (MQ=255) tGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc > 1:21909/1‑137 (MQ=255) aGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGACCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACt > 2:456964/1‑150 (MQ=255) aGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGACCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACt > 2:456962/1‑150 (MQ=255) aCTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGACCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCagc < 1:456964/150‑1 (MQ=255) aCTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGACCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCagc < 1:456962/150‑1 (MQ=255) gCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCagcgaa < 2:451615/150‑1 (MQ=255) ccACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc > 1:740794/1‑117 (MQ=255) ccACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc < 2:740794/117‑1 (MQ=255) ccACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCAAGCCAAGGGGTCAATGCCagcgaat < 2:635183/150‑2 (MQ=255) aCGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCagcgaagc < 2:819797/142‑1 (MQ=255) aCGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCagcgaagc > 1:819797/1‑142 (MQ=255) aCGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGACCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAAc > 2:288931/1‑131 (MQ=255) aCGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGACCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAAc < 1:288931/131‑1 (MQ=255) gCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCCGGGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAg > 1:7820/1‑150 (MQ=255) cACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc > 2:128807/1‑99 (MQ=255) cACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAgc < 1:128807/99‑1 (MQ=255) ctctaTGCCTAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGAcc < 2:732525/146‑1 (MQ=255) tGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAg < 1:165240/150‑1 (MQ=255) gCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAGCGACAAGATc > 1:291114/1‑150 (MQ=255) ttCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAGCGACAAGATCgg > 1:197608/1‑150 (MQ=255) aGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTgg > 1:510089/1‑80 (MQ=255) aGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTgg < 2:510089/80‑1 (MQ=255) gccgccTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAGCGACAAGATCGGCAGCGTa > 1:421622/1‑150 (MQ=255) ccgccTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCt > 1:748821/1‑115 (MQ=255) ccgccTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCt < 2:748821/115‑1 (MQ=255) | TCGGCGGTGACAGAACAGACCAGCGCCGGGGTCAACAATCAGAAGGTCGAAACCGACCAGGTGGCGACTGCGATGAACCAGATGACTGCCACCGTGCACGAAGTGGCACGCAATGCCGAACAGGCTTCGGAAGCCGCCTTGATGGCCGACCAGCAGGCCCGCGAGGGTGATCGCGTGGTGGGCGAGGCAGTGGCACAGATCGAGCGCCTGGCCAGCGAAGTGGTCAACTCCAGCGAAGCGATGAACCTGCTCAAGACCGAAAGCGACAAGATCGGCAGCGTA > NC_002947/5718489‑5718770 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |