Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,417,574 T→A 10.4% intergenic (+6/‑69) 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,5740TA10.4% 122.8 / 7.2 48intergenic (+6/‑69)yejM/proLpredicted hydrolase, inner membrane/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base A (3/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGT  >  minE/1417504‑1417642
                                                                      |                                                                    
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:1386511/71‑1 (MQ=255)
      cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                                <  1:2002820/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)
             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:2065877/1‑71 (MQ=255)
                                                   ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                   >  1:3307358/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:2510541/71‑1 (MQ=255)
                                                                attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCAc      >  1:277274/1‑71 (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)
                                                                    attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc        <  1:2358288/65‑1 (MQ=255)
                                                                    attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc        <  1:1427448/65‑1 (MQ=255)
                                                                      tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt  <  1:671778/69‑1 (MQ=255)
                                                                      tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt  <  1:2612675/69‑1 (MQ=255)
                                                                      |                                                                    
AAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGT  >  minE/1417504‑1417642

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: