Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,201,054 T→C 12.6% G69G (GGT→GGC)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,201,0540TC12.6% 201.0 / 17.9 79G69G (GGT→GGC) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base T (24/45);  new base C (6/4);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

CGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAG  >  minE/1200985‑1201123
                                                                     |                                                                     
cGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTc                                                                      <  1:3245087/71‑1 (MQ=255)
  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:3246124/1‑71 (MQ=255)
                      tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                >  1:1524997/1‑71 (MQ=255)
                       cGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                                                                       >  1:2456999/1‑47 (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:1698745/1‑71 (MQ=255)
                                  aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                    >  1:2807012/1‑71 (MQ=255)
                                  aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                    >  1:28292/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)
                                     ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                                     <  1:3375099/35‑1 (MQ=255)
                                     ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                                     <  1:1181678/35‑1 (MQ=255)
                                     ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                                     <  1:365091/35‑1 (MQ=255)
                                     ggCCACCACACCGATATCGCCATTATTATGGGTCt                                                                     <  1:3632156/35‑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:2613467/1‑67 (MQ=255)
                                      gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                 >  1:36065/1‑67 (MQ=255)
                                        caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAcc                                               <  1:1450100/54‑1 (MQ=255)
                                        caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                              <  1:3356469/71‑1 (MQ=255)
                                        caccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTAt                              <  1:1369986/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:765074/70‑1 (MQ=255)
                                                                atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag       <  1:2999549/70‑1 (MQ=255)
                                                                atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag       <  1:3080664/70‑1 (MQ=255)
                                                                 tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga      <  1:2166205/70‑1 (MQ=255)
                                                                 tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTaga      <  1:2338361/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)
                                                                     |                                                                     
CGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAG  >  minE/1200985‑1201123

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