Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,852 G→T 20.8% A335A (GCG→GCT sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,8520GT20.8% 214.6 / 27.9 81A335A (GCG→GCTsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (7/10);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTG  >  minE/1201783‑1201918
                                                                     |                                                                  
taCCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCg                                                                    >  1:4342796/1‑70 (MQ=255)
  cccGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCgg                                                                   >  1:2522213/1‑69 (MQ=255)
   ccGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCgg                                                                   <  1:3460034/68‑1 (MQ=255)
     gTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAgtt                                                              >  1:3047871/1‑71 (MQ=255)
      ttACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAgttg                                                             <  1:1793481/71‑1 (MQ=255)
         cTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAgttggt                                                           <  1:1758321/70‑1 (MQ=255)
            tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAg                                                                >  1:3451637/1‑62 (MQ=255)
            tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAg                                                                >  1:2194064/1‑62 (MQ=255)
            tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGc                                                        >  1:4301139/1‑70 (MQ=255)
            tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGc                                                        <  1:5313353/70‑1 (MQ=255)
              tGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCa                                                      <  1:1872239/70‑1 (MQ=255)
               ggCAGCGAGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCa                                                      <  1:5188995/69‑1 (MQ=255)
                 cAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAgtt                                                              >  1:3836782/1‑59 (MQ=255)
                   gCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCa                                                      >  1:50090/1‑65 (MQ=255)
                   gCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGa                                                <  1:1280949/71‑1 (MQ=255)
                   gCGGGCGCGATTGGTACATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGcc                                                       >  1:380683/1‑64 (MQ=255)
                        cgcgATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGa                                                <  1:5555514/66‑1 (MQ=255)
                        cgcgATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTgg                                           <  1:3280033/71‑1 (MQ=255)
                             ttGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGtt                                      >  1:3776762/1‑71 (MQ=255)
                             ttGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGtt                                      >  1:786844/1‑71 (MQ=255)
                              tGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGtt                                      >  1:3453460/1‑70 (MQ=255)
                                gTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGcc                                   <  1:5654169/71‑1 (MQ=255)
                                 tGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCt                                  <  1:2979584/71‑1 (MQ=255)
                                  gCATTGTATAAAATGAACGCATCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAg                                              >  1:2370817/1‑58 (MQ=255)
                                   cATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGt                                >  1:2230101/1‑71 (MQ=255)
                                   cATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGt                                >  1:3849571/1‑71 (MQ=255)
                                   cATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGt                                >  1:419170/1‑71 (MQ=255)
                                     ttGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGa                                                <  1:32093/53‑1 (MQ=255)
                                     ttGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTggg                                          <  1:1045454/59‑1 (MQ=255)
                                       gTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAg                                                     >  1:3311333/1‑46 (MQ=255)
                                       gTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                        >  1:2109799/1‑59 (MQ=255)
                                       gTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCaa                            >  1:4111029/1‑71 (MQ=255)
                                        tataAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTc                              <  1:2274963/68‑1 (MQ=9)
                                        tataAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCAAt                           >  1:2170275/1‑71 (MQ=255)
                                        tataAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGa                                                <  1:2739620/50‑1 (MQ=255)
                                        tataAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGa                                                <  1:5571954/50‑1 (MQ=255)
                                         ataAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGAGAAGTGGGTGTTGCCTGTTc                              <  1:1090400/67‑1 (MQ=255)
                                         aaaaaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTc                              <  1:1201554/65‑1 (MQ=255)
                                          taAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAgggtg                                                 <  1:766365/47‑3 (MQ=16)
                                           aaaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAg                                              <  1:3502844/49‑1 (MQ=255)
                                           aaaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgc                        >  1:2639675/1‑71 (MQ=255)
                                            agaTGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGaa                                               >  1:2606903/3‑47 (MQ=16)
                                            agaTGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGaa                                               >  1:4343575/3‑47 (MQ=16)
                                            aaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGc                                                        <  1:5607852/38‑1 (MQ=255)
                                            aaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCa                                                      >  1:3233926/1‑40 (MQ=255)
                                            aaaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAg                                                     >  1:2455315/1‑41 (MQ=255)
                                             aaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCa                                                      >  1:875319/1‑39 (MQ=255)
                                             aaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgctg                      >  1:2878333/1‑71 (MQ=255)
                                             aaTCAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgctg                      <  1:4324716/71‑1 (MQ=255)
                                                gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTtggcg                                        <  1:1311146/50‑3 (MQ=255)
                                                gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTtggcg                                        <  1:3019909/50‑3 (MQ=255)
                                                gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTtggcg                                        <  1:398475/50‑3 (MQ=255)
                                                gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCaa                            >  1:885873/1‑62 (MQ=9)
                                                gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCaa                            >  1:3672582/1‑62 (MQ=9)
                                                       tctATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCt            >  1:1281896/1‑71 (MQ=255)
                                                       tcgaTTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCAATGGCGGCGGCTGGTCt            <  1:3856265/68‑1 (MQ=2)
                                                        ctATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATg                          <  1:5575919/56‑1 (MQ=255)
                                                        ctATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATg                          <  1:2330361/56‑1 (MQ=255)
                                                         tATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCtt           >  1:1033113/1‑70 (MQ=255)
                                                          aTTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGtt                                            >  1:1172022/1‑35 (MQ=16)
                                                          aTTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGgcgtgg                                     <  1:3633851/43‑3 (MQ=16)
                                                          aTTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGc                                    >  1:3698550/1‑44 (MQ=9)
                                                          aTTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg          >  1:423164/1‑70 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTggg                                          >  1:4249059/1‑37 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                        >  1:2060467/1‑39 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                        >  1:3539557/1‑39 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                        >  1:491725/1‑39 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTc                              <  1:2981976/49‑1 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg          <  1:5639978/69‑1 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg          <  1:1414700/69‑1 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGc         <  1:589542/70‑1 (MQ=255)
                                                           tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGc         <  1:1377739/70‑1 (MQ=255)
                                                             tctggTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCAATGGCGGCGGCTGGTCTGGCAg       <  1:2404421/67‑1 (MQ=2)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGc                                    <  1:4917877/41‑1 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGc                                    <  1:5364911/41‑1 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       <  1:5555480/70‑1 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       >  1:1969824/1‑70 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       >  1:3259618/1‑70 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       <  1:2082192/70‑1 (MQ=255)
                                                             tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       >  1:871001/1‑70 (MQ=255)
                                                               tggTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTc                              <  1:2748762/44‑1 (MQ=9)
                                                                 gTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAActg  >  1:4693927/1‑71 (MQ=255)
                                                                  tGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCAAt                           <  1:2913148/45‑1 (MQ=255)
                                                                    cGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg       <  1:4713454/63‑1 (MQ=255)
                                                                     |                                                                  
TACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTG  >  minE/1201783‑1201918

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: