Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,128 G→C 14.0% T427T (ACG→ACC)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1280GC14.0% 191.1 / 14.6 71T427T (ACG→ACC) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (36/24);  new base C (8/2);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCA  >  minE/1202058‑1202195
                                                                      |                                                                   
ttAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACg                                                                     <  1:958400/71‑1 (MQ=255)
 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                    >  1:3798014/1‑71 (MQ=255)
 tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                    >  1:3690374/1‑71 (MQ=255)
    cgccgcGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                  >  1:3460428/1‑70 (MQ=255)
           cgGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAAc                                                          <  1:3704162/71‑1 (MQ=255)
              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)
                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGt                                                                 <  1:1713281/57‑1 (MQ=255)
                  ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAAc                                                          <  1:2093814/64‑1 (MQ=255)
                    ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                                                               >  1:2480085/1‑57 (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)
                           gcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACg                                                                     >  1:2776557/1‑44 (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:2251408/1‑71 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                 >  1:2919357/1‑71 (MQ=255)
                                    cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                 >  1:1981145/1‑71 (MQ=255)
                                     aCCTCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                  <  1:3388292/37‑1 (MQ=255)
                                     aCCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                  <  1:3760320/37‑1 (MQ=255)
                                     aCCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                <  1:2127135/71‑1 (MQ=255)
                                      ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                               >  1:1472404/1‑71 (MQ=255)
                                        gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGaaa                                                           >  1:1866856/1‑41 (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:2223873/3‑71 (MQ=9)
                                                   tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                  >  1:3018705/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:58901/69‑1 (MQ=255)
                                                     gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                  <  1:3059465/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:769621/1‑36 (MQ=255)
                                                           gTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                             >  1:2342920/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:3445904/1‑44 (MQ=255)
                                                               tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                 >  1:3040285/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)
                                                                      |                                                                   
TTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCA  >  minE/1202058‑1202195

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