Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,190,245 | C→G | 10.6% | G142R (GGG→CGG) | PP_5612 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,190,245 | 0 | C | G | 10.6% | 125.1 / 5.7 | 46 | G142R (GGG→CGG) | PP_5612 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (18/23); new base G (2/3); total (20/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTCCGCCTTCTGCCACTCTGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAGGTACGAGTGG > NC_002947/4190102‑4190386 | cTCCGCCTTCTGCCACTCTGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCg > 6:50970/1‑149 (MQ=255) aCTCTGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGAt < 7:235230/138‑1 (MQ=255) aCTCTGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGAt > 8:235230/1‑138 (MQ=255) tctGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCagaag > 7:51751/1‑148 (MQ=255) tGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGaaa > 1:88004/1‑149 (MQ=255) tGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCaa < 8:51751/148‑1 (MQ=255) aaCCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCGAATGGGATAccc > 3:37338/1‑129 (MQ=255) aaCCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCGAATGGGATAccc > 3:41374/1‑129 (MQ=255) aaCCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCGAATGGGATAccc < 4:41374/131‑3 (MQ=255) aaCCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCGAATGGGATAccc < 4:37338/131‑3 (MQ=255) ctaaTCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCTATTCTCGGTTAAGGCAAAGGAAACGCTCGCGCCACGGCTGCTGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGa < 1:215839/146‑1 (MQ=255) aTCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGaaa < 2:115715/148‑1 (MQ=255) aTCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGaaa < 3:219242/148‑1 (MQ=255) cAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCt < 1:275674/148‑1 (MQ=255) gATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGa > 7:223368/1‑136 (MQ=255) gATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGa < 8:223368/136‑1 (MQ=255) aataCGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGAt < 7:248698/108‑1 (MQ=255) aataCGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGAt > 8:248698/1‑108 (MQ=255) tctcGTTGCCGGATCGAGTCTCGGTTAATGCACAGGAAACGGTCGCGCCACGGCTGGGGGTTTTGTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAAc < 6:153641/125‑1 (MQ=255) tctcGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACACAGAAGAAACCAAGTTCGAAAc > 5:153641/1‑125 (MQ=255) ccGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATa > 3:72052/1‑148 (MQ=255) ccGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATAc < 4:72052/149‑1 (MQ=255) ccGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACATGATCAGCACTGACCATTTCCTGATTAATa > 7:63542/1‑148 (MQ=255) gAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTc < 2:88004/149‑1 (MQ=255) aGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTcc < 5:50970/149‑1 (MQ=255) tctcGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGaa > 4:33090/1‑105 (MQ=255) tctcGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGaa < 3:33090/105‑1 (MQ=255) aGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAg < 7:282524/149‑1 (MQ=255) ggCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAAc < 8:24735/97‑1 (MQ=255) ggCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAAc > 7:24735/1‑97 (MQ=255) cacaGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACa > 6:112789/1‑128 (MQ=255) cacaGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACa < 5:112789/128‑1 (MQ=255) cacaGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAggtg < 3:30554/149‑1 (MQ=255) acaGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAggtgg > 8:249014/1‑149 (MQ=255) aGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTc < 5:252888/149‑1 (MQ=255) gCTGCGTGATTTCTTCACATCGCACAAATGGTAGCCCTTTCGGATTCGACTCAGTATAAACCAAGTTCGTAACCTGATCAGCACTGAGCATTTGCTGATAAATAGAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAAttttt < 8:63542/149‑1 (MQ=255) tGCGGGATTTCTTCACATCGCACAAAGGGTATCCCAATCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGATCCTCGTAGTCAAGGTGCTCGAGCCACTTTAGAAtttttt > 4:247583/1‑148 (MQ=255) gggATTTCTTCACATCGCACAAAGGTTATCCGATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTTATTATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGc < 5:140500/149‑1 (MQ=255) ttcACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACttt > 2:78004/1‑127 (MQ=255) ttcACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACttt < 1:78004/127‑1 (MQ=255) cGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAgg > 5:54658/1‑149 (MQ=255) cacaAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGtt < 8:131140/89‑1 (MQ=255) cacaAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGtt > 7:131140/1‑89 (MQ=255) gggTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAGGTACGAgtg < 7:249014/149‑1 (MQ=255) tATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAGGTACGAgtgg > 8:118212/1‑147 (MQ=255) tATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAGGTACGAgtgg < 7:118212/147‑1 (MQ=255) | CTCCGCCTTCTGCCACTCTGAGATGAACCTCATCCAGTCATGATAATACGCGACCTCTCGTTCCCGGATCGAGTCTCGGTTAAGGCACAGGAAACGCTCGCGCCACGGCTGCGGGATTTCTTCACATCGCACAAAGGGTATCCCATTCGGATACGACTCAGAAGAAACCAAGTTCGAAACCTGATCAGCACTGAGCATTTGCTGATAAATACAGCGTTCCTCGTAGTCAAGGTGGTCGAGCCACTTTAGAATTTTTTTAGCCCGACGATCCCCAGGTACGAGTGG > NC_002947/4190102‑4190386 |
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: ‑ |