Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,033 C→T 10.1% T62T (ACC→ACT)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,0330CT10.1% 280.1 / 17.4 99T62T (ACC→ACT) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (42/47);  new base T (6/4);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTT  >  minE/1200965‑1201103
                                                                    |                                                                      
gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                                      >  1:1042811/1‑71 (MQ=255)
gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                                      <  1:3646270/71‑1 (MQ=255)
gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                                      >  1:2029965/1‑71 (MQ=255)
gATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGa                                                                      >  1:3287702/1‑71 (MQ=255)
 aTAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                                                     <  1:228341/71‑1 (MQ=255)
  tAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGata                                                                    <  1:604400/71‑1 (MQ=255)
  tAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGata                                                                    >  1:2492676/1‑71 (MQ=255)
  tAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGata                                                                    <  1:3548074/71‑1 (MQ=255)
  tAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                                                     >  1:778121/1‑70 (MQ=255)
      gTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATc                                                                  <  1:2093260/69‑1 (MQ=255)
       ttACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGcc                                                               <  1:3484305/71‑1 (MQ=255)
       ttACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGcc                                                               <  1:3535816/71‑1 (MQ=255)
        tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCa                                                              >  1:3162831/1‑71 (MQ=255)
        tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCa                                                              >  1:3120536/1‑71 (MQ=255)
        tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCa                                                              >  1:648719/1‑71 (MQ=255)
         aCTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                                             >  1:1257673/1‑71 (MQ=255)
            cgcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCg                                                                       >  1:2423655/1‑58 (MQ=255)
            cgcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCg                                                                       >  1:2053472/1‑58 (MQ=255)
            cgcgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                                                           >  1:85179/1‑70 (MQ=255)
              cgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCTCCAttatta                                                        <  1:24625/71‑1 (MQ=255)
              cgTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCAttatta                                                        <  1:1549943/71‑1 (MQ=255)
               gTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATtattat                                                       <  1:1685925/71‑1 (MQ=255)
                  gCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCAttatta                                                        >  1:1931818/1‑67 (MQ=255)
                    cGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTc                                                  <  1:3245087/71‑1 (MQ=255)
                     gTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGat                                                                     <  1:765252/51‑1 (MQ=255)
                      ttgACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCattat                                                          <  1:407611/59‑1 (MQ=255)
                      tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                                                           >  1:2877505/1‑60 (MQ=255)
                      tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                >  1:2874367/1‑71 (MQ=255)
                      tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                >  1:2573593/1‑71 (MQ=255)
                      tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                >  1:1398921/1‑71 (MQ=255)
                       ggACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATg                                                      <  1:1942937/64‑1 (MQ=255)
                        gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                <  1:2700734/69‑1 (MQ=255)
                        gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                               <  1:2056961/70‑1 (MQ=255)
                        gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                              <  1:3312334/71‑1 (MQ=255)
                          cGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGcgg                                            >  1:1093698/1‑69 (MQ=255)
                           gTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                                                           <  1:21899/55‑1 (MQ=255)
                            tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCCTGCAg                                            <  1:2544885/69‑1 (MQ=255)
                             ttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                          >  1:395702/1‑70 (MQ=255)
                               aTGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                       >  1:3449968/1‑71 (MQ=255)
                                tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                                                           >  1:366244/1‑50 (MQ=255)
                                tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatta                                                           >  1:353146/1‑50 (MQ=255)
                                tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                       <  1:1349149/70‑1 (MQ=255)
                                tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                      <  1:3669786/71‑1 (MQ=255)
                                  gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGc                                                                <  1:3219631/43‑1 (MQ=255)
                                  gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                     <  1:2893648/70‑1 (MQ=255)
                                  gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                     <  1:718191/70‑1 (MQ=255)
                                   ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                   <  1:2552819/71‑1 (MQ=255)
                                   ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                   <  1:260318/71‑1 (MQ=255)
                                    tCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                               <  1:2189234/58‑1 (MQ=255)
                                     cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                                             >  1:953794/1‑43 (MQ=255)
                                     cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCat                                                             >  1:1075389/1‑43 (MQ=255)
                                     cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                 >  1:3847740/1‑71 (MQ=255)
                                        tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                               >  1:3577907/1‑70 (MQ=255)
                                          tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCatt                                                            >  1:447872/1‑39 (MQ=255)
                                          tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                      <  1:3503918/61‑1 (MQ=255)
                                          tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                            >  1:1524997/1‑71 (MQ=255)
                                          tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                            >  1:3246124/1‑71 (MQ=255)
                                           cGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                                                   >  1:2456999/1‑47 (MQ=255)
                                            gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                           <  1:1479701/70‑1 (MQ=255)
                                            gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                          <  1:2339332/71‑1 (MQ=255)
                                            gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                          <  1:191605/71‑1 (MQ=255)
                                                 ccggTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGat                     <  1:2709983/69‑1 (MQ=255)
                                                 cGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGat                     <  1:2131114/71‑1 (MQ=255)
                                                   ggTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                   >  1:1541547/1‑71 (MQ=255)
                                                    gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                   >  1:3463086/1‑70 (MQ=255)
                                                    gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                  >  1:3458611/1‑71 (MQ=255)
                                                    gTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                   >  1:3475882/1‑70 (MQ=255)
                                                     tAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCg                 <  1:3424658/71‑1 (MQ=255)
                                                     tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                 <  1:1739371/55‑1 (MQ=255)
                                                     tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATc                  >  1:2437087/1‑70 (MQ=255)
                                                      aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                 >  1:3160451/1‑70 (MQ=255)
                                                      aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                >  1:2807012/1‑71 (MQ=255)
                                                      aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                >  1:1698745/1‑71 (MQ=255)
                                                      aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                >  1:28292/1‑71 (MQ=255)
                                                       aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                          <  1:2795167/60‑1 (MQ=255)
                                                       aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                        >  1:221149/1‑62 (MQ=255)
                                                       aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                >  1:2524524/1‑70 (MQ=255)
                                                       aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc               <  1:1115253/71‑1 (MQ=255)
                                                       aaGGCCACCACACCGATATCCCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                          <  1:1499202/60‑1 (MQ=255)
                                                        aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt              <  1:2122150/71‑3 (MQ=255)
                                                        aGGCCACCACACCGATATCGCCA‑TATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg             <  1:1121595/71‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                 <  1:3375099/35‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                 <  1:1181678/35‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                 <  1:365091/35‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                 <  1:3632156/35‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                               >  1:2891181/1‑53 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                 <  1:693914/67‑1 (MQ=255)
                                                         ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg             <  1:242644/71‑1 (MQ=255)
                                                          gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc             >  1:2613467/1‑67 (MQ=255)
                                                          gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc             >  1:36065/1‑67 (MQ=255)
                                                            caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                           <  1:1450100/54‑1 (MQ=255)
                                                            caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt          >  1:658597/1‑71 (MQ=255)
                                                            caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt          <  1:3356469/71‑1 (MQ=255)
                                                            caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt          <  1:1369986/71‑1 (MQ=255)
                                                              ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                 <  1:1035295/62‑1 (MQ=255)
                                                                acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                   >  1:850225/1‑42 (MQ=255)
                                                                 cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa           >  1:2320488/1‑65 (MQ=255)
                                                                 cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAtt         <  1:762724/67‑1 (MQ=255)
                                                                 cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATAGACAGt            >  1:1564110/1‑64 (MQ=255)
                                                                    cGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt  <  1:2998993/71‑1 (MQ=255)
                                                                    |                                                                      
GATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTT  >  minE/1200965‑1201103

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