Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,933 A→G 15.4% S440G (AGT→GGT)  yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,449,9330AG15.4% 257.2 / 36.2 104S440G (AGT→GGT) yifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base G (10/6);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 107 total aligned reads displayed.

AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGATCCA  >  minE/2449864‑2450003
                                                                     |                                                                      
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                       >  1:1320283/1‑71 (MQ=255)
aaaCAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAg                                                                       >  1:156309/1‑71 (MQ=255)
 aacaacGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGt                                                                      <  1:2049611/71‑1 (MQ=255)
    aacGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCg                                                                   >  1:2868102/1‑71 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCgg                                                                  >  1:2768149/1‑70 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCgg                                                                  >  1:3299594/1‑70 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGa                                                                 <  1:3639391/71‑1 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGa                                                                 <  1:2499677/71‑1 (MQ=255)
      cgcgTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGa                                                                 <  1:2243007/71‑1 (MQ=255)
            aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACg                                                           >  1:2536499/1‑71 (MQ=255)
            aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACg                                                           >  1:1605259/1‑71 (MQ=255)
            aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACg                                                           >  1:3176972/1‑71 (MQ=255)
            aTGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACg                                                           >  1:2844170/1‑71 (MQ=255)
             tGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAAt                                                              <  1:1114141/67‑1 (MQ=255)
              ggCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGc                                                          >  1:3129889/1‑70 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGcc                                                         >  1:815965/1‑70 (MQ=255)
               gccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGACAGTCCGGAAATACGcc                                                         >  1:1886214/1‑70 (MQ=255)
                ccgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCaa                                                       <  1:2825200/71‑1 (MQ=255)
                 cgccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAt                                                      <  1:1801135/71‑1 (MQ=255)
                  gccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGa                                                                 <  1:3336662/59‑1 (MQ=255)
                   ccAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaa                                                                <  1:3884862/59‑1 (MQ=255)
                     agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                     >  1:2839224/1‑68 (MQ=255)
                     agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAAtt                                                     >  1:3937567/1‑68 (MQ=255)
                     agagagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTg                                                  <  1:896620/71‑1 (MQ=255)
                        gagCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTg                                                  >  1:844821/1‑68 (MQ=255)
                         agCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCaaa                                              >  1:1370017/1‑71 (MQ=255)
                         agCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCaaa                                              >  1:1266917/1‑71 (MQ=255)
                           cgcCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACt                                            >  1:2814344/1‑71 (MQ=255)
                            gcCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           <  1:2773648/71‑1 (MQ=255)
                             ccAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaaa                                                               <  1:1648160/50‑1 (MQ=255)
                             ccAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTg                                          <  1:1018421/71‑1 (MQ=255)
                              cAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           <  1:1745012/69‑1 (MQ=255)
                              cAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           >  1:3793019/1‑69 (MQ=255)
                              cAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           >  1:3788143/1‑69 (MQ=255)
                                aTTTAAGCGGCAGAATAAACTGAA‑‑‑‑TGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                    >  1:1323374/1‑70 (MQ=255)
                                 tttAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGaaa                                                               <  1:2847142/46‑1 (MQ=255)
                                 tttAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGagag                                      >  1:3373582/1‑71 (MQ=255)
                                   tAAGCGGCAGAATAAACTGAATGCCTGGCTGGATa                                                                        >  1:2835218/1‑35 (MQ=255)
                                   tAAGCGGCAGAATAAACTGAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTga                                         >  1:2507420/1‑66 (MQ=255)
                                     aGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                                  <  1:3793194/71‑1 (MQ=255)
                                      gCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                    >  1:35075/1‑68 (MQ=255)
                                      gCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAg                                 <  1:678320/71‑1 (MQ=255)
                                       cGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                                <  1:3301652/71‑1 (MQ=255)
                                        ggCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:2094757/71‑1 (MQ=255)
                                        ggCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:535652/71‑1 (MQ=255)
                                        ggCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:56450/71‑1 (MQ=255)
                                        ggCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:1240123/71‑1 (MQ=255)
                                         gCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           <  1:3742470/58‑1 (MQ=255)
                                         gCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:1337716/70‑1 (MQ=255)
                                          cAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                             <  1:2633790/71‑1 (MQ=255)
                                          cAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGc                             <  1:1711824/71‑1 (MQ=255)
                                            gAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGa                           >  1:2187535/1‑71 (MQ=255)
                                             aaTAAACTGAATGCCTGGCTGGATAGTCCGGAAAt                                                              <  1:543577/35‑1 (MQ=255)
                                              aTAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtg                         >  1:1868069/1‑71 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                         >  1:2245683/1‑54 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                         >  1:3038312/1‑54 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                        >  1:3608977/1‑71 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                        >  1:395288/1‑71 (MQ=255)
                                               tAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgt                        >  1:463356/1‑71 (MQ=255)
                                                aaaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAtgtg                       >  1:1302766/1‑71 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACtt                                           <  1:1344115/50‑1 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                    <  1:2929964/57‑1 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:682569/62‑1 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAt                          >  1:3643311/1‑67 (MQ=255)
                                                 aaCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTgg                      >  1:1147316/1‑71 (MQ=255)
                                                  aCTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                         >  1:4062354/1‑51 (MQ=255)
                                                   cTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGtt                    >  1:136739/1‑71 (MQ=255)
                                                    tGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                                  <  1:1286044/56‑1 (MQ=255)
                                                    tGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                                  <  1:1241316/56‑1 (MQ=255)
                                                     gAATGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   >  1:939182/1‑70 (MQ=255)
                                                     gAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   >  1:347627/1‑70 (MQ=255)
                                                     gAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   >  1:1230715/1‑70 (MQ=255)
                                                     gAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   >  1:2743023/1‑70 (MQ=255)
                                                     gAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   >  1:432959/1‑70 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg                  <  1:632241/70‑1 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg                  >  1:3330705/1‑70 (MQ=255)
                                                      aaTGCCTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTg                   <  1:3678144/69‑1 (MQ=255)
                                                        tGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaag               >  1:3088143/1‑71 (MQ=255)
                                                        tGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaag               >  1:543631/1‑71 (MQ=255)
                                                        tGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTgg                  <  1:650273/68‑1 (MQ=255)
                                                          ccTGGCTGGATGGTCCGGAAATACGCCAATTCTGCAAAc                                             >  1:1175857/1‑39 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGTAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg              <  1:3594535/70‑1 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg              >  1:137892/1‑70 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg              >  1:272299/1‑70 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg              >  1:1185618/1‑70 (MQ=255)
                                                          ccTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaagg              >  1:3396066/1‑70 (MQ=255)
                                                           ctggctggATGGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTggaaggaa            <  1:1758075/71‑1 (MQ=255)
                                                              gctggATAGTCCGGAAATACGCCAATTCTGCAAACTTga                                         >  1:1350333/1‑39 (MQ=255)
                                                              gctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                                >  1:3444286/1‑48 (MQ=255)
                                                              gctggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAaca          >  1:2430758/1‑70 (MQ=255)
                                                                 ggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                    <  1:532320/41‑1 (MQ=255)
                                                                 ggATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCg                                    <  1:835333/41‑1 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGa                                >  1:1007322/1‑44 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTg       <  1:424856/69‑1 (MQ=255)
                                                                  gATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTg       <  1:1411487/69‑1 (MQ=255)
                                                                   aTAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGAt                               <  1:422952/44‑1 (MQ=255)
                                                                   aTAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGAt                          >  1:3692367/1‑49 (MQ=255)
                                                                   aTAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGATc    >  1:2358853/1‑71 (MQ=255)
                                                                     aGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGaa                                  >  1:1551162/1‑39 (MQ=255)
                                                                     aGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGATCCa  >  1:3456261/1‑71 (MQ=255)
                                                                     |                                                                      
AAACAACGCGTAATGGCCGCCAGAGAGCGCCAATTTAAGCGGCAGAATAAACTGAATGCCTGGCTGGATAGTCCGGAAATACGCCAATTCTGCAAACTTGAGAGCGAAGATGCGATGTGGTTGGAAGGAACACTGATCCA  >  minE/2449864‑2450003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: