Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,740,474 | C→A | 20.6% | intergenic (+39/+32) | hutF → / ← PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,740,474 | 0 | C | A | 20.6% | 70.1 / 14.1 | 34 | intergenic (+39/+32) | hutF/PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Reads supporting (aligned to +/- strand): ref base C (10/17); new base A (4/3); total (14/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01 |
GCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGG > NC_002947/5740329‑5740622 | gcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgc > 2:264660/1‑149 (MQ=255) aTGTGATGGTGAATGGGCAAGGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa < 7:149962/149‑1 (MQ=255) gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCGGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc < 8:32335/149‑2 (MQ=38) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa > 6:141752/1‑139 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa < 5:141752/139‑1 (MQ=255) aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac < 3:118301/148‑1 (MQ=255) aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac < 2:175761/148‑1 (MQ=255) gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAAcc < 1:264660/149‑1 (MQ=255) gCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct > 2:342683/1‑136 (MQ=38) gCAATGGGCGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct < 1:342683/138‑3 (MQ=38) cAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGTTTCAGCCCATAGCCGGGAACCCCCCGCAGGGGTGGaacccaaccc > 2:103198/1‑149 (MQ=255) aTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCCTGGAAACCCCCGCACGGGTTCAACCCAAcccc > 4:122100/1‑148 (MQ=255) aTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAAcccc > 6:72948/1‑148 (MQ=255) gtggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaa > 5:367018/1‑148 (MQ=255) gcgcCTGGGGCGGAATGCCCAGTAAGATGAACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCa < 6:4135/149‑1 (MQ=255) ggcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGAcgcccgcc < 3:122100/149‑1 (MQ=255) aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg < 6:367018/149‑1 (MQ=255) tgtTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtgt < 5:72948/149‑1 (MQ=255) gggCTGCGCGTACCGGGGTTGCGTTGCAGCCCTGCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc < 1:350514/148‑1 (MQ=255) gCTGAGCGTACCGGGGTTGCCTTTGAGCCCATAGCGGGCCAGCCCCCCCACGGTTGGAACCCAACCCCCAAAAACACAGTGAAATTTTT‑‑ATATTTTCTGATAGGACGCCCCCCACATCAAACGGCGGGGGTAAAAACCCAGCTGCcgcg > 4:278985/1‑149 (MQ=255) ggAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc < 7:261022/148‑1 (MQ=255) aGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAACCTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGccc < 7:238151/148‑1 (MQ=255) cGGGGTTGGGTTGCAGCCCTTCGGGGTCAAGCCCGCGTAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTCTCTTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca < 1:113920/149‑1 (MQ=255) cGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc < 3:278985/149‑1 (MQ=255) gggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag < 1:328196/149‑1 (MQ=255) gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc < 4:243008/125‑1 (MQ=255) gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc > 3:243008/1‑125 (MQ=255) gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTAGTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc > 3:192431/1‑149 (MQ=255) tgCCCTTCGCGTGCAAGCCCGCGTAGGGTTGCAACGCTACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTg < 2:52790/70‑1 (MQ=255) aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc > 2:363359/3‑149 (MQ=38) cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg > 8:168061/1‑149 (MQ=39) ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGa < 1:329400/86‑1 (MQ=37) ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGa > 2:329400/1‑86 (MQ=37) aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGAGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg > 5:320634/1‑149 (MQ=255) ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCgg > 7:160291/1‑149 (MQ=255) | GCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTA‑‑CTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGG > NC_002947/5740329‑5740622 |
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: ‑ |