Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,073,563 | A→G | 12.1% | L433P (CTG→CCG) | gatA ← | glutamyl‑tRNA(Gln) amidotransferase subunit A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,073,563 | 0 | A | G | 12.1% | 73.1 / 3.6 | 33 | L433P (CTG→CCG) | gatA | glutamyl‑tRNA(Gln) amidotransferase subunit A |
Reads supporting (aligned to +/- strand): ref base A (16/13); new base G (2/2); total (18/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.52e-01 |
AGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC > NC_002947/1073414‑1073712 | aGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCa < 1:790700/150‑1 (MQ=255) ggggCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGcc > 1:374988/1‑150 (MQ=255) gcgGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGtt > 2:706813/1‑150 (MQ=255) aTGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGGGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGt < 1:502194/150‑1 (MQ=255) gCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGa > 1:10187/1‑150 (MQ=255) gCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGa < 1:586082/150‑1 (MQ=255) ggTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCGGGCCTGCCAggc < 1:116995/134‑2 (MQ=255) ggTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCGGGCCTGCCAggc > 2:116995/1‑133 (MQ=255) gTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTa < 2:62899/138‑1 (MQ=255) gTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTa > 1:62899/1‑138 (MQ=255) gTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggc > 1:197592/1‑150 (MQ=255) gTGCGCGACGCTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggc < 2:508373/150‑1 (MQ=255) gcgcGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAggcgg < 1:706813/150‑1 (MQ=255) gtggAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCgg > 1:247232/4‑150 (MQ=255) gCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGGGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCgc > 2:616582/1‑150 (MQ=255) ggcggcCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGGCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCgctgc < 2:247232/150‑1 (MQ=255) gcCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTGGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCtt > 2:157738/1‑150 (MQ=255) ctcccTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTgg < 1:606762/147‑1 (MQ=255) cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATg > 1:698545/1‑104 (MQ=255) cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATg < 2:698545/104‑1 (MQ=255) cTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTgg > 1:441235/1‑150 (MQ=255) cTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCGCTTGgcgtc > 1:336663/1‑148 (MQ=255) cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc > 1:325063/1‑143 (MQ=255) cagcagTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTcc < 2:325063/143‑1 (MQ=255) ttGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGTCGGCACCCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGTCGGGTTTgg < 2:10187/150‑1 (MQ=255) gCAGGCCGTCGACGAAGCTGGCAGGCATCGACAGGGCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGTGGCAATCGGGTCGCTACTCTTGAGGCCGAGCTTCCAGGGCGGGGGTAGCGTGGGCGGGCCgg > 2:155721/1‑149 (MQ=255) ggCCGTCTACGAAGCCCGCAGGCATCGACAGGCATGGCAGGGCTGCCTGTTTGCCGGTGATGGTGTAGACGTCTTCCATATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGaa < 1:616582/150‑1 (MQ=255) tCGACGAAGCCCGCAGGAGTAGATAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAgg < 1:683899/150‑1 (MQ=255) gacgaAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTc > 2:233260/1‑150 (MQ=255) ccGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGTGAACCGGGTCGCTGCGCTTGTCGTCAAGCTTCCAGGGCGGGTTGGGCGTGGTCGGGCCGAGAATGAGCTCGACATCGt > 1:60417/1‑150 (MQ=255) aaccTGGCGGGCCTGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGAt < 2:161547/49‑1 (MQ=25) aaccTGGCGGGCCTGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGAt > 1:161547/3‑51 (MQ=25) gCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGt > 2:452223/1‑125 (MQ=255) gCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGt < 1:452223/125‑1 (MQ=255) aGGCCCGCCAGGTTGGCGGTGACGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGTCCGGGTTTGGCGTGGGCGGGCCGAGAATCAGGTCGACATTGTTGAAGGCGGCCATGAAGTCGTTc < 1:155721/150‑1 (MQ=255) | AGCCGTTGGGGGCGCGGGTATGCCAGTCGGTCACCTGTTGGTAGCGGTGCGCGACGTTGAGCAGGCGGCCTTCCTGGAAGTACGGCGCCAGCAGTTGCACGCCGACCGGCAGGCCGTCGACGAAGCCCGCAGGCATCGACAGGCCTGGCAGGCCCGCCAGGTTGGCGGTGATGGTGTAGACGTCTTCGAGATAGGCGGCAACCGGGTCGCTGCTCTTGGCGCCAAGCTTCCAGGCCGGGTTTGGCGTGGTCGGGCCGAGAATCAGGTCGACATCGTTGAAGGCGGCCATGAAGTCGTTC > NC_002947/1073414‑1073712 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |