Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,146 G→A 7.4% 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,1460GA 7.4% 305.9 / 10.0 95K433K (AAG→AAA) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base A (3/4);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAT  >  minE/1202076‑1202213
                                                                      |                                                                   
ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                     <  1:3473038/71‑1 (MQ=255)
ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                     >  1:3828835/1‑71 (MQ=255)
ctctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                     <  1:1463917/71‑1 (MQ=255)
  ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                                                   <  1:712785/71‑1 (MQ=12)
  ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                   >  1:2208185/1‑71 (MQ=255)
             cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                        <  1:2452681/71‑1 (MQ=255)
              aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                        <  1:1034642/70‑1 (MQ=255)
               gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                >  1:2898195/1‑61 (MQ=255)
               gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                >  1:2380527/1‑61 (MQ=255)
               gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                      <  1:2126291/71‑1 (MQ=255)
               gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                       >  1:489676/1‑70 (MQ=255)
                  cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                   >  1:2019384/1‑71 (MQ=255)
                  cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                   <  1:2214042/71‑1 (MQ=255)
                  cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                                                   >  1:1360214/1‑71 (MQ=255)
                  cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAAAc                                                   >  1:2187312/1‑71 (MQ=255)
                    ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                    >  1:2645408/1‑68 (MQ=255)
                    ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                                  <  1:3132007/70‑1 (MQ=255)
                    ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcg                                                 >  1:1763813/1‑71 (MQ=255)
                      gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                       >  1:375368/1‑63 (MQ=255)
                       cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                         >  1:326770/1‑60 (MQ=255)
                       cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                              <  1:2960531/71‑1 (MQ=255)
                            ctcGCTGGATAAGGTCATCGAATCGATGTACGAAACCGGTAAGGa                                                                   <  1:2101549/45‑1 (MQ=255)
                            ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                   <  1:3929155/45‑1 (MQ=255)
                            ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                                                 >  1:2464782/1‑47 (MQ=255)
                              cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                                       >  1:3246103/1‑71 (MQ=255)
                               gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAc                                      <  1:2278689/71‑1 (MQ=255)
                                cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                     <  1:1667805/71‑1 (MQ=255)
                                cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                     <  1:901792/71‑1 (MQ=255)
                                  ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                   <  1:4065938/39‑1 (MQ=255)
                                   gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                                             >  1:2266177/1‑44 (MQ=255)
                                   gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                   <  1:190977/70‑1 (MQ=255)
                                   gATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                     <  1:3972827/68‑1 (MQ=255)
                                     tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                                >  1:461654/1‑69 (MQ=9)
                                     tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCggcg                                <  1:382465/71‑3 (MQ=9)
                                      aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                               <  1:355144/39‑1 (MQ=255)
                                      aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                           >  1:1684957/1‑59 (MQ=255)
                                      aaGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                               >  1:2299572/1‑71 (MQ=255)
                                       aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                             <  1:1213629/56‑1 (MQ=255)
                                       aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                               <  1:969779/70‑1 (MQ=255)
                                         gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                        <  1:4137919/43‑1 (MQ=255)
                                          ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                        <  1:2694751/40‑1 (MQ=17)
                                          ttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                        <  1:310864/40‑1 (MQ=17)
                                          tCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTg                           <  1:2297868/71‑1 (MQ=255)
                                            aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCt                            >  1:3159503/1‑68 (MQ=255)
                                             tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                        <  1:3073080/39‑1 (MQ=255)
                                             tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                                         >  1:3004455/1‑54 (MQ=255)
                                             tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                    <  1:142465/59‑1 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                      >  1:298657/1‑39 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                            >  1:2223297/1‑49 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaaa                                            >  1:2227986/1‑49 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                           >  1:2673881/1‑50 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                      >  1:4121969/1‑71 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                      >  1:3157823/1‑71 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                      >  1:3602236/1‑71 (MQ=255)
                                               gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACAGCGAAAc                                           >  1:3177867/1‑50 (MQ=255)
                                               gAAACCATGTACGAAACAGGTAAAGATATGAACGcc                                                         >  1:1676812/1‑36 (MQ=12)
                                               gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                      >  1:743552/1‑71 (MQ=17)
                                                 aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    <  1:2023539/71‑1 (MQ=255)
                                                  aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAAc                                           >  1:2836474/1‑47 (MQ=255)
                                                  aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                         >  1:562780/1‑65 (MQ=255)
                                                  aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaa                   <  1:1724198/71‑1 (MQ=255)
                                                  aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    >  1:3562572/1‑70 (MQ=255)
                                                  aCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    >  1:1993506/1‑70 (MQ=255)
                                                   cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    >  1:4007118/1‑69 (MQ=255)
                                                   cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    >  1:2792915/1‑69 (MQ=255)
                                                    gATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCt                                         >  1:3024897/1‑47 (MQ=255)
                                                     aTGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTg                           <  1:294420/60‑1 (MQ=255)
                                                      tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc               <  1:1014421/71‑1 (MQ=255)
                                                       gTACGAAACCGGTAAGGACATGAACGCCAAATAcc                                                  <  1:2176499/35‑1 (MQ=255)
                                                       gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                       <  1:3760122/62‑1 (MQ=255)
                                                        tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTcc              <  1:3625702/70‑1 (MQ=255)
                                                         aCGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg            >  1:1288104/1‑71 (MQ=255)
                                                          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt           >  1:3444359/1‑71 (MQ=255)
                                                          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt           >  1:3431877/1‑71 (MQ=255)
                                                          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg            >  1:738033/1‑70 (MQ=255)
                                                          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGATCTGGCAATCa                    <  1:2211954/62‑1 (MQ=255)
                                                          agAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt           <  1:2425418/70‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAgtgt         <  1:432932/71‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt         <  1:2996298/71‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt         <  1:2920582/71‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt         <  1:318213/71‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt         <  1:1728372/71‑1 (MQ=255)
                                                            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt         <  1:2243705/71‑1 (MQ=255)
                                                             aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                          >  1:2512439/1‑53 (MQ=255)
                                                             aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                         >  1:3434326/1‑54 (MQ=255)
                                                             aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc               <  1:937274/64‑1 (MQ=255)
                                                                cGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt     <  1:1164676/71‑1 (MQ=255)
                                                                 ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                         <  1:38333/50‑1 (MQ=255)
                                                                  gTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                     <  1:269596/37‑1 (MQ=255)
                                                                  gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                    <  1:941229/54‑1 (MQ=255)
                                                                  gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa   <  1:3817928/71‑1 (MQ=255)
                                                                   tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAt  <  1:3268370/71‑1 (MQ=255)
                                                                   tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAt  <  1:2972439/71‑1 (MQ=255)
                                                                   tAAGGACATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACt     <  1:1528519/68‑1 (MQ=255)
                                                                    aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa             >  1:1630602/1‑59 (MQ=255)
                                                                    aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa             >  1:1467427/1‑59 (MQ=255)
                                                                    aaGGACATGAAAGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa             >  1:2892022/1‑59 (MQ=255)
                                                                     aGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGAc      <  1:177403/65‑1 (MQ=255)
                                                                      |                                                                   
CTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAAT  >  minE/1202076‑1202213

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