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 | 93.1% | 87.6 / 0.4 | 29 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (13/14); total (14/15) | |||||||||||
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 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC > NC_002947/4585917‑4586199 | aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCg < 2:35505/149‑1 (MQ=39) gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCgg > 5:366781/1‑149 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCGTGGCAAGGTGATGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCtt < 5:259280/148‑1 (MQ=38) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGaa < 6:222156/149‑1 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAAt > 6:372586/1‑149 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAAt < 8:281742/149‑1 (MQ=39) gCACCGGCGTTATCATCTGATCTGCGCTGATCCATTGAAAGGTTGAGGATGACGGGGTGCTTGTCAAGGTGAATATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGTC‑C‑GTGCCTGCCTTGAATg‑cgg < 8:242595/149‑1 (MQ=25) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 6:366781/149‑1 (MQ=39) atcTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa < 1:361553/147‑1 (MQ=38) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa < 2:107204/147‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCATTGAATG‑CGGCGCTGCCCCTGGTa > 1:107204/1‑147 (MQ=38) cgcTGCTCCATTGAAAGGTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 5:66224/146‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGc < 5:372586/144‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg < 5:350049/146‑1 (MQ=38) gATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGt < 7:60157/78‑1 (MQ=18) gATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGt > 8:60157/1‑78 (MQ=18) gggTACACCCGCTCCTCCCCGGG‑ACC‑GGGCCGGCCTTGAATG‑CGGCGCTGCCCCTGTTAGGAGGGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc < 4:209127/129‑1 (MQ=38) gggTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc > 3:209127/2‑130 (MQ=255) aCCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc > 4:277474/1‑149 (MQ=255) aCCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACACGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc > 2:58424/1‑149 (MQ=255) cccGCTCCTACCAGGG‑G‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt > 2:284217/1‑149 (MQ=18) cccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAACGGGTTTACCGGCGAAGAGGCCGGCCGCGCCAGCACATAACCCCGCAGATTCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc > 2:282660/1‑148 (MQ=255) ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc < 7:99717/122‑1 (MQ=255) ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc > 8:99717/1‑122 (MQ=255) ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg < 2:158143/143‑1 (MQ=255) ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg > 1:158143/1‑143 (MQ=255) ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACtt > 7:292352/1‑149 (MQ=255) ccAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt > 5:108085/1‑149 (MQ=255) aGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGc > 1:343609/1‑149 (MQ=255) | AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC > NC_002947/4585917‑4586199 |
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: ‑ |