Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,542,008 | G→C | 100% | S19S (TCC→TCG) | sixA ← | phosphohistidine phosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,542,008 | 0 | G | C | 75.0% | 57.4 / 17.8 | 40 | S19S (TCC→TCG) | sixA | phosphohistidine phosphatase |
Reads supporting (aligned to +/- strand): ref base G (10/0); new base C (0/30); total (10/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTTCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATA > minE/1541948‑1542065 | tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGTACGAACCGAATc < 1:2014585/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACTAACCGAATc < 1:166276/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:3348598/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:2120826/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:24710/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:2491299/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:2581633/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:3124389/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:3124779/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:556840/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:682859/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:889118/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:1182104/66‑1 (MQ=255) tttCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:1077675/66‑1 (MQ=255) aGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCg > 1:1090851/1‑69 (MQ=255) aGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCg > 1:2835688/1‑69 (MQ=255) aGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCg > 1:928325/1‑69 (MQ=255) gCGAGATTCGTCACAACCGTTAGTAGTCAGAGGACGAACCGAATc < 1:2101223/45‑1 (MQ=38) aTTCGTCACAACCATTAGTGGTCATAGGACGAACGGAATCACTGGCGGCATCGAGGCCTGCGTCGCCGt > 1:1757564/1‑69 (MQ=255) aTTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGt > 1:2295796/1‑69 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:1472700/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2218561/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:1753059/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:810840/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:721416/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:1845703/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:530441/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:37995/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:195696/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2005606/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2818331/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2612827/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2049844/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2568257/37‑1 (MQ=255) acaACCGTTAGTAGTCAGAGGACGAACCGAATCACTg < 1:2087985/37‑1 (MQ=255) cATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATa > 1:1978453/1‑68 (MQ=255) cATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATa > 1:3271550/1‑68 (MQ=255) cATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATa > 1:1899813/1‑68 (MQ=255) cATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATa > 1:2472273/1‑68 (MQ=255) agagGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCata > 1:2874775/1‑69 (MQ=255) | TTTCAGCCAGTTCGCCATCAGGCGAGATTCGTCACAACCATTAGTGGTCAGAGGACGAACGGAATCACTGGCGGCATCGAGGGCTGCGTCGCCGTGACGCATGATAAAAACTTGCATA > minE/1541948‑1542065 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |