Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,140 C→A 20.9% T431T (ACC→ACA sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1400CA20.9% 209.7 / 25.8 76T431T (ACC→ACAsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base A (13/3);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 8.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGT  >  minE/1202072‑1202206
                                                                    |                                                                  
aTGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAgg                                                                  >  1:2502641/1‑71 (MQ=12)
 tGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCgg                                                                  >  1:1529851/1‑70 (MQ=255)
   gcgctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTa                                                                >  1:3611106/4‑70 (MQ=12)
   gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTaa                                                               <  1:3319701/71‑1 (MQ=255)
   gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTaa                                                               <  1:2327117/71‑1 (MQ=255)
      ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAgg                                                             >  1:1259692/1‑70 (MQ=255)
         cgccgcACCAGTGCACCGCGCTTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                         <  1:820052/71‑1 (MQ=255)
         cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                         >  1:2304702/1‑71 (MQ=255)
           ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                       <  1:2055974/71‑1 (MQ=255)
             gcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTaa                                                               <  1:271335/61‑1 (MQ=255)
                 cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                 <  1:1226728/71‑1 (MQ=255)
                   gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                               >  1:2548114/1‑71 (MQ=255)
                     gagccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTa                                             >  1:1773298/4‑71 (MQ=9)
                      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                            >  1:1981145/1‑71 (MQ=255)
                      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                            >  1:2919357/1‑71 (MQ=255)
                      cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                            >  1:2251408/1‑71 (MQ=255)
                       aCCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                           <  1:2127135/71‑1 (MQ=255)
                        ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                                          >  1:1472404/1‑71 (MQ=255)
                           cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGa                                       <  1:3355599/71‑1 (MQ=9)
                           cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                       <  1:3422741/71‑1 (MQ=255)
                              gTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                            <  1:3237355/47‑1 (MQ=255)
                               tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGAc                                                           >  1:3901268/1‑47 (MQ=255)
                               tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAcc                                   >  1:3313604/1‑71 (MQ=255)
                                 tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                              >  1:1142902/1‑42 (MQ=255)
                                  cGCTGGATAAGGTCATCGAAACGATCTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                     >  1:499541/1‑66 (MQ=255)
                                     tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                             >  1:3018705/3‑71 (MQ=9)
                                     tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                             >  1:2223873/3‑71 (MQ=9)
                                     tGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                             >  1:1956194/1‑71 (MQ=255)
                                      ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                  >  1:1951702/1‑49 (MQ=255)
                                      cgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGaaa                                     >  1:852263/2‑62 (MQ=14)
                                       gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                            >  1:1606738/1‑70 (MQ=255)
                                       gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                             <  1:3059465/69‑1 (MQ=255)
                                       gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                             <  1:58901/69‑1 (MQ=255)
                                       gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                           >  1:1593923/1‑71 (MQ=255)
                                       gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                           <  1:2097776/71‑1 (MQ=255)
                                        aTAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggt                          >  1:3418903/1‑71 (MQ=255)
                                         tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                 >  1:902578/1‑47 (MQ=255)
                                         tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATCAACGCCAAATACCGCGAAACCTCACGCggtg                         >  1:1263373/1‑71 (MQ=255)
                                         tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                         >  1:263893/1‑69 (MQ=9)
                                            ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                           <  1:1367076/66‑1 (MQ=255)
                                            ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                      >  1:935991/1‑71 (MQ=255)
                                            ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                      <  1:626062/71‑1 (MQ=255)
                                             gttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCt                     <  1:2906174/68‑1 (MQ=12)
                                             gTCATCGAAACTATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                     <  1:2756947/71‑1 (MQ=255)
                                             gTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                        >  1:2342920/1‑36 (MQ=255)
                                             gTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                        >  1:769621/1‑36 (MQ=255)
                                             gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                     <  1:3562811/71‑1 (MQ=255)
                                             gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                     <  1:1826045/71‑1 (MQ=255)
                                              tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                             >  1:2438567/1‑46 (MQ=255)
                                              tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                           >  1:3059730/1‑48 (MQ=255)
                                               taTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                                         >  1:1111565/2‑49 (MQ=17)
                                                aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                           >  1:909735/1‑62 (MQ=255)
                                                aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                   >  1:2350951/1‑70 (MQ=255)
                                                aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                   <  1:1054546/70‑1 (MQ=255)
                                                 tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                            >  1:3040285/1‑44 (MQ=255)
                                                 tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                            >  1:3445904/1‑44 (MQ=255)
                                                   gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                >  1:3831990/1‑70 (MQ=255)
                                                    aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt               >  1:1703542/1‑70 (MQ=255)
                                                     aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                      <  1:2609792/62‑1 (MQ=255)
                                                     aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa             <  1:2277316/71‑1 (MQ=255)
                                                       cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa             >  1:121230/1‑69 (MQ=255)
                                                        caTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt               <  1:2067284/65‑1 (MQ=17)
                                                         aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                                         >  1:1203697/1‑39 (MQ=17)
                                                         aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg              >  1:2598076/1‑65 (MQ=17)
                                                         aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg              >  1:2686122/1‑65 (MQ=17)
                                                          tGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgcg                                        >  1:1984711/1‑39 (MQ=17)
                                                           gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc        >  1:718376/1‑70 (MQ=255)
                                                            tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATc              <  1:3846907/63‑1 (MQ=255)
                                                            tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa             <  1:2684182/64‑1 (MQ=255)
                                                            tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa      <  1:491799/71‑1 (MQ=255)
                                                             aCGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg                         <  1:1528441/51‑1 (MQ=255)
                                                              cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt    >  1:2844133/1‑71 (MQ=255)
                                                               gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtg   >  1:2327289/1‑71 (MQ=255)
                                                                aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt  <  1:45108/71‑1 (MQ=255)
                                                                aaaCAGGTAAAGATATGAACGCCAAGTACCGCGaa                                      >  1:2588943/1‑35 (MQ=9)
                                                                 aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt  <  1:3421238/70‑1 (MQ=255)
                                                                    |                                                                  
ATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGT  >  minE/1202072‑1202206

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: