Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,236,903 G→A 27.7% 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,9030GA27.7% 108.4 / 45.9 65A150T (GCT→ACT) purDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base A (7/11);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGC  >  minE/2236842‑2236973
                                                             |                                                                      
aCCTGCGCTGGCGTATCTGCGTGAGAAAGTCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa                                                               <  1:3042397/71‑1 (MQ=255)
aCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa                                                               <  1:1098481/71‑1 (MQ=255)
 ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa                                                               >  1:2167761/1‑70 (MQ=255)
 ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa                                                               >  1:787939/1‑70 (MQ=255)
          gCGTATCTGCGTGAGAAAGGCGGGCCAATCGTCATTAAGGCGGACGGTCTGGCTGCCGGGaaa                                                             >  1:669681/1‑63 (MQ=255)
          gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa                                                             >  1:3328702/1‑63 (MQ=255)
          gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                     >  1:3436195/1‑71 (MQ=255)
          gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                     >  1:2042696/1‑71 (MQ=255)
          gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt                                                     >  1:3482935/1‑71 (MQ=255)
                 tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGtt                                                       <  1:331035/62‑1 (MQ=255)
                 tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGtt                                                       <  1:412297/62‑1 (MQ=255)
                 tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               >  1:2671843/1‑70 (MQ=255)
                 tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGc                                               >  1:435384/1‑70 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               <  1:841690/68‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               <  1:1143134/68‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               <  1:2311075/68‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               <  1:2202056/68‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc                                               <  1:2165099/68‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                            <  1:2239992/71‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                            <  1:2155433/71‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                            <  1:3714704/71‑1 (MQ=255)
                   cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt                                            <  1:3935053/71‑1 (MQ=255)
                     tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                          >  1:1958875/1‑71 (MQ=255)
                     tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                          >  1:2657215/1‑71 (MQ=255)
                     tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa                                          >  1:1436898/1‑71 (MQ=255)
                      gagaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGa                                          <  1:2479587/70‑1 (MQ=255)
                           aGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTgg                                    <  1:3385566/71‑1 (MQ=255)
                           aGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTgg                                    <  1:1800926/71‑1 (MQ=255)
                              cgcgcCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaag                                 <  1:2965709/71‑1 (MQ=255)
                                cgcCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaaga                                <  1:1701084/70‑1 (MQ=255)
                                     aTCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa                                                             <  1:2567473/36‑1 (MQ=255)
                                      tCGTCATTAAAGCTGACGGTCTGGCTGCCGGGaaa                                                             <  1:3482463/35‑1 (MQ=255)
                                      tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt                                                     <  1:3051856/43‑1 (MQ=255)
                                       cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaag                        >  1:3292206/1‑71 (MQ=255)
                                       cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaa                         >  1:926289/1‑70 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt                                                     <  1:1046668/41‑1 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCg                                                   >  1:3634237/1‑43 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCg                                                   >  1:1544179/1‑43 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg                                           >  1:3217465/1‑51 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg                                           >  1:1468003/1‑51 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg                                           >  1:235601/1‑51 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg                                           >  1:2405440/1‑51 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc                       >  1:279236/1‑71 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc                       >  1:2046206/1‑71 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc                       >  1:2338066/1‑71 (MQ=255)
                                        gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGAGTTATCg                                                   >  1:2035120/1‑43 (MQ=255)
                                        gTCATAAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg                                           >  1:2108892/1‑51 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:1215150/38‑1 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:868098/38‑1 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:2899263/38‑1 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:4021/38‑1 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:37305/38‑1 (MQ=255)
                                         tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt                                                       <  1:3814224/38‑1 (MQ=255)
                                          cATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg                      >  1:4118916/1‑70 (MQ=255)
                                          cATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg                     >  1:3002319/1‑71 (MQ=255)
                                              aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt                 >  1:2496986/1‑71 (MQ=255)
                                                    gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa           <  1:2042139/71‑1 (MQ=255)
                                                    gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa           <  1:2487637/71‑1 (MQ=255)
                                                    gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa           <  1:3126265/71‑1 (MQ=255)
                                                    gACGGTATGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa           >  1:531738/1‑71 (MQ=255)
                                                      cGGTCTGGCTGCCGGGAAAGGCGTTATCGTGTCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGata         >  1:3164264/1‑71 (MQ=255)
                                                           tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg    >  1:1865823/1‑71 (MQ=255)
                                                           tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg    >  1:448438/1‑71 (MQ=255)
                                                           tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg    >  1:3237202/1‑71 (MQ=255)
                                                           tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg    >  1:2600266/1‑71 (MQ=255)
                                                             gCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt     <  1:808172/68‑1 (MQ=255)
                                                             gCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGc  >  1:1515742/1‑71 (MQ=255)
                                                             |                                                                      
ACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGC  >  minE/2236842‑2236973

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: