Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,074 C→A 11.2% A76D (GCC→GAC) ‡ sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,0740CA11.2% 248.1 / 11.0 80A76D (GCC→GAC) ‡sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (29/42);  new base A (5/4);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCT  >  minE/1201005‑1201140
                                                                     |                                                                  
tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                                    >  1:3577907/1‑70 (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: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)
             tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATc                                                       >  1:2437087/1‑70 (MQ=255)
              aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                                      >  1:3160451/1‑70 (MQ=255)
              aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                     >  1:28292/1‑71 (MQ=255)
              aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                     >  1:1698745/1‑71 (MQ=255)
              aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                     >  1:2807012/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)
                aGGCCACCACACCGATATCGCCATATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                                  <  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:658597/1‑71 (MQ=255)
                    caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                               <  1:1369986/71‑1 (MQ=255)
                    caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                                               <  1:3356469/71‑1 (MQ=255)
                      ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                                      <  1:1035295/62‑1 (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:765074/70‑1 (MQ=255)
                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                        <  1:2999549/70‑1 (MQ=255)
                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                        <  1:3080664/70‑1 (MQ=255)
                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                        <  1:1646315/70‑1 (MQ=255)
                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                        <  1:1986121/70‑1 (MQ=255)
                                            atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag                        >  1:2038524/1‑70 (MQ=255)
                                             tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga                       <  1:2338361/70‑1 (MQ=255)
                                             tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga                       <  1:2166205/70‑1 (MQ=255)
                                             tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga                       <  1:1195346/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:3573672/1‑70 (MQ=255)
                                                    tGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                >  1:3569103/1‑70 (MQ=255)
                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                               >  1:3881240/1‑54 (MQ=255)
                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                <  1:1174304/69‑1 (MQ=255)
                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg               <  1:1485302/70‑1 (MQ=255)
                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGcgcg               <  1:681754/70‑1 (MQ=255)
                                                     gCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATT‑‑CGACGTAGAAGAGCGCGAAc            >  1:2948877/1‑71 (MQ=255)
                                                      cAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg                           >  1:3610725/1‑57 (MQ=255)
                                                      cAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGa              <  1:453335/70‑1 (MQ=255)
                                                       aGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAAc            >  1:2086589/1‑71 (MQ=255)
                                                        ggTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACg           <  1:2962713/71‑1 (MQ=255)
                                                         gTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                      <  1:288050/59‑1 (MQ=255)
                                                          tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                      <  1:2480684/58‑1 (MQ=255)
                                                          tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa                      <  1:66606/58‑1 (MQ=255)
                                                          tAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgcgc                <  1:2145602/64‑1 (MQ=255)
                                                           aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                                       >  1:1307446/1‑40 (MQ=255)
                                                           aaCGAACCTGCCACCGTGGATATCGACAGTATTCCCGGtt                                       >  1:2614726/1‑40 (MQ=255)
                                                             cGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGt          >  1:3659357/1‑67 (MQ=255)
                                                             cGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgc      >  1:1764377/1‑71 (MQ=255)
                                                             cGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgc      >  1:147160/1‑71 (MQ=255)
                                                              gAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgct     <  1:2278587/71‑1 (MQ=255)
                                                              gAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgct     <  1:697878/71‑1 (MQ=255)
                                                              gAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgct     <  1:1674090/71‑1 (MQ=255)
                                                                aCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                                   <  1:99089/39‑1 (MQ=255)
                                                                aCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgc   >  1:50924/1‑71 (MQ=255)
                                                                aCCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgc   <  1:3620426/71‑1 (MQ=255)
                                                                 ccTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgct  >  1:622907/1‑71 (MQ=255)
                                                                     |                                                                  
TGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCT  >  minE/1201005‑1201140

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: