Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,586 C→T 8.2% L247L (CTG→TTG)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,5860CT 8.2% 179.6 / 5.1 61L247L (CTG→TTG) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (23/33);  new base T (2/3);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGC  >  minE/1201516‑1201656
                                                                      |                                                                      
gATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGc                                                                        <  1:173178/71‑1 (MQ=255)
 aTGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGc                                                                        <  1:831789/70‑1 (MQ=255)
 aTGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCt                                                                       <  1:3514200/71‑1 (MQ=255)
    tccacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTgg                                                                     <  1:3125876/68‑1 (MQ=255)
    tccacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTgg                                                                     <  1:651696/68‑1 (MQ=255)
    tccacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTgg                                                                     <  1:587583/68‑1 (MQ=255)
    aaCACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTgg                                                                     <  1:1086417/70‑1 (MQ=255)
     ccacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGt                                                                    <  1:2011758/69‑1 (MQ=255)
     acacCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGt                                                                    <  1:36764/70‑1 (MQ=255)
      cacCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGtt                                                                   <  1:1832130/70‑1 (MQ=255)
      cacCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGtt                                                                   <  1:1876368/70‑1 (MQ=255)
      cacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGtc                                                                   <  1:1336907/70‑2 (MQ=255)
      cacCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGAACTACGCCGGATGCTGGtc                                                                   <  1:1807799/70‑2 (MQ=255)
          gAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTg                                                                      >  1:132956/1‑63 (MQ=255)
          gAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCATGCGCCGGATGCTGGTTTCCAg                                                              >  1:3455150/1‑71 (MQ=255)
           aaGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGTTGGTTTCCAGc                                                             >  1:3017506/1‑71 (MQ=255)
            aGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGTTGGTTTCCAGCg                                                            <  1:2978890/71‑1 (MQ=255)
                   cTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGTTGGTTTCCAGCGATAAAc                                                      <  1:1920546/70‑1 (MQ=255)
                       cAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCt                                                 <  1:1322180/71‑1 (MQ=255)
                        aGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCt                                                 >  1:643109/1‑70 (MQ=255)
                           cccGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAAc                                              >  1:2834599/1‑70 (MQ=255)
                             cGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGAt                                           <  1:755175/71‑1 (MQ=255)
                               cTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                         >  1:2774868/1‑71 (MQ=255)
                                tGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                         >  1:1219337/1‑70 (MQ=255)
                                tGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                         >  1:2335195/1‑70 (MQ=255)
                                 gcgcGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGa                                       >  1:1940334/1‑71 (MQ=255)
                                 gcgcGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGa                                       >  1:1500858/1‑71 (MQ=255)
                                 gcgcGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGa                                       >  1:219763/1‑71 (MQ=255)
                                   gcgTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGAt                                      <  1:3025978/70‑1 (MQ=255)
                                    cgTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                         <  1:2381226/66‑1 (MQ=255)
                                    cgTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATcc                                         <  1:2859988/66‑1 (MQ=255)
                                    cgTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATg                                     >  1:551647/1‑70 (MQ=255)
                                      tGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTaa                                               <  1:1816742/58‑1 (MQ=255)
                                       gCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATg                                 >  1:249885/1‑71 (MQ=255)
                                        ccACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTc                                                  <  1:2606962/53‑1 (MQ=255)
                                         cACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATa                                                         <  1:2657589/45‑1 (MQ=255)
                                          aCGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCa                              <  1:3554379/71‑1 (MQ=255)
                                          aCGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCa                              <  1:1351198/71‑1 (MQ=255)
                                            gtcgtGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCAtt                            <  1:1363936/71‑1 (MQ=255)
                                               gtGCGTCTGCCCTGCGCCGTATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTg                           <  1:1989439/69‑1 (MQ=255)
                                               gtGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAAt                                  >  1:1484559/1‑62 (MQ=255)
                                               gtGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTg                           <  1:2332688/69‑1 (MQ=255)
                                                 gCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTc                               >  1:2040877/1‑63 (MQ=255)
                                                 gCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTg                       >  1:1448094/1‑71 (MQ=255)
                                                  cGTCTGCCCTGCGCCGGATGTTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTgg                      >  1:1623180/1‑71 (MQ=255)
                                                      tGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTa                   >  1:1054020/1‑70 (MQ=255)
                                                        cccTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGAt                                      <  1:595398/49‑1 (MQ=255)
                                                          cTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACAt              <  1:1578951/71‑1 (MQ=255)
                                                           tGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTa                   <  1:3927092/65‑1 (MQ=255)
                                                           tGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             >  1:2954775/1‑71 (MQ=255)
                                                           tGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             >  1:3116987/1‑71 (MQ=255)
                                                           tGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             >  1:1333574/1‑71 (MQ=255)
                                                           tGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             >  1:436612/1‑71 (MQ=255)
                                                            gcgcCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGAt                                           >  1:3287486/1‑40 (MQ=255)
                                                             cgcCGGATGTTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGt            <  1:979659/70‑1 (MQ=255)
                                                             cgcCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGtt           <  1:1352262/71‑1 (MQ=255)
                                                               ccGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGttt          <  1:426458/70‑1 (MQ=255)
                                                                cGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTgc        <  1:3561109/71‑1 (MQ=255)
                                                                  gATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGa                                       >  1:1645287/1‑38 (MQ=255)
                                                                   aTGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATTAATGTCATTg                           >  1:1190443/1‑49 (MQ=255)
                                                                   aTGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCt     <  1:2210159/71‑1 (MQ=255)
                                                                    tGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             <  1:1288763/62‑1 (MQ=255)
                                                                    tGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATg             <  1:1255958/62‑1 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa                                   >  1:2232326/1‑38 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGaa                                   >  1:1050231/1‑38 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGa                                    <  1:897377/37‑1 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGtt           >  1:446935/1‑62 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGtt           >  1:2955572/1‑62 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTgg   <  1:2661900/70‑1 (MQ=255)
                                                                      cTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGc  <  1:3498042/71‑1 (MQ=255)
                                                                      |                                                                      
GATGAACACCGAAGGTGTACTGCCAGGCCCGCTGCGCGTGCCACGTCGTGCGTCTGCCCTGCGCCGGATGCTGGTTTCCAGCGATAAACTGTCTAACGATCCGATGAATGTCATTGACTGGGTAAACATGTTTGCGCTGGC  >  minE/1201516‑1201656

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