Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,176 A→T 8.7% S443S (TCA→TCT)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,1760AT 8.7% 221.0 / 11.8 80S443S (TCA→TCT) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base A (29/44);  new base T (5/2);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCG  >  minE/1202109‑1202246
                                                                   |                                                                      
tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                                     >  1:3018705/3‑71 (MQ=9)
tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                                                     >  1:2223873/3‑71 (MQ=9)
tGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                                     >  1:1956194/1‑71 (MQ=255)
  gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                                                    >  1:1606738/1‑70 (MQ=255)
  gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                                     <  1:58901/69‑1 (MQ=255)
  gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                                                     <  1:3059465/69‑1 (MQ=255)
  gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                                   <  1:2097776/71‑1 (MQ=255)
  gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                                   >  1:1593923/1‑71 (MQ=255)
   aTAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggt                                                                  >  1:3418903/1‑71 (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)
        gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                                                             <  1:3562811/71‑1 (MQ=255)
        gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                                                             <  1:1826045/71‑1 (MQ=255)
           aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                                                   >  1:909735/1‑62 (MQ=255)
           aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                                           <  1:1054546/70‑1 (MQ=255)
           aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                                           >  1:2350951/1‑70 (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)
                    aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                                      >  1:2598076/1‑65 (MQ=17)
                    aTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                                                      >  1:2686122/1‑65 (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)
                            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)
                                         aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                             <  1:1213692/54‑1 (MQ=255)
                                          tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                  >  1:2101311/1‑64 (MQ=255)
                                           gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAc                                 <  1:2774320/64‑1 (MQ=255)
                                            aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                             <  1:3383363/51‑1 (MQ=255)
                                            aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc                          <  1:2159856/70‑1 (MQ=255)
                                            aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt                         <  1:3650946/71‑1 (MQ=255)
                                             aCGCCAAATACCGCGAACCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTc                        >  1:3297575/1‑71 (MQ=255)
                                             aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                             <  1:124844/50‑1 (MQ=255)
                                             aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt                         >  1:2964595/1‑70 (MQ=255)
                                               gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt                         >  1:269641/1‑68 (MQ=255)
                                               gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGc                      >  1:1159513/1‑71 (MQ=255)
                                                 cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc                          <  1:2135566/65‑1 (MQ=255)
                                                  aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                    >  1:2556333/1‑54 (MQ=255)
                                                  aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa                   >  1:2252232/1‑71 (MQ=255)
                                                   aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct                             <  1:772815/60‑1 (MQ=255)
                                                   aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct                             <  1:768172/60‑1 (MQ=255)
                                                   aaTACCGCGAAACCTCACGCGGTGGTATGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAt                  <  1:2852716/71‑1 (MQ=255)
                                                     tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtctt                            <  1:1434611/59‑1 (MQ=255)
                                                       ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTg               <  1:1847171/70‑1 (MQ=255)
                                                       ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTg               <  1:696948/70‑1 (MQ=255)
                                                       ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGc              <  1:2065192/71‑1 (MQ=255)
                                                        cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCa             <  1:469610/71‑1 (MQ=255)
                                                          cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                        >  1:1216221/1‑42 (MQ=255)
                                                          cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAAc           >  1:905063/1‑71 (MQ=255)
                                                           gAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACg          <  1:2691118/71‑1 (MQ=255)
                                                            aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACgg         <  1:29767/71‑1 (MQ=255)
                                                              aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt       <  1:1148823/71‑1 (MQ=255)
                                                              aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt       <  1:3529015/71‑1 (MQ=255)
                                                                cTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATgg     <  1:14523/71‑1 (MQ=255)
                                                                 tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGc              <  1:3208140/61‑1 (MQ=255)
                                                                  cACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt                               <  1:2563721/43‑1 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg    >  1:3865182/1‑69 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg    <  1:2209835/69‑1 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg    >  1:1617158/1‑69 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGc   >  1:1671026/1‑70 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCg  <  1:3951400/71‑1 (MQ=255)
                                                                   aCGCGGTGGTCTGGCAATCAAAGTCAAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATggg    >  1:2082567/1‑69 (MQ=255)
                                                                   |                                                                      
TGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCG  >  minE/1202109‑1202246

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