Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,146 G→A 18.7% K433K (AAG→AAA)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1460GA18.7% 187.6 / 29.8 75K433K (AAG→AAA) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (35/26);  new base A (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.46e-01

CTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAC  >  minE/1202078‑1202215
                                                                    |                                                                     
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)
           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)
                        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:626062/71‑1 (MQ=255)
                                      ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTc                               >  1:935991/1‑71 (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: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: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)
                                                                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)
                                                                    |                                                                     
CTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAC  >  minE/1202078‑1202215

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