Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,790,919 G→A 7.6% E255E (GAG→GAA 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,9190GA 7.6% 187.1 / 4.2 66E255E (GAG→GAAgabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/3);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT  >  minE/1790854‑1790987
                                                                 |                                                                    
tCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGc                                                                     >  1:4934978/1‑67 (MQ=255)
 cGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGc                                                                >  1:1426253/1‑71 (MQ=255)
 cGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGc                                                                >  1:3364956/1‑71 (MQ=255)
 cGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGc                                                                >  1:3560672/1‑71 (MQ=255)
  gACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTggg                                                                 <  1:3363270/69‑1 (MQ=255)
   aCCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCgg                                                              >  1:3107810/1‑71 (MQ=255)
    ccGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCgg                                                              >  1:733336/1‑70 (MQ=255)
       aaaTTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAAc                                                          <  1:3805694/71‑1 (MQ=255)
          ttGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCt                                                                    <  1:1766586/58‑1 (MQ=255)
            ggCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTgg                                                                  >  1:5009314/1‑58 (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:3306941/70‑1 (MQ=255)
             gccgccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGcc                                                     <  1:3916302/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)
                   aGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTaa                                                           >  1:3259785/1‑58 (MQ=255)
                     ttAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                   <  1:1010572/64‑1 (MQ=255)
                      tAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTaa                                                           >  1:4520204/1‑55 (MQ=255)
                      tAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTaa                                                           >  1:129902/1‑55 (MQ=255)
                       aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCg                                                               >  1:3797872/1‑50 (MQ=255)
                       aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                   >  1:2294923/1‑62 (MQ=255)
                       aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGt                                                   >  1:4697145/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)
                         tGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTg                                                                   <  1:4856354/44‑1 (MQ=255)
                          ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCgg                                                              <  1:5561807/48‑1 (MQ=255)
                          ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCg                                                    >  1:3167787/1‑58 (MQ=255)
                          ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                  >  1:1057846/1‑60 (MQ=255)
                          ggAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGtt                                                  >  1:924634/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:3327581/71‑1 (MQ=255)
                                 tGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGc                                <  1:3320921/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:5101960/1‑54 (MQ=255)
                                               cAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGa                                   >  1:1855420/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)
                                                                 |                                                                    
TCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT  >  minE/1790854‑1790987

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: