| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,672,724 | C→T | 7.3% | A180T (GCG→ACG) | hisS ← | histidyl tRNA synthetase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,672,724 | 0 | C | T | 7.3% | 175.1 / 3.4 | 55 | A180T (GCG→ACG) | hisS | histidyl tRNA synthetase |
| Reads supporting (aligned to +/- strand): ref base C (21/30); new base T (2/2); total (23/32) | |||||||||||
| 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.75e-01 | |||||||||||
AGTGTACATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGAAA > minE/1672656‑1672792 | aGTGTACATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCa < 1:3953437/70‑1 (MQ=255) tACATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGTATcgcg < 1:2352537/71‑1 (MQ=255) cATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCgg > 1:1870763/1‑70 (MQ=255) cATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGt < 1:1869364/71‑1 (MQ=255) cggcgTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAAtt > 1:1181128/1‑71 (MQ=255) tttGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAAtt < 1:3093698/66‑1 (MQ=255) tttGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGgcgc > 1:2492563/1‑71 (MQ=255) tCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATcgc < 1:1668976/52‑1 (MQ=255) tCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTa < 1:271322/56‑1 (MQ=255) tCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGgc < 1:3317518/62‑1 (MQ=255) ttCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAgcgc < 1:12320/45‑1 (MQ=255) ttCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGt < 1:2993907/53‑1 (MQ=255) cGTCCAGCTTTTCTTTATGCTGCTCAATGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTcc < 1:1535456/68‑1 (MQ=255) cGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGc > 1:15185/1‑71 (MQ=255) ccAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGaa > 1:371497/1‑71 (MQ=255) aGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGt < 1:587620/46‑1 (MQ=255) aGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAAcc < 1:3740049/71‑1 (MQ=255) aGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAAcc < 1:3406461/71‑1 (MQ=255) gCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGt < 1:2581129/45‑1 (MQ=255) gCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGt < 1:3995487/45‑1 (MQ=255) tttCTTTATGCTGCTCAAGGAATGCCACCAGCGTATCGCGGTaa > 1:910899/1‑44 (MQ=255) tttCTTTATGCTGCTCAAGGAATGCCACCAGCGTATCGCGGTaa > 1:3370625/1‑44 (MQ=255) tCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAAtt > 1:3088516/1‑44 (MQ=255) tCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAAtt > 1:3083943/1‑44 (MQ=255) cTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATagag > 1:3831230/1‑71 (MQ=255) ttATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGtt > 1:3817229/1‑71 (MQ=255) ttATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGt > 1:527170/1‑70 (MQ=255) tATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGtt < 1:2197144/70‑1 (MQ=255) tgctgcTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAg > 1:199717/1‑71 (MQ=255) ctgctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATagag < 1:1696376/64‑1 (MQ=255) ctgctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATagag < 1:2871874/64‑1 (MQ=255) ctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTc < 1:3437813/69‑1 (MQ=255) ctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTc < 1:1599842/69‑1 (MQ=255) ctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCaa < 1:1080726/71‑1 (MQ=255) ctCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCa < 1:1399421/70‑1 (MQ=255) tCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAg < 1:2414101/47‑1 (MQ=255) cAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCaa > 1:1777962/1‑69 (MQ=255) cAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCaa > 1:1968053/1‑69 (MQ=255) ggAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCa < 1:1687287/65‑1 (MQ=255) ggAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGtt < 1:841439/71‑1 (MQ=255) gAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTa < 1:2978059/71‑1 (MQ=255) aaTGCCACCAGCGTATCGCGGTAATTGGCGCGTGCTTCCATCGaa < 1:174973/45‑1 (MQ=38) aaTGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAg > 1:595537/1‑41 (MQ=255) aaTGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAag > 1:2896776/1‑65 (MQ=255) aaTGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTAc < 1:3801406/71‑1 (MQ=255) tGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCaa < 1:1240294/62‑1 (MQ=255) tGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAaga > 1:3016818/1‑64 (MQ=255) tGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAaga > 1:3294436/1‑64 (MQ=255) gCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCaa > 1:2027078/1‑61 (MQ=255) gCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGtt < 1:2481346/66‑1 (MQ=255) ccaccaGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGTACCGATAGAGTTCAGCTCAAGAGTTACGTGc > 1:912164/1‑71 (MQ=255) caccaGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCt < 1:1348801/71‑1 (MQ=255) accaGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCg < 1:3164505/37‑1 (MQ=255) accaGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTc > 1:3540691/1‑71 (MQ=255) ccaGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCg > 1:1471174/1‑71 (MQ=255) aGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGa < 1:3939526/71‑1 (MQ=255) gcgcATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGaa < 1:1007674/71‑1 (MQ=255) gcgcATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGa < 1:278370/70‑1 (MQ=255) cgcATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGaaa < 1:3397957/71‑1 (MQ=255) cATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATa > 1:1965404/1‑38 (MQ=255) | AGTGTACATGCGGCGTTTGCAGTCTTCGTCCAGCTTTTCTTTATGCTGCTCAAGGAATGCCACCAGCGCATCGCGGTAATTGGCGCGTGCTTCCAGCGAACCGATAGAGTTCAGCTCAAGAGTTACGTGCTCGGAAA > minE/1672656‑1672792 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |