Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,358,189 G→T 10.7% S311R (AGC→AGA yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,358,1890GT10.7% 178.9 / 4.7 55S311R (AGC→AGAyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base T (3/3);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACC  >  minE/1358119‑1358258
                                                                      |                                                                     
gTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACg                                                                       >  1:1066621/1‑71 (MQ=255)
gTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACg                                                                       >  1:3946783/1‑71 (MQ=255)
    cagcaTATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGa                                                                   <  1:3572516/71‑1 (MQ=255)
      gcaTATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACa                                                                 <  1:1155151/71‑1 (MQ=255)
      gcaTATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACa                                                                 <  1:3749805/71‑1 (MQ=255)
             tGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATc                                                               >  1:197242/1‑66 (MQ=255)
             tGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGa                                                           <  1:473836/70‑1 (MQ=255)
                 aCCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGc                                                      <  1:316650/71‑1 (MQ=255)
                 aCCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGc                                                      >  1:3047392/1‑71 (MQ=255)
                  ccGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTg                                                       <  1:1746270/69‑1 (MQ=255)
                    gCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGa                                                                   >  1:581310/1‑55 (MQ=255)
                     cTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACTCTGACATCAAGATCTGCTcca                                                  >  1:752372/1‑71 (MQ=255)
                     cTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAg                                                            <  1:2560840/61‑1 (MQ=255)
                        cATAAGGCACCAGATAGGTGGTTACCGAGCCTTCCGCCTCTTTCACGCTGACATCAAGATCTGCTccac                                                 <  1:936378/69‑1 (MQ=255)
                        cATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTccacc                                                <  1:2219587/70‑1 (MQ=255)
                        cATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTccac                                                 <  1:3855713/69‑1 (MQ=255)
                            aGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATc                                                               <  1:2647077/51‑1 (MQ=255)
                            aGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGAt                                                          >  1:2438771/1‑56 (MQ=255)
                                aCCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACTCTGACATCAAGATCTg                                                       >  1:1271053/1‑55 (MQ=255)
                                 ccAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTccacc                                                <  1:1389868/61‑1 (MQ=255)
                                  cAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACg                                                                       <  1:1879226/37‑1 (MQ=255)
                                    gATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAg                                                            <  1:825813/46‑1 (MQ=255)
                                    gATAGGTGGTTACCGAGCCGTCCGCCGCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACt                                     >  1:2696352/1‑69 (MQ=255)
                                          tggtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAg                                                            <  1:776050/40‑1 (MQ=255)
                                          tggtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCt                             <  1:2836164/71‑1 (MQ=255)
                                           ggtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGc                                                      <  1:1840594/45‑1 (MQ=255)
                                            gtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTg                            >  1:2628100/1‑70 (MQ=255)
                                            gtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGt                           <  1:2277863/71‑1 (MQ=255)
                                              tACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCt                                                     >  1:1879513/1‑43 (MQ=255)
                                               aCCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGa                                                           <  1:1643702/36‑1 (MQ=255)
                                               aCCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAAt                        >  1:2425598/1‑71 (MQ=255)
                                                 cGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATc                                                         >  1:3473983/1‑36 (MQ=255)
                                                 cGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATcg                      <  1:3353314/71‑1 (MQ=255)
                                                  gAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTc                                                    <  1:168934/40‑1 (MQ=255)
                                                  gAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCaa                                 <  1:2723204/59‑1 (MQ=255)
                                                  gAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCaa                                 <  1:4020911/59‑1 (MQ=255)
                                                   aGCCGTCCGCCTCTTTCACTCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATcgcg                    <  1:179673/71‑1 (MQ=255)
                                                       gTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCaa                                 <  1:2758002/54‑1 (MQ=255)
                                                          cGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAAc                                      >  1:1879574/1‑46 (MQ=255)
                                                          cGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGa                   >  1:2653491/1‑65 (MQ=255)
                                                          cGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGa                   >  1:2330014/1‑65 (MQ=255)
                                                          cGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGa                   >  1:1476923/1‑65 (MQ=255)
                                                          cGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGc             >  1:4051702/1‑71 (MQ=255)
                                                            ccTCTTTCACTCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGc             >  1:588987/1‑69 (MQ=255)
                                                             ctctTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCag          >  1:197820/1‑71 (MQ=255)
                                                             ctctTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCag          <  1:4095044/71‑1 (MQ=255)
                                                                tttCACTCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATc                              <  1:865280/48‑1 (MQ=255)
                                                                tttCACTCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATc                              <  1:3943705/48‑1 (MQ=255)
                                                                tttCACGCTGACATCAAGATCTGCTCCACCACCGGCCaa                                       >  1:310526/1‑39 (MQ=255)
                                                                tttCACGCTGACATCAAGATCTGCTCCACCACCGGCCaa                                       >  1:2632717/1‑39 (MQ=255)
                                                                tttCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTg                            <  1:2295198/50‑1 (MQ=255)
                                                                 ttCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCaggag       >  1:478678/1‑70 (MQ=255)
                                                                 ttCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCa           >  1:3999930/1‑66 (MQ=255)
                                                                  tCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCaggagg      <  1:3967905/70‑1 (MQ=255)
                                                                  tCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCaggagg      >  1:2005272/1‑70 (MQ=255)
                                                                  tCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCaggagg      <  1:1722985/70‑1 (MQ=255)
                                                                  tCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCaggagg      <  1:103451/70‑1 (MQ=255)
                                                                     cGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAAcc  <  1:434091/71‑1 (MQ=255)
                                                                      |                                                                     
GTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACC  >  minE/1358119‑1358258

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: