Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,003 C→A 16.6% T283T (ACG→ACT)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,0030CA16.6% 161.2 / 13.5 54T283T (ACG→ACT) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base A (3/6);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATG  >  minE/2093935‑2094072
                                                                    |                                                                     
ggACGACGGAACTTATCGTAATAAGCCAGCACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                                                     <  1:4042543/71‑1 (MQ=255)
   cgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAc                                                                  >  1:1670121/1‑71 (MQ=255)
      cgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAc                                                               >  1:4040782/1‑71 (MQ=255)
       taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                                               <  1:2549982/68‑1 (MQ=255)
           cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcc                                                          >  1:3169939/1‑71 (MQ=255)
              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                       >  1:1998707/1‑71 (MQ=255)
              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                                                       >  1:3507484/1‑71 (MQ=255)
                   aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCa                                                     >  1:30859/1‑68 (MQ=255)
                   aataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                  >  1:2630179/1‑71 (MQ=255)
                    ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCt                                                  >  1:868048/1‑70 (MQ=255)
                    ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTg                                                 >  1:1189465/1‑71 (MQ=255)
                           aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc                                          <  1:262743/71‑1 (MQ=255)
                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTTCGTTttctt                                        <  1:682388/71‑1 (MQ=255)
                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                                               <  1:4121814/64‑1 (MQ=255)
                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        <  1:291785/71‑1 (MQ=255)
                             tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        <  1:4096712/71‑1 (MQ=255)
                              aCTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt                                        >  1:224093/1‑70 (MQ=255)
                               cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                      <  1:327164/71‑2 (MQ=255)
                               cTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTctcttcg                                      <  1:2716984/71‑2 (MQ=255)
                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCa                                      >  1:2928687/1‑71 (MQ=255)
                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACAACGCGCCCGCCAGCTGCGTTTTCTTCa                                      >  1:1124931/1‑71 (MQ=255)
                                 gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                     >  1:361529/1‑70 (MQ=255)
                                 gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc                                     >  1:3214810/1‑70 (MQ=255)
                                  ccGGAATAATGCCTCACGCACCGTTAGTCGGTGCCGTTACCa                                                                >  1:3215485/1‑42 (MQ=255)
                                  ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACt                                    >  1:1250371/1‑70 (MQ=255)
                                  ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                   <  1:139694/71‑1 (MQ=255)
                                   cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa                                  <  1:2098593/71‑1 (MQ=255)
                                       ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                              >  1:3465786/1‑71 (MQ=255)
                                       ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc                              >  1:3190474/1‑71 (MQ=255)
                                           tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAcc                                >  1:3695737/1‑65 (MQ=255)
                                           tGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGTCAgc                          <  1:424534/71‑1 (MQ=255)
                                             ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                         >  1:915853/1‑68 (MQ=255)
                                             ccGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtg                         >  1:2807830/1‑68 (MQ=255)
                                              cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                       <  1:1148017/71‑2 (MQ=255)
                                              cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgca                       <  1:3406803/71‑2 (MQ=255)
                                              cgcacgcaCCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGTTCTCTTCGTTCACTGCCAgtgc                        >  1:3438212/1‑70 (MQ=255)
                                               gcacgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGcgcg                       >  1:2686903/1‑70 (MQ=255)
                                                  cgcacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACa                   >  1:3285942/1‑71 (MQ=255)
                                                   gcacCGTTAGTCGGCGCAGGCACCACGCGACCGCCAGCAGCGTTctct                                         <  1:3950309/48‑1 (MQ=16)
                                                    cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg                 <  1:520067/71‑1 (MQ=255)
                                                    cacCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATg                 <  1:2201729/71‑1 (MQ=255)
                                                     acCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGt                <  1:3827713/71‑1 (MQ=255)
                                                      ccGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct                                         <  1:4082012/45‑1 (MQ=255)
                                                         ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            >  1:1401770/1‑71 (MQ=255)
                                                         ttAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            <  1:2463316/71‑1 (MQ=255)
                                                          tAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            >  1:3042194/1‑70 (MQ=255)
                                                           aGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTa                                   <  1:2193282/46‑1 (MQ=255)
                                                            gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATcc          <  1:2943441/70‑1 (MQ=255)
                                                            gTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCa         <  1:3478794/71‑1 (MQ=255)
                                                             tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGcc                             >  1:4110362/1‑50 (MQ=255)
                                                             tCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAg        <  1:2559309/71‑1 (MQ=255)
                                                                gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGAt            <  1:2510980/64‑1 (MQ=255)
                                                                gTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCg     <  1:3659028/71‑1 (MQ=255)
                                                                 tGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGa    >  1:1423232/1‑71 (MQ=255)
                                                                   ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGAt   <  1:1911099/70‑1 (MQ=255)
                                                                   ccGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATg  <  1:565092/71‑1 (MQ=255)
                                                                    cGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGt                >  1:3777975/1‑56 (MQ=255)
                                                                    |                                                                     
GGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGCCAGCGCGTACATGTTGATCCAGTCGATG  >  minE/2093935‑2094072

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