Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,846 C→T 17.5% S333S (TCC→TCT)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,8460CT17.5% 193.2 / 27.6 80S333S (TCC→TCT) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (7/7);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAG  >  minE/1201780‑1201913
                                                                  |                                                                   
cTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTg                                                                 >  1:1593093/1‑71 (MQ=255)
   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:398475/50‑3 (MQ=255)
                                                   gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTtggcg                                   <  1:3019909/50‑3 (MQ=255)
                                                   gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCaa                       >  1:3672582/1‑62 (MQ=9)
                                                   gAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTGGCGTGCTCaa                       >  1:885873/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:491725/1‑39 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                   >  1:2060467/1‑39 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                   >  1:3539557/1‑39 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTc                         <  1:2981976/49‑1 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg     <  1:1414700/69‑1 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg     <  1:5639978/69‑1 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGc    <  1:1377739/70‑1 (MQ=255)
                                                              tttCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGc    <  1:589542/70‑1 (MQ=255)
                                                                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:2082192/70‑1 (MQ=255)
                                                                tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg  >  1:871001/1‑70 (MQ=255)
                                                                tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAg  >  1:3259618/1‑70 (MQ=255)
                                                                  |                                                                   
CTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAG  >  minE/1201780‑1201913

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