Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,149 C→T 18.6% D434D (GAC→GAT)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1490CT18.6% 188.3 / 28.7 75D434D (GAC→GAT) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (11/3);  total (46/29)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  minE/1202081‑1202219
                                                                    |                                                                      
cgccgcACCAGTGCACCGCGCTTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                      <  1:820052/71‑1 (MQ=255)
cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                      >  1:2304702/1‑71 (MQ=255)
  ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    <  1:2055974/71‑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:2919357/1‑71 (MQ=255)
             cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                         >  1:2251408/1‑71 (MQ=255)
             cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                         >  1:1981145/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)
                      tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGAc                                                                        >  1:3901268/1‑47 (MQ=255)
                      tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAcc                                                >  1:3313604/1‑71 (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:1054546/70‑1 (MQ=255)
                                       aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                >  1:2350951/1‑70 (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:2686122/1‑65 (MQ=17)
                                                aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                           >  1:2598076/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)
                                                             gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                   >  1:3270597/1‑61 (MQ=255)
                                                             gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                   >  1:198561/1‑61 (MQ=255)
                                                             gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt           <  1:412999/69‑1 (MQ=255)
                                                              tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAGTGTGACt           <  1:830853/68‑1 (MQ=255)
                                                              tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa         >  1:457045/1‑70 (MQ=255)
                                                                aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc      >  1:1262430/1‑71 (MQ=255)
                                                                aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc      <  1:2996347/71‑1 (MQ=255)
                                                                  gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt    <  1:1824390/71‑1 (MQ=255)
                                                                  gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt    >  1:3215150/1‑71 (MQ=255)
                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:1774333/58‑1 (MQ=255)
                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct  >  1:1631395/1‑71 (MQ=255)
                                                                    |                                                                      
CGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  minE/1202081‑1202219

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