Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,708,973 | T→C | 12.9% | A685A (GCT→GCC) | rsxC → | SoxR [2Fe‑2S] reducing system protein RsxC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,708,973 | 0 | T | C | 12.9% | 89.3 / 4.2 | 43 | A685A (GCT→GCC) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base T (27/8); new base C (2/3); total (32/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGC > NC_000913/1708829‑1709118 | gCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGACCAGCAACCGGCTAATGAGGAGCCAGAAGACCACGTAGATCAGCGCAAAGCTGcaga > 1:1613339/1‑147 (MQ=14) cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGc > 1:1379197/1‑150 (MQ=14) gCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGACGCCGCTATTGCCCGCGCCAACGCCCGCAAACCGGAACAGCACCCGGCTAATGCGGAGCCAGAAGAACAGGTAGAACAGCGAAACGCAGCAGTAGAAGCGGCTATAGCAcg > 1:298261/1‑150 (MQ=14) ccAGAACACCAGGCCGATCCTCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAcgc < 1:2031/150‑1 (MQ=14) caacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaaa > 2:594923/1‑150 (MQ=255) gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGTTCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCAAAAGCACGAAAg > 2:966717/1‑150 (MQ=255) cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTGGACGCGGGGAATGCACGCGGCAACGCCCGCAAGCTgg > 2:410957/1‑150 (MQ=17) ccgcAACCCCGCCGGCCACGCCGCCCATGCCCGAGCCAAAGCTCGCAAACTGGAACAGCAACAGGCTAATGCGGAGCCAGAATAACAGGTTGATCCGCGCAAAGCTTCCGTCGAAGCGGCTTTTGCACGCGCCAAAGCACGCAAGCTGGa < 2:1475537/149‑1 (MQ=14) ccgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAACCAGAAGAACAGGTTGAACACCGCAAAGCTGACGTAGAAGAGGCTATTGCAAGAGCCAAAGCACGCAAGCTGGAACAGaaaa > 2:220602/1‑149 (MQ=14) cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg > 2:974813/1‑101 (MQ=255) cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg < 1:974813/101‑1 (MQ=255) cgccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACa < 1:594923/150‑1 (MQ=255) cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGACCAGCAACAGGCTaa > 2:1401240/1‑150 (MQ=255) cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTCGAACAGCAACAGGCaaa > 1:636646/1‑150 (MQ=255) cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAACAAGCAATCGCATCGTCTATTGAAAGCGACACAGCTCGAAAGCAGCAACAGCGACCGGCTaa > 2:1345346/1‑150 (MQ=18) gCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGACCAGCAACAGGCTAATGCg > 1:1218890/1‑150 (MQ=255) tGCCCTATCCAAAGCGCGCAAACGGGAACACCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAATCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCaga < 1:220602/150‑1 (MQ=255) cGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGACCAGCAACAGGCTAAAGCGGTACAcgaagac > 2:747123/1‑149 (MQ=18) cGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCACAAGCCCGCAAGCACGACCAGCAACAGGCTAAAGCCGAACCagaagac > 1:668305/1‑149 (MQ=25) cAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAATGCACGCGCCAAAGCACGCAAGAAGGAACAGCAACAGGCTAATGAGGGACCAGAAGACCAgct > 1:1607111/1‑148 (MQ=17) aGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGa > 1:103393/1‑150 (MQ=255) gcgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCCCGCGCCAAAGCACACAAACACGCACAGCACAAGGCTAAAGCGGAACAAGAAGACCAGGTAGATcc > 2:495325/1‑150 (MQ=21) cgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcg > 1:613863/1‑150 (MQ=255) cAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCACCAGGCTAAAGCGGTACCAGAAGACCAGGTAGAtacccc > 1:506444/1‑145 (MQ=18) gggAACAGCAACCGGCTAAAGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTCCCCAAGAACCGCCTATTGCCCCCGCCAAAGCACGCAAGCTCGAACAGCAACAGGCTAAAGCCGTACCAGAACAACACGATGCTCCGCGCATAg > 2:901225/1‑150 (MQ=255) ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCGAAAGCACGCAAGCTGGAACAGCAGCAGGCTAATGAGGTGCCAGAAGACCAGGTGGATCCCCCCAAAGc > 1:380395/1‑150 (MQ=17) cAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCa > 1:1547212/1‑150 (MQ=255) aaCCGGCTAATGCGTATCCAGAAGAACAGGTCGTTCCCCGCAAAGCTTCCGTCTAAGCTGCTATTGCACGCGCCAAAGCGCGCTAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGTGGCAGTTg < 1:1007113/150‑1 (MQ=14) ccGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCAAAAGCAAGCAACCACGAACAGCAACAGGCTAATGAGGTACCAGAAGAACAGGTTGATACGCGCAAAGCGGCAGTTGcc > 2:1319361/1‑150 (MQ=38) ggCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACAAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgc > 1:114302/1‑150 (MQ=255) cTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAACTGCcgctt > 2:804968/1‑148 (MQ=255) gCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCACAGCGGCAGTTGCCGCGGCTAtt > 2:1246020/1‑150 (MQ=255) gCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcg > 2:837543/1‑150 (MQ=255) agaagaACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCgcgc > 1:1488410/1‑150 (MQ=255) atcccGCGCAAAGCCGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaaa < 1:247576/147‑1 (MQ=255) tccccccAAAGCCGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgc < 1:169535/144‑1 (MQ=255) cgcgCAAAGCTGCCGTCTAAGCGGCTATTGCACGAGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccg < 1:647919/150‑1 (MQ=255) cgcgCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccg < 2:196269/150‑1 (MQ=255) cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAACAACAGGTTGATCAGCGCAAAGCGGCAGTTGCAGCGGCTATTGCACGCGCGCAGGCCaaaaaaagat > 2:247576/1‑146 (MQ=14) cgcAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccgcc > 2:461645/1‑150 (MQ=255) aaGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagc > 2:172092/1‑150 (MQ=255) aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg > 2:851840/1‑105 (MQ=11) aaGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg < 1:851840/105‑1 (MQ=255) | GCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGC > NC_000913/1708829‑1709118 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |