Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,790,934 G→A 8.6% A260A (GCG→GCA gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,790,9340GA 8.6% 163.0 / 3.9 58A260A (GCG→GCAgabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (2/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCG  >  minE/1790866‑1791002
                                                                    |                                                                    
ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                    >  1:1020817/1‑71 (MQ=255)
ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                    >  1:535697/1‑71 (MQ=255)
ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                    >  1:1129395/1‑71 (MQ=255)
 gccgccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACgcgcag                                                                   <  1:3959860/71‑3 (MQ=255)
 gccgccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                    <  1:3916302/70‑1 (MQ=255)
 gccgccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                    <  1:3306941/70‑1 (MQ=255)
    gccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                                 >  1:4179409/1‑70 (MQ=255)
    gccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                                 >  1:306722/1‑70 (MQ=255)
     ccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTa                                                               >  1:4385803/1‑71 (MQ=255)
         ttAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                                  <  1:1010572/64‑1 (MQ=255)
           aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                                  >  1:4697145/1‑62 (MQ=255)
           aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                                  >  1:2294923/1‑62 (MQ=255)
           aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCt                                                         >  1:3460380/1‑71 (MQ=255)
            aTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                   <  1:4261234/60‑1 (MQ=255)
            aTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                   <  1:2753019/60‑1 (MQ=255)
            aTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                   <  1:3118648/60‑1 (MQ=255)
            aTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGACCTGGGCGGTAACGCGCCg                                                                   <  1:5350857/60‑1 (MQ=255)
              ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                   >  1:3167787/1‑58 (MQ=255)
              ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                                 >  1:924634/1‑60 (MQ=255)
              ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                                 >  1:1057846/1‑60 (MQ=255)
               gAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGa                                                     >  1:3594024/1‑71 (MQ=255)
               gAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGa                                                     >  1:3262891/1‑71 (MQ=255)
                 aCAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGa                                                     <  1:3821496/69‑1 (MQ=255)
                 aCAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACg                                                   <  1:1512211/71‑1 (MQ=255)
                     tGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGc                                               <  1:3320921/71‑1 (MQ=255)
                     tGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGc                                               <  1:3327581/71‑1 (MQ=255)
                      cggcgAAAAACATCACCAAAGTGTGTCTGGAATTGGGGGGTAAAGCACCAGCTATCGTAATTGACGATGcc                                              >  1:358294/3‑71 (MQ=255)
                       cgcgAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGcc                                              >  1:3318511/1‑70 (MQ=255)
                       cgcgAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCg                                             <  1:229707/71‑1 (MQ=255)
                       cgcgAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCg                                             <  1:1844167/71‑1 (MQ=255)
                        gcgAAAGACATCAAGAATGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                   <  1:1402171/48‑1 (MQ=255)
                           aaaaACATCACCAAAGTGTGTCTGGAATTGGGGGGTAAAGCACCAGCTATCGTAATGGACGATGCCGatc                                          >  1:4523038/1‑68 (MQ=255)
                           aaaaACATCACCAAAGTGTGTCTGGAATTGGGGGGTAAAGCACCAGCTATCGTAATGGACGATGCCGa                                            <  1:5694601/68‑1 (MQ=255)
                             aGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGa                                                  >  1:3733683/1‑60 (MQ=255)
                             aGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGa                                                  >  1:3740376/1‑60 (MQ=255)
                                cATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGAc                                     <  1:65041/70‑1 (MQ=255)
                                  tCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATg                                                >  1:2660258/1‑57 (MQ=255)
                                   cAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTc                                                          >  1:3943463/1‑46 (MQ=255)
                                   cAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGa                                                  >  1:1855420/1‑54 (MQ=255)
                                   cAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGa                                                  >  1:5101960/1‑54 (MQ=255)
                                    aaGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGAc                                           <  1:1764110/60‑1 (MQ=255)
                                      gAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCg                              >  1:1911374/1‑71 (MQ=255)
                                       aaaGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGt                             >  1:86721/1‑71 (MQ=255)
                                       aaaGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGt                             <  1:2485705/71‑1 (MQ=255)
                                         aGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTgg                           <  1:2463317/71‑1 (MQ=255)
                                               gCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGgcg                     <  1:3464272/71‑1 (MQ=255)
                                                atGGAACTGGGCGGTAAATCACCGCTGATCGTTTTCGATGATGCGGATCTCg                                       <  1:4872575/51‑1 (MQ=255)
                                                 tGGAACTGGGCGGTAAATCACCGCTGATCGTTTTCGATGATGCGGAtctc                                        <  1:1412338/50‑1 (MQ=255)
                                                   gAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCt                 <  1:3867037/71‑1 (MQ=255)
                                                   gAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCt                 >  1:2574603/1‑71 (MQ=255)
                                                           cGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCg          >  1:1222230/1‑70 (MQ=255)
                                                           cGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGa         <  1:60543/71‑1 (MQ=255)
                                                           cGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGa         <  1:3943606/71‑1 (MQ=255)
                                                            ggTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGaa        <  1:5534347/71‑1 (MQ=255)
                                                                aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTg                            <  1:1990434/47‑1 (MQ=255)
                                                                aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAg                        >  1:352267/1‑51 (MQ=255)
                                                                aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAg                        >  1:4139746/1‑51 (MQ=255)
                                                                aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAg                        >  1:577037/1‑51 (MQ=255)
                                                                aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt     <  1:2484708/70‑1 (MQ=255)
                                                                   cgcCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCg  >  1:4024487/1‑70 (MQ=255)
                                                                    gcCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt     <  1:1960514/66‑1 (MQ=255)
                                                                    |                                                                    
GGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATTCCG  >  minE/1790866‑1791002

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: