Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,478,883 | G→A | 16.1% | intergenic (+93/+352) | PP_2172 → / ← PP_2173 | hypothetical protein/AraC family transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,478,883 | 0 | G | A | 16.1% | 64.6 / 6.3 | 31 | intergenic (+93/+352) | PP_2172/PP_2173 | hypothetical protein/AraC family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base G (16/10); new base A (2/3); total (18/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.25e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01 |
TTCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAG > NC_002947/2478739‑2479018 | ttCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTg > 1:731937/1‑150 (MQ=255) cTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCa > 1:260766/1‑150 (MQ=255) cGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGtgttt < 2:143524/150‑1 (MQ=255) tGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgtgtggggg < 2:466636/150‑6 (MQ=38) cttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACca < 2:209928/150‑1 (MQ=255) ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcat < 2:532202/150‑1 (MQ=255) ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCGGCTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCCGGGTGCTTTGTTTTGGAGACcat > 1:480775/1‑150 (MQ=255) tcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcatc < 2:260766/150‑1 (MQ=255) aaaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAc > 2:482529/1‑150 (MQ=255) gatgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt < 2:833205/117‑1 (MQ=38) gatgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt > 1:833205/1‑117 (MQ=38) atgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTc > 1:243068/1‑150 (MQ=255) aCCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt > 1:11380/1‑118 (MQ=255) aCCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt < 2:11380/118‑1 (MQ=255) aTCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATTATCACGAGCGAGTCAGCACGTGGCggg > 2:833689/1‑150 (MQ=255) gCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgatcgat < 2:537836/150‑1 (MQ=255) gCTGTAGCGTGCACGAGAGCTCCGCGTCAGCCCGACCAGGGTAACGTGTTACTCCCCACACCATGCGAGCGCTGAGCCCTCTAGCTGACATTTTCTTGGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGCACCTATGTCGGTAAc > 1:619558/1‑150 (MQ=255) aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga < 2:232344/89‑1 (MQ=255) aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga > 1:232344/1‑89 (MQ=255) cTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg > 2:846496/1‑150 (MQ=255) tCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATgg > 2:658065/1‑150 (MQ=255) tAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAgg < 2:243068/150‑1 (MQ=255) ccGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACg > 1:425998/1‑150 (MQ=255) aTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg < 2:118312/131‑1 (MQ=255) aTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg > 1:118312/1‑131 (MQ=255) cGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCGGGCATGGCCAAGGTACGTTGGTCCAGGCgtt < 2:480775/150‑3 (MQ=255) ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGCCCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaa > 1:745517/1‑150 (MQ=255) cccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGCGTCAGCACGTGCCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTGCAAGCGGTGGTTCGCaag > 1:318836/1‑150 (MQ=255) cccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag > 2:111937/1‑150 (MQ=255) aaaacacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag > 2:784630/4‑150 (MQ=38) ccTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGc < 2:297931/140‑1 (MQ=255) | TTCGAAAGGCCTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAG > NC_002947/2478739‑2479018 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |