Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,093,727 A→G 7.7% H375H (CAT→CAC)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,093,7270AG 7.7% 255.5 / 11.2 91H375H (CAT→CAC) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base A (47/37);  new base G (4/3);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTG  >  minE/2093658‑2093797
                                                                     |                                                                      
aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                                        <  1:3972130/70‑1 (MQ=255)
aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                                       >  1:2390965/1‑71 (MQ=255)
aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                                       >  1:3120283/1‑71 (MQ=255)
aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                                       >  1:3337251/1‑71 (MQ=255)
 aTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                                      >  1:4009992/1‑71 (MQ=255)
    gcgttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTc                                                                   >  1:69646/4‑71 (MQ=255)
    gCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTc                                                                   <  1:2627347/71‑1 (MQ=255)
     cATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTc                                                                   >  1:1579481/1‑70 (MQ=255)
       ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCt                                                                    <  1:2258390/67‑1 (MQ=255)
         aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                              >  1:699692/1‑71 (MQ=255)
         aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAtgg                                                              >  1:2228876/1‑69 (MQ=255)
          cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                              >  1:1117471/1‑70 (MQ=255)
            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                                      >  1:2011541/1‑60 (MQ=255)
            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                                      >  1:858876/1‑60 (MQ=255)
            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGa                                                           <  1:2503273/71‑1 (MQ=255)
             tCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAt                                                          >  1:2310371/1‑71 (MQ=255)
              cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                                         >  1:2129286/1‑71 (MQ=255)
              cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                                         >  1:4093872/1‑71 (MQ=255)
                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                      <  1:1536656/71‑1 (MQ=255)
                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                      >  1:2888998/1‑71 (MQ=255)
                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                                      >  1:2745009/1‑71 (MQ=255)
                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                                       >  1:1670312/1‑70 (MQ=255)
                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                    >  1:369025/1‑71 (MQ=255)
                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                    >  1:3242848/1‑71 (MQ=255)
                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                                    >  1:2305667/1‑71 (MQ=255)
                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccg                                                     <  1:397361/70‑1 (MQ=255)
                    aCGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgcc                                                   >  1:3321422/1‑71 (MQ=255)
                     cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg                                                  >  1:3809978/1‑71 (MQ=255)
                     cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg                                                  <  1:2485952/71‑1 (MQ=255)
                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCAt                                                                >  1:2908358/1‑55 (MQ=255)
                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                                                <  1:1040543/71‑1 (MQ=255)
                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                                                <  1:859819/71‑1 (MQ=255)
                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                                                >  1:2310296/1‑71 (MQ=255)
                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                                                <  1:1789189/71‑1 (MQ=255)
                        gATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                                               >  1:2266082/1‑71 (MQ=255)
                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:2655062/71‑1 (MQ=255)
                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:392822/71‑1 (MQ=255)
                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:3079567/71‑1 (MQ=255)
                          tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:1633957/70‑1 (MQ=255)
                           ctGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgc                                                 <  1:101106/65‑1 (MQ=12)
                            tGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgc                                                 <  1:766565/65‑1 (MQ=12)
                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                                       >  1:4073972/1‑43 (MQ=255)
                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                                       >  1:2712877/1‑59 (MQ=255)
                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              >  1:2051032/1‑68 (MQ=255)
                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                                           >  1:3514626/1‑71 (MQ=255)
                             gTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCa                                           >  1:1487172/1‑70 (MQ=255)
                             gTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                                           >  1:3698576/1‑70 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCcgccgc                                                 <  1:3578752/63‑1 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            <  1:769223/68‑1 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            <  1:3268612/68‑1 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            <  1:1918185/68‑1 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            <  1:3447487/68‑1 (MQ=255)
                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATa                                         <  1:2384791/71‑1 (MQ=255)
                               aCCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              >  1:3439989/1‑65 (MQ=255)
                                ccTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATAcc                                       >  1:570289/1‑71 (MQ=255)
                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                                        >  1:53049/1‑35 (MQ=255)
                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            >  1:503312/1‑63 (MQ=255)
                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            >  1:1330896/1‑63 (MQ=255)
                                    tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcg                                   >  1:1264869/1‑71 (MQ=255)
                                       ggCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                                       <  1:225504/48‑1 (MQ=255)
                                       ggCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                <  1:2758217/71‑1 (MQ=255)
                                        gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                                                <  1:3575496/54‑1 (MQ=255)
                                         cAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGAGATTTc                                                       >  1:4093129/1‑46 (MQ=255)
                                          aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc                             >  1:2078979/1‑71 (MQ=255)
                                           aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                                          <  1:2454734/57‑1 (MQ=255)
                                            ccGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATATCTGCGCCGGACTg                           >  1:2669204/1‑71 (MQ=255)
                                            ccGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCAGCCGCATTGCATACCTGCGCCGGACTg                           >  1:2219078/1‑71 (MQ=255)
                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:4117087/51‑1 (MQ=255)
                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                                 <  1:2440865/64‑1 (MQ=255)
                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgc                          >  1:335540/1‑71 (MQ=255)
                                              ggATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:2891229/50‑1 (MQ=255)
                                                aTCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                                            >  1:1735763/1‑50 (MQ=255)
                                                  cGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                                               >  1:829668/1‑45 (MQ=255)
                                                  cGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                                               >  1:2929045/1‑45 (MQ=255)
                                                    cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                   >  1:3170170/1‑71 (MQ=255)
                                                    cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag                   <  1:667957/71‑1 (MQ=255)
                                                       gTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                                              <  1:3381146/41‑1 (MQ=255)
                                                       gTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                             <  1:1682587/42‑1 (MQ=255)
                                                       gTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtag                <  1:3128044/71‑1 (MQ=255)
                                                       gTCAGCCCAAGGTTATGCACCATCGCGATTTCcgccgc                                                 >  1:3508017/1‑38 (MQ=255)
                                                       gTCAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgcga                                               >  1:37820/1‑38 (MQ=255)
                                                        tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt               >  1:693515/1‑71 (MQ=255)
                                                         cAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt               >  1:3078577/1‑70 (MQ=255)
                                                         cAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt              >  1:1853429/1‑71 (MQ=255)
                                                         cAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgc                                                 <  1:1617252/36‑1 (MQ=9)
                                                          aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt              <  1:2031998/70‑1 (MQ=255)
                                                              cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATAcc                                       >  1:3815246/1‑41 (MQ=255)
                                                              cAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGtt         <  1:3653623/71‑1 (MQ=255)
                                                               aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa        <  1:2616165/71‑1 (MQ=255)
                                                               aaGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTa        <  1:791864/71‑1 (MQ=255)
                                                                aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa       <  1:791155/71‑1 (MQ=255)
                                                                aGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTaa       <  1:4024254/71‑1 (MQ=255)
                                                                 ggTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCg                                <  1:575663/45‑1 (MQ=255)
                                                                     aTGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTGAGTTAACCCTg  >  1:1042683/1‑71 (MQ=255)
                                                                     |                                                                      
AATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTTCAGTTAACCCTG  >  minE/2093658‑2093797

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: ‑