Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,246,455 | T→C | 36.6% | I56I (ATT→ATC) | thiC → | thiamin (pyrimidine moiety) biosynthesis protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,246,455 | 0 | T | C | 36.6% | 39.3 / 43.5 | 71 | I56I (ATT→ATC) | thiC | thiamin (pyrimidine moiety) biosynthesis protein |
Reads supporting (aligned to +/- strand): ref base T (14/31); new base C (8/18); total (22/49) | |||||||||||
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.84e-01 |
TATATCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTC > minE/2246387‑2246517 | tataTCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAAtt > 1:1091433/1‑69 (MQ=255) tataTCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAAtt > 1:3345218/1‑69 (MQ=255) cacacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg < 1:451683/69‑1 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:1747805/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:3096274/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:3344853/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:2897262/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:2594491/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:776149/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:2365402/1‑69 (MQ=255) cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg > 1:104328/1‑69 (MQ=255) acacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg < 1:1835611/68‑1 (MQ=255) acacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg < 1:1694545/68‑1 (MQ=255) cacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg < 1:184800/67‑1 (MQ=255) gcgtgcgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGc < 1:2246164/69‑1 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa > 1:2225698/1‑56 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa > 1:2038279/1‑56 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa < 1:2359927/56‑1 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa < 1:3135629/56‑1 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa > 1:3316878/1‑56 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa > 1:1252622/1‑56 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa < 1:1165627/56‑1 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa > 1:3117512/1‑68 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa > 1:3020654/1‑68 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa > 1:1321960/1‑68 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa > 1:2777417/1‑68 (MQ=255) cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa > 1:1056653/1‑68 (MQ=255) gtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa < 1:433743/53‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:315064/69‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:3171746/69‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:2861871/69‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:496848/69‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:2596257/69‑1 (MQ=255) agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:1348332/69‑1 (MQ=255) aTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa < 1:2827732/67‑1 (MQ=255) ccAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTcc > 1:575560/1‑69 (MQ=255) ccAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTc > 1:2874530/1‑68 (MQ=255) tAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCt > 1:1804548/1‑68 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACTAAGAAAACGAAGCGATTCCGGTCTAc < 1:1117983/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:1716686/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:3235438/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:1235182/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:1965124/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2824673/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2115288/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:665987/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2651602/67‑1 (MQ=255) cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2384376/67‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:660826/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:503640/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:445378/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:685838/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:151068/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:2276408/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:1747325/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:1004102/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:270696/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:1753570/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:2773711/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:2800052/45‑1 (MQ=255) ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:1950466/45‑1 (MQ=255) ccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:141613/66‑1 (MQ=255) cGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:139209/44‑1 (MQ=255) cGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:134344/65‑1 (MQ=255) gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:784606/64‑1 (MQ=255) gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:566420/64‑1 (MQ=255) gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:1944916/64‑1 (MQ=255) gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:136276/64‑1 (MQ=255) gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2328650/64‑1 (MQ=255) gCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaaga > 1:3318277/1‑39 (MQ=255) gCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTc > 1:940574/1‑69 (MQ=255) cTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc < 1:2889076/60‑1 (MQ=255) aTTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa < 1:3222601/36‑1 (MQ=255) | TATATCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTC > minE/2246387‑2246517 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |