Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,032,846 | Δ1 bp | 10.8% | intergenic (+54/+82) | dppA‑IV → / ← PP_5284 | dipeptide transport protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,032,846 | 0 | C | . | 10.8% | 108.0 / 8.8 | 37 | intergenic (+54/+82) | dppA‑IV/PP_5284 | dipeptide transport protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (19/14); new base . (2/2); total (21/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.25e-01 |
TGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCA > NC_002947/6032713‑6032995 | tgctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTAc < 1:692191/150‑1 (MQ=255) atcattggtgatggggagccgttgtcgatggtcgagcatgggtggattgtggtctgagggagtactgagagccttatagggtcagtacgggcctcTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTAc > 1:29873‑M1/96‑150 (MQ=255) ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑t‑t > 2:463653/1‑135 (MQ=255) ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑t‑t < 1:463653/135‑1 (MQ=255) cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAg > 2:789130/1‑150 (MQ=255) cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAg > 2:129183/1‑150 (MQ=255) gTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCgg < 1:513136/150‑1 (MQ=255) gatggtcgagcatgggtggattgtggtctgagggagtactgagagccttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGc > 2:344898‑M1/70‑121 (MQ=255) gatggtcgagcatgggtggattgtggtctgagggagtactgagagccttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGc < 1:344898‑M1/52‑1 (MQ=255) ggTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa < 1:129183/150‑1 (MQ=255) aTTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAgg < 1:328483/150‑1 (MQ=255) acggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTccc > 2:33502/1‑150 (MQ=255) ggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTAccc > 1:621951/1‑105 (MQ=255) ggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTAccc < 2:621951/105‑1 (MQ=255) tgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCACGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAg > 1:3771/1‑150 (MQ=255) ttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGc > 2:462736‑M1/23‑74 (MQ=255) ttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGc < 1:462736‑M1/52‑1 (MQ=255) ttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTAc > 1:439065‑M1/23‑77 (MQ=255) ttatagggtcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTAc < 2:439065‑M1/55‑1 (MQ=255) tcagtacgggccctTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAAC‑TGCCAG‑GCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAc < 2:505194‑M1/136‑1 (MQ=255) ccTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCag < 2:547126/150‑2 (MQ=255) cGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCGTCAGACCACACTCCACCCCTGCt > 1:99443/1‑150 (MQ=255) gCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTgg > 2:521005/1‑47 (MQ=255) gCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTgg < 1:521005/47‑1 (MQ=255) ccTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGa < 2:692812/100‑1 (MQ=255) ccTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGa > 1:692812/1‑100 (MQ=255) tctTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGAcc > 2:784561/1‑150 (MQ=255) TCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAAC‑TGCCAGGC‑TGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATc > 1:505194‑M1/1‑150 (MQ=255) ACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCa < 1:101096‑M1/112‑1 (MQ=255) ACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCa > 2:101096‑M1/1‑112 (MQ=255) CCGCTCCTACAGGGGTAGCGTCAAC‑TGCCAG‑GCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAAt < 2:782019‑M1/119‑1 (MQ=255) CCGCTCCTACAGGGGTAGCGTCAAC‑TGCCAG‑GCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAAt > 1:782019‑M1/1‑119 (MQ=255) ctGGGTTGCGCAAGCCTGGCTG‑T‑TGAGGCTTGCCCGTTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCaa > 2:443707‑M1/3‑150 (MQ=255) GGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACcaatag < 1:443707‑M1/150‑3 (MQ=255) GCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat > 1:789130‑M1/1‑150 (MQ=255) GCGCCAGCCTCGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat > 1:33502‑M1/1‑150 (MQ=255) ccAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGAtgtg > 2:151162/1‑150 (MQ=255) cTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCa < 2:3771/150‑1 (MQ=255) | TGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAG‑T‑TGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATCCA > NC_002947/6032713‑6032995 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |