Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,242,004 | C→A | 18.2% | E157* (GAG→TAG) | bfr‑II ← | bacterioferritin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,242,004 | 0 | C | A | 18.2% | 78.1 / 9.9 | 33 | E157* (GAG→TAG) | bfr‑II | bacterioferritin |
Reads supporting (aligned to +/- strand): ref base C (14/13); new base A (3/3); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01 |
AGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGT > NC_002947/1241886‑1242145 | aGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg > 1:288978/1‑141 (MQ=255) aGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg < 2:288978/141‑1 (MQ=255) gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa < 1:334369/150‑1 (MQ=255) ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCt > 1:333526/1‑150 (MQ=255) aTCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc < 2:470850/150‑1 (MQ=255) aTCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc > 1:718690/1‑150 (MQ=255) aGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATc > 1:68529/1‑150 (MQ=255) aCGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc < 1:672050/150‑1 (MQ=255) ttCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg > 2:595298/1‑98 (MQ=255) ttCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAgg < 1:595298/98‑1 (MQ=255) tGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGttt > 1:620327/1‑131 (MQ=255) tGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGttt < 2:620327/131‑1 (MQ=255) gATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGctct > 1:117552/1‑150 (MQ=255) ttCCGGCAAGCATAAAAAAACCCGCGAGCAGCGCATGTTTGGATGGGGCGAATAGTTACTCTTCGCCCATCTGCGAATGCAGTTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGGTTCCAGCCAGTCGATGTGCTCTTCCt > 1:496017/1‑150 (MQ=255) gCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct > 1:690344/1‑85 (MQ=255) gCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct < 2:690344/85‑1 (MQ=255) gCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGAc > 2:552931/1‑150 (MQ=255) gCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACACTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGAc > 2:552730/1‑150 (MQ=255) gCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGa < 2:588960/150‑1 (MQ=255) tAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGa > 2:513930/1‑150 (MQ=255) aaaaaaaCCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGat > 1:78009/1‑150 (MQ=255) aCCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc < 2:78009/150‑1 (MQ=255) ttttttATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGc < 1:513930/150‑1 (MQ=255) tATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAg < 2:333526/150‑1 (MQ=255) ggggCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCg > 2:300842/1‑108 (MQ=255) ggggCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCg < 1:300842/108‑1 (MQ=255) gCGAAGAGTAACTATTCGCCCCTCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc < 2:555557/66‑1 (MQ=25) gCGAAGAGTAACTATTCGCCCCTCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc > 1:555557/1‑66 (MQ=25) gCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCg > 2:27636/1‑105 (MQ=255) gCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCg < 1:27636/105‑1 (MQ=255) gCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc < 2:581035/122‑1 (MQ=255) gCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc > 1:581035/1‑122 (MQ=255) tAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGt < 2:68529/150‑1 (MQ=255) | AGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCGATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGT > NC_002947/1241886‑1242145 |
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: ‑ |