Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 860,911 G→A 37.7% R330R (CGC→CGT trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE860,9110GA37.7% 43.9 / 43.9 74R330R (CGC→CGTtrpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base G (15/31);  new base A (12/16);  total (27/47)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATA  >  minE/860849‑860961
                                                              |                                                  
tCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGCCGACCGCGTGGGCGTGtt                                              >  1:1950235/1‑69 (MQ=255)
tCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGtt                                              >  1:2212732/1‑69 (MQ=255)
tCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGtt                                              >  1:2477674/1‑69 (MQ=255)
       cAATACGGCTGTCGAGGTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:984551/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:1041054/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:896199/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:526465/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:2546400/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:2775642/69‑1 (MQ=255)
       cAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCaa                                       <  1:2948501/69‑1 (MQ=255)
           aCGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCg                                    >  1:2189140/1‑68 (MQ=255)
           aCGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCg                                    >  1:2196463/1‑68 (MQ=255)
           aCGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCg                                    >  1:2358693/1‑68 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:2882868/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:811797/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:2165148/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:643421/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:549794/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:1414630/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:494994/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:1439586/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:2801693/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:1837426/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:1812409/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:1793868/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:3180724/69‑1 (MQ=255)
              gCTGTCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAAGCGGGTa                                <  1:441548/69‑1 (MQ=255)
              gCTGGCGAGGTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:2471902/69‑1 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:2699539/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:1115292/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:2167220/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:115005/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:1435493/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:153102/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:1750680/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:1859136/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:1863052/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:2178465/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtctc                          >  1:2138887/1‑69 (MQ=255)
                    gaggTCTCTGTCCAGCGAACCATCGGCACGACGACCGCGTGGACGTGTTCCGGCAATCGGGTAAAtct                           >  1:3126606/1‑68 (MQ=255)
                    gagaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCgg                                   >  1:283142/1‑60 (MQ=255)
                    gagaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGAtctc                          >  1:965557/1‑69 (MQ=255)
                    gagaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGAtctc                          >  1:954054/1‑69 (MQ=255)
                    gagaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGAtctc                          >  1:318677/1‑69 (MQ=255)
                      gaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCa                         <  1:2860483/68‑1 (MQ=255)
                      gaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCa                         <  1:816996/68‑1 (MQ=255)
                      gaTCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCa                         <  1:935797/68‑1 (MQ=255)
                               ccAGCGAACCATCGGCACGAAGACCGCGTGGACGTGTTCCGGCAATCGGGTa                                <  1:743003/52‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:3360696/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:129112/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:572135/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:1357536/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:77310/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2052773/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:3393977/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2014867/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:1966513/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2581901/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:1970897/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2789706/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2504660/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2434804/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:1668081/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2991669/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:3248748/69‑1 (MQ=255)
                                        cATCGGCGCGACGACCGCGTGGGCGTGTTCCGCCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:2700703/69‑1 (MQ=255)
                                          tCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcatca    <  1:128841/69‑1 (MQ=255)
                                          tCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcatca    <  1:1193683/69‑1 (MQ=255)
                                            ggCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATa  >  1:135535/1‑69 (MQ=255)
                                            ggCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATa  >  1:1586457/1‑69 (MQ=255)
                                            ggCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATa  >  1:1946589/1‑69 (MQ=255)
                                            ggCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATa  >  1:3410176/1‑69 (MQ=255)
                                            ggCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATa  >  1:981437/1‑69 (MQ=255)
                                               gcgACGACCGCGTGGGCGTGCTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGcat      <  1:3341159/62‑1 (MQ=255)
                                                              |                                                  
TCCAGTTCAATACGGCTGTCGAGATCTCTGTCCAGTGAACCATCGGCGCGACGACCGCGTGGGCGTGTTCCGGCAATCGGGTAGATCTCAATCTGGCGGCTGGTGGCATCATA  >  minE/860849‑860961

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: