Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,790,931 C→A 8.4% N259K (AAC→AAA 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,9310CA 8.4% 186.0 / 3.4 59N259K (AAC→AAAgabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base A (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATT  >  minE/1790861‑1790999
                                                                      |                                                                    
aaaTTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAAc                                                                      <  1:3805694/71‑1 (MQ=255)
     ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                 >  1:1020817/1‑71 (MQ=255)
     ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                 >  1:1129395/1‑71 (MQ=255)
     ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                                 >  1:535697/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:3118648/60‑1 (MQ=255)
                 aTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                                <  1:2753019/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:3262891/1‑71 (MQ=255)
                    gAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGa                                                  >  1:3594024/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:2574603/1‑71 (MQ=255)
                                                        gAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCt              <  1:3867037/71‑1 (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:577037/1‑51 (MQ=255)
                                                                     aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAg                     >  1:4139746/1‑51 (MQ=255)
                                                                     aCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAAtt  <  1:2484708/70‑1 (MQ=255)
                                                                      |                                                                    
AAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCTGGCCTCGAAATT  >  minE/1790861‑1790999

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: