Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,117,194 A→G 54.8% intergenic (+22/‑110) aroH → / → ydiE 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,117,1940AG54.8% 7.7 / 40.7 53intergenic (+22/‑110)aroH/ydiE3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base G (12/17);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGC  >  minE/1117135‑1117253
                                                           |                                                           
tCGAAAAACTCGCCTCTGCGGTAGATACTCGCTTCTGAATGCGTGCCCATTCCTGACGGGATGGGCAtt                                                    >  1:274560/1‑69 (MQ=255)
tCGAAAAACTCGCCTCTGCGGTAGATACTCGCTTCTGAATGCGTGCCCATTCCTGACGGGATGGGCAtt                                                    >  1:1793705/1‑69 (MQ=255)
tCGAAAAACTCGCCTCTGCGGTAGATACTCGCTTCTGAATGCGTGCCCATTCCTGACGGGATGGGCAtt                                                    >  1:686790/1‑69 (MQ=255)
             ctctGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                        <  1:311834/68‑1 (MQ=255)
                    gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTtgttgt                                   <  1:997972/66‑1 (MQ=255)
                    gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:2933069/69‑1 (MQ=255)
                    gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:2218834/69‑1 (MQ=255)
                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:1547699/67‑1 (MQ=255)
                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:1040529/67‑1 (MQ=255)
                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:3287445/67‑1 (MQ=255)
                      aGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:95036/67‑1 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCAtt                                                    >  1:2760900/1‑44 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:2648179/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:422511/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:2439024/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:2763845/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:2060984/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTca                            >  1:2150223/1‑68 (MQ=255)
                         tACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCAGTTCTGCGCAACTtgttgt                                   >  1:3107819/1‑61 (MQ=255)
                                ttCTGAATGCGTGCCAATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:391345/57‑1 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:3035213/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:7210/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:2915423/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:3180966/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:722004/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:3391082/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:264578/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:472965/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:2380434/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:2244428/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:2229450/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:1069338/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:2096253/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:1722916/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:121977/68‑2 (MQ=255)
                                    gAATGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:1339015/68‑2 (MQ=255)
                                     aaTGCGTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACta                 <  1:3057227/67‑2 (MQ=255)
                                      aTGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACtt                                        <  1:2481625/43‑1 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCaa                                           >  1:1347529/1‑36 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:1068442/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:889225/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:847072/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:679305/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:1533606/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:154806/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:1619587/1‑69 (MQ=255)
                                          gTGCCCATTCCTGACGGGATGGGCATTTCTGAGCAACTTGTTGTCTTCTCAACAAATTACTACTTGctc          >  1:2575769/1‑69 (MQ=255)
                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTg        >  1:674164/1‑62 (MQ=255)
                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGc  >  1:828017/1‑68 (MQ=255)
                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGc  >  1:1833577/1‑68 (MQ=255)
                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGc  >  1:1895893/1‑68 (MQ=255)
                                                   ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGc  >  1:2210956/1‑68 (MQ=255)
                                                    cTGACGGAATGGGCATTTCTGCGCAACTTGTTGtctt                                <  1:1543716/37‑1 (MQ=255)
                                                           |                                                           
TCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGC  >  minE/1117135‑1117253

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: