Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,057 T→G 14.1% L70L (CTT→CTG)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,0570TG14.1% 141.2 / 22.3 71L70L (CTT→CTG) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base T (21/40);  new base G (6/4);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGC  >  minE/1200987‑1201126
                                                                      |                                                                     
tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                                       >  1:2874367/1‑71 (MQ=255)
tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                                       >  1:2573593/1‑71 (MQ=255)
tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                                       >  1:1398921/1‑71 (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)
      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCCTGCAg                                                                   <  1:2544885/69‑1 (MQ=255)
       ttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                                                 >  1:395702/1‑70 (MQ=255)
         aTGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                                              >  1:3449968/1‑71 (MQ=255)
          tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAAc                                                              <  1:1349149/70‑1 (MQ=255)
          tGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                             <  1:3669786/71‑1 (MQ=255)
            gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                            <  1:718191/70‑1 (MQ=255)
            gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                            <  1:2893648/70‑1 (MQ=255)
             ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                          <  1:260318/71‑1 (MQ=255)
             ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                          <  1:2552819/71‑1 (MQ=255)
              tCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                                                      <  1:2189234/58‑1 (MQ=255)
               cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                        >  1:3847740/1‑71 (MQ=255)
                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                      >  1:3577907/1‑70 (MQ=255)
                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                             <  1:3503918/61‑1 (MQ=255)
                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                   >  1:1524997/1‑71 (MQ=255)
                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                   >  1:3246124/1‑71 (MQ=255)
                      gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                                                  <  1:1479701/70‑1 (MQ=255)
                      gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                 <  1:191605/71‑1 (MQ=255)
                      gCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                                 <  1:2339332/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:28292/1‑71 (MQ=255)
                                aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                       >  1:1698745/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)
                                   ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                      >  1:2891181/1‑53 (MQ=255)
                                   ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                        <  1:693914/67‑1 (MQ=255)
                                   ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                    <  1:242644/71‑1 (MQ=255)
                                    gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                    >  1:36065/1‑67 (MQ=255)
                                    gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                    >  1:2613467/1‑67 (MQ=255)
                                      caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                                                  <  1:1450100/54‑1 (MQ=255)
                                      caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                 <  1:1369986/71‑1 (MQ=255)
                                      caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                 <  1:3356469/71‑1 (MQ=255)
                                      caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                 >  1:658597/1‑71 (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)
                                                     gCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcg                  >  1:2261783/1‑71 (MQ=255)
                                                      ccATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTcag                 <  1:923037/71‑3 (MQ=255)
                                                         ttattaTGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGAc              <  1:2480922/71‑1 (MQ=255)
                                                             tatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                <  1:910709/65‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          >  1:2038524/1‑70 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          <  1:1986121/70‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          <  1:1646315/70‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          <  1:3080664/70‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          <  1:2999549/70‑1 (MQ=255)
                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag          <  1:765074/70‑1 (MQ=255)
                                                               tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga         <  1:1195346/70‑1 (MQ=255)
                                                               tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga         <  1:2166205/70‑1 (MQ=255)
                                                               tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga         <  1:2338361/70‑1 (MQ=255)
                                                                  gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag     <  1:24227/71‑1 (MQ=255)
                                                                   tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag     <  1:2159079/70‑1 (MQ=255)
                                                                     ttGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc  >  1:1185062/1‑71 (MQ=255)
                                                                      tGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc  >  1:3569103/1‑70 (MQ=255)
                                                                      tGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc  >  1:3573672/1‑70 (MQ=255)
                                                                      |                                                                     
TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGC  >  minE/1200987‑1201126

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