Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,236,903 G→A 45.9% A150T (GCT→ACT)  purD → phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,236,9030GA45.9% 19.2 / 81.4 61A150T (GCT→ACT) purDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base A (18/10);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT  >  minE/2236837‑2236970
                                                                  |                                                                   
gTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGc                                                                 >  1:2264255/1‑71 (MQ=255)
  aGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCg                                                               >  1:3185880/1‑71 (MQ=255)
     aCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa                                                            <  1:777273/71‑1 (MQ=255)
      ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGaa                                                           <  1:3185941/71‑1 (MQ=255)
      ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa                                                            >  1:153489/1‑70 (MQ=255)
          cgcTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAg                                                         <  1:4009243/69‑1 (MQ=255)
               gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                    <  1:1725880/69‑1 (MQ=255)
                cGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATc                                                 <  1:1248767/71‑1 (MQ=255)
                cGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATc                                                 <  1:2476028/71‑1 (MQ=255)
                   aTCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGAGTGCCGGGa                                                            <  1:366115/57‑1 (MQ=255)
                    tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAgg                                                        <  1:2672240/60‑1 (MQ=255)
                    tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg                                             >  1:3117701/1‑71 (MQ=255)
                    tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg                                             >  1:3839979/1‑71 (MQ=255)
                    tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg                                             >  1:1775064/1‑71 (MQ=255)
                    tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa                                                            <  1:670178/56‑1 (MQ=255)
                      tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCg                                                >  1:1497393/1‑66 (MQ=255)
                       gCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa                                                            <  1:747367/53‑1 (MQ=255)
                        cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                            <  1:1016544/68‑1 (MQ=255)
                        cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                            <  1:2485555/68‑1 (MQ=255)
                            agaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGa                                          >  1:3731541/1‑66 (MQ=255)
                              aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCt                                   >  1:790524/1‑71 (MQ=255)
                              aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTTCCGGGAAAGGCg                                                      >  1:1987108/1‑52 (MQ=255)
                              aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCg                                                      >  1:1086079/1‑52 (MQ=255)
                                  gcgcgcCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa                                                          >  1:812703/1‑44 (MQ=255)
                                         aaTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcg                         >  1:2270006/1‑70 (MQ=255)
                                          aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                  >  1:1954187/1‑44 (MQ=255)
                                          aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                  >  1:1426525/1‑44 (MQ=255)
                                          aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                  >  1:3061269/1‑44 (MQ=255)
                                           tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                         <  1:1535155/52‑1 (MQ=255)
                                           tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                         <  1:1530242/52‑1 (MQ=255)
                                           tCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTg                                  <  1:4051871/59‑1 (MQ=255)
                                           tCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTg                                  <  1:3760827/59‑1 (MQ=255)
                                            cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:536630/1‑53 (MQ=255)
                                            cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:2039352/1‑53 (MQ=255)
                                              tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg                                     <  1:516141/53‑1 (MQ=255)
                                              tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg                                     <  1:967498/53‑1 (MQ=255)
                                              tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg                                     <  1:2863338/53‑1 (MQ=255)
                                              tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg                   <  1:2685062/71‑1 (MQ=255)
                                               cATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg                  >  1:2017247/1‑71 (MQ=255)
                                               cATTAAAGCGGAAGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg                  >  1:1598905/1‑71 (MQ=255)
                                                aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGc                 >  1:741351/1‑71 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaagaa                            <  1:858603/58‑1 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg               >  1:1039613/1‑71 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg               >  1:1710888/1‑71 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg               >  1:617218/1‑71 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg               >  1:3245023/1‑71 (MQ=255)
                                                  tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg               >  1:914242/1‑71 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:4250096/1‑46 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:269238/1‑46 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:1806436/1‑46 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                       >  1:3829756/1‑46 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt              >  1:2567528/1‑71 (MQ=255)
                                                   aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt              >  1:3874952/1‑71 (MQ=255)
                                                    aaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGtt             <  1:690621/71‑1 (MQ=255)
                                                         gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa        >  1:524953/1‑71 (MQ=255)
                                                         gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTAGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa        >  1:649557/1‑71 (MQ=255)
                                                              tCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc   >  1:2859478/1‑71 (MQ=255)
                                                              tCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc   >  1:1514422/1‑71 (MQ=255)
                                                              gcTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc   <  1:2267709/70‑1 (MQ=255)
                                                               cTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt  >  1:4005696/1‑71 (MQ=255)
                                                               cTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt  >  1:3107326/1‑71 (MQ=255)
                                                                  |                                                                   
GTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT  >  minE/2236837‑2236970

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: