Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,417,571 A→T 10.0% intergenic (+3/‑72) yejM → / → proL predicted hydrolase, inner membrane/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,417,5710AT10.0% 135.4 / 7.0 50intergenic (+3/‑72)yejM/proLpredicted hydrolase, inner membrane/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base T (3/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGA  >  minE/1417504‑1417640
                                                                   |                                                                     
aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                    <  1:967533/71‑1 (MQ=255)
  aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtat                                                                  <  1:945412/71‑1 (MQ=255)
  aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGAAAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                   >  1:3840357/1‑70 (MQ=255)
   aaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtat                                                                  >  1:1197518/1‑70 (MQ=255)
    aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                <  1:328293/71‑1 (MQ=255)
      cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                              <  1:2002820/71‑1 (MQ=255)
      cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                              <  1:1386511/71‑1 (MQ=255)
       acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                             <  1:1887043/71‑1 (MQ=255)
       acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                             <  1:3531895/71‑1 (MQ=255)
          gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                          <  1:3860134/71‑1 (MQ=255)
           ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                         <  1:901218/71‑1 (MQ=255)
             aaGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaat                                                                     <  1:2169344/57‑1 (MQ=255)
             aaGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                       >  1:1400153/1‑71 (MQ=255)
                 ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                    >  1:1425157/1‑70 (MQ=255)
                  tGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGc                                                      >  1:1748764/1‑67 (MQ=255)
                  tGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                   <  1:2827750/70‑1 (MQ=255)
                   gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                          >  1:771757/1‑62 (MQ=255)
                   gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                          >  1:3120557/1‑62 (MQ=255)
                    ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAat                                                 <  1:1048797/70‑1 (MQ=255)
                    ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAat                                                 <  1:1781982/70‑1 (MQ=255)
                      aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                               >  1:299380/1‑54 (MQ=255)
                      aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                               <  1:764674/70‑1 (MQ=255)
                      aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                              >  1:465759/1‑71 (MQ=255)
                        tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                             >  1:254282/1‑70 (MQ=255)
                           aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                         >  1:264495/1‑71 (MQ=255)
                             gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTa                                                        >  1:210634/1‑54 (MQ=21)
                               gCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                      >  1:1388463/1‑70 (MQ=255)
                                   aCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                           >  1:3713103/1‑45 (MQ=255)
                                   aCAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                 >  1:2503860/1‑71 (MQ=255)
                                    cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                <  1:1862576/55‑1 (MQ=255)
                                    cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgg                                >  1:977796/1‑71 (MQ=255)
                                    cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgg                                >  1:1854972/1‑71 (MQ=255)
                                     aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTa                                                        >  1:141979/1‑46 (MQ=255)
                                     aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                               <  1:2495695/71‑1 (MQ=255)
                                     aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                               <  1:3126211/71‑1 (MQ=255)
                                      gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaaa                                                  <  1:2755347/51‑1 (MQ=255)
                                       aCGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                             >  1:2997372/1‑71 (MQ=255)
                                           gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                    <  1:1560001/44‑1 (MQ=255)
                                           gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaa                         <  1:64277/71‑1 (MQ=255)
                                           gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGa                          >  1:275958/1‑70 (MQ=255)
                                             aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaa                       >  1:2994099/1‑71 (MQ=255)
                                              gCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                          <  1:1205660/51‑1 (MQ=255)
                                                 ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                   <  1:611956/55‑1 (MQ=255)
                                                 ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                           <  1:3149538/63‑1 (MQ=255)
                                                   ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                 >  1:3307358/1‑71 (MQ=255)
                                                   ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                 >  1:2065877/1‑71 (MQ=255)
                                                           aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCg                                 <  1:2583984/47‑1 (MQ=18)
                                                           aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCg                                 <  1:2802486/47‑1 (MQ=18)
                                                           aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTag           >  1:1951321/1‑69 (MQ=255)
                                                             cTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc       >  1:3895461/1‑71 (MQ=38)
                                                             cTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc          >  1:1046394/1‑68 (MQ=37)
                                                             cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc       <  1:3920153/71‑1 (MQ=255)
                                                             cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcag        <  1:2848811/70‑1 (MQ=255)
                                                                attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt              <  1:689402/61‑1 (MQ=255)
                                                                attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc      <  1:158244/69‑1 (MQ=255)
                                                                attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc    >  1:277274/1‑71 (MQ=255)
                                                                attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc    <  1:2510541/71‑1 (MQ=255)
                                                                 ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                 <  1:2222780/41‑1 (MQ=255)
                                                                 ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                               <  1:1779831/43‑1 (MQ=255)
                                                                 ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc    <  1:3022077/70‑1 (MQ=255)
                                                                  taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaa                       >  1:330659/1‑50 (MQ=255)
                                                                   aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa  <  1:2108266/70‑1 (MQ=255)
                                                                   |                                                                     
AAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGA  >  minE/1417504‑1417640

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: