| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,586,056:1 | +C | 100% | intergenic (+139/+76) | PP_4061 → / ← PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,586,056 | 1 | . | C | 86.2% | 78.1 / 7.3 | 29 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
| Reads supporting (aligned to +/- strand): ref base . (3/1); new base C (12/13); total (15/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT > NC_002947/4585917‑4586197 | aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCg < 5:57061/149‑1 (MQ=39) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATg‑cg < 5:184361/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTGCATTGTAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTTAAGATCTTGGATTGATCGGGCCGGCCTCTTCGGG‑GGTAAACCCGCTCCTACGT‑GGGACC‑GTGCCGGCGTTGAATG‑CGGCGCTGCCCCTg < 6:236786/149‑1 (MQ=255) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt > 1:80456/1‑130 (MQ=38) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt < 2:80456/130‑1 (MQ=38) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAg < 5:59285/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAg > 4:225076/1‑148 (MQ=39) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGa < 2:333616/147‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCgg < 3:225076/146‑1 (MQ=38) ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTAATCGGGCCGGCCTC‑TTAGCGGGTAAAACCCCTCCTAACAGGGGC‑C‑CTGGCGGGCGTGAATT‑AGGCGCTGCCCCGGG‑TAGGAGCGGGgtt > 1:98486/1‑148 (MQ=255) ctCCCTTGAAAGGTTGAGGCTGACCGGGCGCTTGGCAAGGTGAAGATCTTGCATTGCTCGGGCAGGACTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTAc < 7:216536/147‑1 (MQ=37) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTTAAACCGCGCCTAACACGGGC‑C‑GTTGCGGGCGTGTATT‑AGGGGCGGGCCCAGGTTTGGCGCGGGTTTAcc > 3:319391/1‑149 (MQ=255) ccATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TGCGCGGGTAAACCCGCTCCTACCA‑GGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGG‑TAGGAGCGGGTTTACCcgc < 1:214946/149‑1 (MQ=14) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCgg > 1:229419/1‑100 (MQ=25) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCgg < 2:229419/100‑1 (MQ=25) ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAggcc < 5:82887/149‑1 (MQ=38) tGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTAGCGGGTAAACCCGCACCTAACA‑GGGACC‑GTGCCGGCCGTGACTG‑CGGGGCTGGCCCGGG‑GATGGTCGGGGTTTCCGGCGACGAggcgggc > 4:221871/1‑149 (MQ=255) cTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGA‑ccgt > 6:189711/1‑70 (MQ=18) ggTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt < 7:320659/78‑1 (MQ=18) ggTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt > 8:320659/1‑78 (MQ=18) ggccTC‑TTCGCGGGTAAACCCGCTACTACCAGGGGT‑A‑GCGCC‑GCATTCAAGGCCGGCACGGTCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGGGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGa > 6:351747/1‑149 (MQ=9) gcgggTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCa > 4:54105/4‑149 (MQ=255) gTAAACCCGCCCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg > 2:128520/1‑149 (MQ=255) aaaCCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc < 5:49490/140‑1 (MQ=255) aaaCCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc > 6:49490/1‑140 (MQ=255) cTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCAGCAGATGAAGACGCAATCTCCGCCACACTCATCTCACGCATCCGGAACTTCTGCAc > 3:158936/1‑149 (MQ=39) tACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 8:148490/1‑149 (MQ=255) ccA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACACGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt > 5:196913/1‑149 (MQ=255) cA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCa > 7:118619/1‑92 (MQ=255) cA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCa < 8:118619/92‑1 (MQ=255) | AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT > NC_002947/4585917‑4586197 |
| 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: ‑ |