Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,648,260 | G→C | 15.0% | intergenic (+16/+8) | elaA → / ← PP_4957 | acetyltransferase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,648,260 | 0 | G | C | 15.0% | 97.9 / 9.5 | 39 | intergenic (+16/+8) | elaA/PP_4957 | acetyltransferase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (14/19); new base C (2/4); total (16/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
GTTGGCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATGT > NC_002947/5648111‑5648407 | gttgGCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCggg < 1:1353490/150‑1 (MQ=255) gggTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacagg < 1:586035/150‑3 (MQ=255) tGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAgg < 1:1286383/150‑1 (MQ=255) tGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAgg < 1:149683/150‑1 (MQ=255) tGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacagggtg < 1:170441/150‑6 (MQ=255) gCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAggg < 1:677411/150‑1 (MQ=255) gCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcg > 1:647457/1‑150 (MQ=255) cTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGCATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc < 2:291204/150‑1 (MQ=255) aTTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacagggtg < 1:1301867/123‑6 (MQ=255) aTTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacagggtg > 2:1301867/1‑118 (MQ=255) aGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGt > 2:1333006/1‑150 (MQ=255) gCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacaggg > 2:951003/1‑110 (MQ=255) gCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGCCCGcacaggg < 1:951003/113‑4 (MQ=255) tttATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCTACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGga < 1:446572/150‑1 (MQ=255) tttATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCAACCTGTGGGGGCCGCCCCAGGGGTGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGga < 2:49802/150‑1 (MQ=255) aTGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc > 1:19286/1‑126 (MQ=255) aTGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTcgc < 2:19286/126‑1 (MQ=255) ccGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCa < 1:65755/150‑1 (MQ=255) gtggGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGGGGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGccccc > 1:1170247/1‑150 (MQ=255) gtggGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCNCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGccccc < 1:1039515/150‑1 (MQ=255) cGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTc > 2:893761/1‑150 (MQ=255) ggagTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATg < 2:450654/150‑1 (MQ=255) tttttGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATg < 2:336053/150‑1 (MQ=255) ttGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCACCAGGTCGGGGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGa < 2:749741/150‑1 (MQ=255) gaggaCGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATaa < 1:1333006/150‑1 (MQ=255) gaggaCGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATaa > 1:378959/1‑150 (MQ=255) gCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGTTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCgg > 2:126275/1‑150 (MQ=255) tCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGAt < 1:862909/115‑1 (MQ=255) tCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGAt > 2:862909/1‑115 (MQ=255) cgTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGTTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGt < 1:126275/150‑1 (MQ=255) gTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCt > 2:1335433/1‑104 (MQ=255) gTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCt < 1:1335433/104‑1 (MQ=255) tCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTtg > 2:1272197/1‑150 (MQ=255) aTTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCCCAACAGGGGTAACCCAGCACTTCGCGGATCGGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTtgcgt > 1:388191/1‑150 (MQ=255) aaGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCctgc > 2:202497/1‑150 (MQ=255) aTGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTc < 2:412914/150‑1 (MQ=255) cccACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTc > 1:660888/1‑136 (MQ=255) cccACCCTGTGCGGGCCGCACAGCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTc < 2:660888/136‑1 (MQ=255) ccACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGt > 1:69181/1‑122 (MQ=255) ccACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGt < 2:69181/122‑1 (MQ=255) gggCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATGt > 1:640493/1‑150 (MQ=255) ggCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATg < 2:493525/91‑1 (MQ=255) ggCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATg > 1:493525/1‑91 (MQ=255) | GTTGGCCGGGGGTGCCGGTTTATCTGTCGGCTCCGGCGCATTTACAGGCGTTTTATGCCCGCCATGGGTTTGCGGTGGTGGGCGAGGAGTTTTTGGAGGACGGCATCGCGTCCATTGGAATGCGCAAGGGATGAGCCCACCCTGTGCGGGCCGCACATCAGGGGTAACCCAGCACTTCGCGGATCTGCCGCAAGTGACGGATGATCCACTGCTTGTCGATCGCCCCCCAGTCATGAATGCGATAATGCCCGGCGTGGTTGCGTGCGCCTTCCTGCTGTTCGAACTCGCAGACGATGT > NC_002947/5648111‑5648407 |
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: ‑ |